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Abstract 

A  direct  analytical  beam  formulation  is  developed  for  predicting  the  effective  elastic  stiffnesses  and 
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generally  predicting  beam  deformations  within  10%  of  detailed  finite  element  solutions.  The  importance  of 
three  non-dassical  structural  phenomenon  is  systematically  investigated  for  coupled  composite  beams. 
Torsion-related  out-of-plane  warping  can  substantially  influence  torsion  and  coupled  torsion  deformations; 
twist  of  a  symmetric  layup  box-beam  under  a  tip  bending  load  can  increase  up  to  200%  due  to  warping. 
Couplings  associated  with  transverse  shear  deformations  can  significantly  alter  the  elastic  response  of 
tailored  composite  box-beams;  effective  bending  stiffness  of  highly  coupled  anti-symmetric  layup  beams  can 
be  reduced  more  than  30%.  Two-dimensional  in-piane  elastic  behavior  of  the  plies  is  also  very  important 
to  the  accuracy  of  composite  box-beam  analysis;  load  deflection  results  for  anti-symmetric  layup  beams  can 
be  altered  by  30-100%  by  not  accounting  for  this  elastic  behavior. 


NOTE:  This  abstract  is  being  substituted  for  the  paper  by  R.  Barrett.  It  will  be  presented  Tuesday,  March 
13, 1990  from  4:00-4:30  p.m.  by  Edward  C.  Smith. 
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Abstract 

for 

ARO/Duke  Workshop  on  Dynamics  and  Aeroelastic  Stability 
Modeling  of  Rotorcraft  Systems,  March  12-14,  1990 

The  establishment  of  the  US  Army  Aviation  Research  and 
Technology  Activity,  the  ARO  Rotorcraft  Centers  of  Excellence 
Program,  the  NASA  Rotorcraft  Program,  and  the  expansion  of 
rotorcraft  applications  in  the  military  and  civil  fields  have  all 
contributed  to  tremendous  expansion  in  rotorcraft  basic  research 
over  the  last  twenty  years.  In  view  of  the  challenging  technical 
complexity  of  rotorcraft,  it  is  not  surprising  that  many  areas  of 
research  still  need  considerable  attention.  Nevertheless,  some 
topics  inevitably  receive  substantial  attention,  while  other 
topics  could  benefit  from  more  emphasis.  The  purpose  of  this 
presentation  is  to  present  observations  gained  from  the 
government  research  perspective  with  the  hope  of  stimulating  more 
attention  for  selected  topics  and  suggesting  alternative  methods 
or  approaches  to  achieving  the  end  goals  of  rotorcraft  research. 
The  goals  themselves  warrant  some  discussion  as  well. 

The  presentation  will  include  specific  examples  to  address 
these  topics  in  the  context  of  rotorcraft  aeroelastic  stability, 
rotor  loads,  experimental  investigations  and  analysis  validation, 
and  the  role  of  comprehensive  analysis  in  rotorcraft  R&D. 
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HELICOPTER  ROTOR  DYNAMICS  AND  AEROELASTICITY: 
SOME  KEY  IDEAS  AND  INSIGHTS 


1'erelz  P.  Pricdmann* 

Mechanical,  Aerospace  and  Nuclear  Engineering  Department 
University  of  California,  Los  Angeles,  CA  90024,  U.SA. 


Abstract  -  The  purpose  of  this  paper  is  to  provide  a  detailed  discussion  of  four  important  current 
topics  in  helicopter  rotor  dynamics  and  acroelasticity.  These  topics  arc:  (1)  the  role  of 
geometric  nonlinearities  in  rotary-wing  acroelasticity;  (2)  structural  modeling,  free  vibration  and 
aeroelastic  analysis  of  composite  rotor  blades;  (3)  modeling  of  coupled  rotor/fuselagc 
aeromcchanical  problems  and  their  active  control;  and  (4)  use  of  higher  harmonic  control  (HHC) 
for  vibration  reduction  in  helicopter  rotors  in  forward  (light.  Hopefully  this  discussion  will  pro¬ 
vide  an  improved  fundamental  understanding  of  the  current  state  of  the  art.  Titus  future  research 
can  be  focused  on  problems  which  remain  to  be  solved  instead  of  producing  marginal  improve¬ 
ments  on  problems  which  arc  already  understood. 
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NOTATION 

=  two  dimensional  lift  curve  slope 
=  semi-chord 

=  first  order  system  and  control  matrices 
=  profile  drag  coefficient 

=  unit  vector  associated  with  the  undeformed  blade,  Fig.  1. 

=  triad  ex,  cy,  e2  after  deformation,  Fig.  1 
=  blade  root  offset 

=  unit  vector  in  the  direction  x,y,  and  z  respectively.  Fig.  1 

=  cost  functional 
=  flap,  lag  and  torsion  springs 
=  length  of  flexible  portion  of  blade 
=  nondimensional  fuselage  mass 
=  blade  radius 

=  transformation  matrix  between  deformed  and  undeformed  triads  of  unit  vectors,  Eq.  (1) 

=  elements  of  matrix  [S] 

=  time 

=  n  x  m  HHC  transfer  matrix 

=  components  of  the  displacement  of  a  point  on  the  elastic  axis  of  the  blade.  Fig.  1 
=  diagonal  weighting  matrix  on  vibrations 
=  diagonal  weighting  matrix  on  control  amplitudes 
=  diagonal  weighting  on  control  rate  of  change 
=  offset  between  blade  aerodynamic  center  (A.C)  and  clastic  axis  (E,A) 

=  offset  between  blade  center  of  mass  (C.G)  and  elastic  axis 
=  coordinates  shown  in  Fig.  1 
=  system  state  and  control  vectors 
=  n  x  1  vector  of  vibration  amplitudes 
=  n  x  1  vector  of  baseline  vibrations 
=  preconc  angle 

=  basis  for  order  of  magnitude,  associated  with  typical  elastic  blade  slope 
=  HHC  control  angle 
=  0(i)  — 00-1) 

=  total  pilch  angle 
=  collective  pitch 
-  cyclic  pitch  components 
=  pitch  of  the  kth  blade 

=  amplitudes  of  HHC  sine  input  in  collective,  longitudinal,  and  lateral  control  degrees  of  freedom 
=  amplitudes  of  HHC  cosine  input  in  collective,  longitudinal,  and  lateral  control  degrees  of  free¬ 
dom 

=  inflow  ratio,  and  its  cyclic  components 
=  advance  ratio 
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=  roiaiion  of  cross  scciion  of  blade  around  the  elastic  axis 
=  azimuth  angle  of  blade,  y  =  fit 
=  speed  of  rotation 
=  HHC  frequency 

=  inplane,  torsional  and  llapwisc  nondimensional  fundamental  frequency 
=  nondimensionalized.  quantity,  frequencies  by  fi,  length  by  / 

=  t>/0q/ 


1.  INTRODUCTION 

During  the  last  two  decades,  since  Locwy’s  11]  comprehensive  review  of  rotary-wing  dynamic  and  acroelastic 
problems  was  published  a  vast  body  of  published  research,  on  these  topics  has  appeared  in  the  literature.  This  substan 
tial  body  of  research  has  also  been  discussed  in  a  considerable  number  of  survey  papers  which  have  emphasized  vari¬ 
ous  aspects  of  the  rotary-wing  acroelastic  stability  and  response  problem  together  with  the  related  vibration  problems 
and  their  control  by  active  and  passive  means.  For  convenience  these  survey  papers  [1-13]  are  cited  in  chronological 
order.  Loewy’s  survey  [I]  was  followed  by  a  more  restrictive  survey  by  Dat  [2]  which  discussed  unsteady  aero¬ 
dynamic  and  vibration  problems  in  forward  flight.  Hohenemscr  [3]  discussed  some  aeromcchanical  stability  problems 
as  part  of  the  broader  field  of  flight  mechanics.  Friedmann  [4,5]  presented  a  detailed  chronological  discussion  of 
rotary-wing  aeroelasticity  emphasizing  the  role  of  geometrical  nonlmearidcs,  due  to  moderate  blade  deflections, 
unsteady  aerodynamics  and  forward  flight.  A  similar  discussion  which  was  restricted  to  the  case  of  hover  and  hinge- 
less  and  bearingless  rotors  was  also  provided  by  Ormiston  [6].  In  addition  to  these  papers  which  have  emphasized  pri¬ 
marily  acroelastic  stability,  two  other  surveys  [7,8]  have  dealt  exclusively  with  the  vibration  problem  and  its  active  and 
passive  control  in  rotorcraft. 

More  recently  Johnson  [9,10]  has  published  a  comprehensive  review  paper  which  described  both  the  acroelastic 
stability  and  rotorcraft  vibration  problems  in  the  context  of  dynamics  of  advanced  rotor  systems.  Friedmann  [1 1]  has 
described  the  main  developments  between  1983-87  emphasizing  new  methods  for  formulating  acroelastic  problems, 
treatment  of  the  forward  flight  problem,  coupled  rotor/fuselage  analyses,  structural  modeling  and  structural  optimiza¬ 
tion,  and  the  use  of  active  controls  for  vibration  reduction  and  stability  augmentation.  A  more  restrictive  review 
emphasizing  some  practical  design  aspects  capable  of  alleviating  aeromechanical  problems  was  presented  by  Miao 
( 12]  Finally  it  should  be  noted  that  a  very  comprehensive  research  report  [13]  has  been  published  recently  which  con¬ 
tains  a  detailed  review  of  research  carried  out  under  Army/NASA  sponsorship,  between  1967  and  1987. 

The  purpose  of  this  paper  is  not  to  present  another  review  of  the  literature  in  this  field.  Instead  the  main  objec¬ 
tive  here  is  to  discuss  four  specific  and  important  topics  in  rotor  dynamics  and  aeroelasticity  where  the  body  of  avail¬ 
able  research  has  reached  a  sufficient  level  of  maturity  for  keys  ideas  to  emerge.  The  emergence  of  these  key  concepts 
provides  good  insight  on  the  course  of  research  which  needs  to  be  conducted  on  these  topics  in  the  future. 

The  four  topics  which  will  be  discussed  in  this  paper  are  briefly  described  below: 

1  The  role  of  geometric  nonlinearities  in  rotary-wing  aeroelasticity.  Methods  for  the  effective  formulation  of  the 
equations  governing  the  rotary-wing  aeroelastic  problem  arc  described  and  the  importance  of  third  and  higher 
order  geometrically  nonlinear  terms  is  discussed. 

2  Structural  modeling,  free  vibration  and  aeroelastic  analysis  of  composite  rotor  blades.  Available  structural  and 
structural  dynamic  models  for  composite  rotor  blades  are  discussed  including  geometrically  nonlinear  terms 
and  general  cross-secuonal  geometries. 

3  Modeling  of  coupled  rotorlfuselage  aeromechanical  problems  and  their  active  control.  Various  approaches  for 
formulating  such  equations  in  forward  flight  are  discussed.  Various  control  approaches  for  stabilizing  air  reso¬ 
nance  in  hover  and  forward  flight  are  considered.  The  active  control  of  die  flap-lag  problem  and  coupled  flap- 
lag-torsional  problem  are  also  briefly  discussed. 

4  Use  of  higher  harmonic  control  for  vibration  reduction  in  helicopter  rotors  in  forward  flight.  Available  active 
control  techniques  for  vibration  reduction  in  forward  flight  are  discussed.  A  new  aeroelastic  simulation  capa¬ 
bility  for  higher  harmonic  control  is  presented  together  with  some  important  conclusions  obtained  from  this 
research. 


It  is  expected  that  the  discussion  of  the  four  topics  described  above  will  provide  an  improved  understanding  of 
the  state  of  the  art.  Thus  future  research  can  be  focused  on  problems  which  remain  to  be  solved. 


2.  THE  ROLE  OF  GEOMETRIC  NONLINEARITIES  IN  ROTARY-WING  AEROELASTICITY 
During  the  last  twenty  years  it  has  been  firmly  established  that  geometrical  nonlinearities,  due  to  moderate 
^  deflections,  play  an  important  role  in  the  aeroelastic  stability  analysis  of  hingeless  and  bearingless  rotor  blades 
(4-6,9- 11,13]  The  role  of  these  terms  for  articulated  blades  is  slightly  less  important  but  still  significant  [1,2,12-13], 


_  As  described  in  Refs.  5,  9  11,  and  13  a  number  of  beam  theories  accounting  for  moderate  as  well  as  large 
uctiections  of  blades  have  been  developed  and  are  available.  These  theories  can  be  divided  into  two  basic  categories: 
a)  those  which  arc  based  on  ordering  schemes  and  arc  valid  for  moderate  deflections,  and  (b)  those  which  are  valid  for 
large  deflections,  and  arc  not  based  on  ordering  schemes.  The  moderate  deflection  theories  have  been  developed  pri¬ 
marily  between  1970-80  and  the  large  deflection  theories  have  been  developed  since  1980.  It  is  therefore  quite 
relevant  to  try  and  determine  the  significance  of  higher  order  geometrical  nonlinearities,  such  as  third  order  and  higher 
order  nonlinear  terms,  on  the  acroelastic  stability  of  hingclcss  and  bcannglcss  rotor  blades  in  hover  and  forward  (light. 


T  he  source  and  structure  of  the  geometrically  nonlinear  terms  associated  with  moderate  blade  deflections  is 
conveniently  illustrated  by  the  transformation  between  the  triad  of  unit  vectors  describing  the  dcfoimed  and  unde- 
formed  state  of  a  hingclcss  rotor  blade,  sec  Figs,  la  and  lb.  Such  a  transformation,  based  on  the  assumption  of  small 
strains  and  finite  rotations  (slopes)  has  the  following  mathematical  form  [5] 
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uherc  ilie  elements  of  the  transformation  matrix  |S|  determine  the  accuracy  or  order  of  tl>;  theory.  I  or  example  for  a 
second  order  type  theory  S,|  ‘s  arc  given  by: 

S,|  =  1;  S,2=v„;  S|3=w.x 

Si)  =  ~(v,x +ou'.»);  S22  =  i;  S23  =  0 

S31  =-(w,xx-<Kx);  S32  =-«t  +  v.xw,x):  S33  =  l 

Transfomrations  of  this  type  combined  with  the  Euler  Bernoulli  assumption  have  been  used  as  the  basis  for 
moderate  deflection  beam  theories  which  are  suitable  for  the  aeroelastic  stability  and  response  analysis  of  isotropic 
hmgelcss  and  bcaringlcss  rotor  blades  (5,1 1,13]. 

When  such  transformations  are  incorporated  in  the  derivation  of  the  inertia  and  aerodynamic  operators  associ¬ 
ated  with  the  rotary  wing  aeroelastic  problem  many,  relatively  small,  nonlinear  terms  emerge.  Such  terms  clearly 
represent  considerable  complication  in  the  equations  of  motion  from  an  algebraic  point  of  view.  The  final  equations  of 
motion  can  be  obtained  using  two  different  approaches  (11). 

One  approach  is  based  on  generating  the  equations  of  motion  in  explicit  form.  The  explicit  approach  is  fre¬ 
quently  combined  with  an  ordering  scheme  (5]  to  "manage"  the  large  number  of  terms  associated  with  geometric  non- 
linearities  The  purpose  of  the  ordering  scheme  is  to  provide  a  rational  basis  for  neglecting  higher  order  terms  in  a  con¬ 
sistent  manner  A  suitable  example  for  an  ordering  scheme  used  in  deriving  the  equations  of  motion  for  a  hingeless 
rotor  blade  in  forward  flight  is  given  below.  Orders  of  magnitude  are  assigned  to  the  various  parameters  of  the  prob¬ 
lem  in  terms  of  elastic  blade  slopes  which  arc  assumed  to  be  moderate,  i.e.,  slopes  arc  of  order  E,  with  0.10  <  C  £  0.20, 
and  the  assumption 


0(1)  +  0(£2)  =  0(1) 


(3) 


is  used  when  deriving  die  equations. 

For  the  coupled  flap-lag-torsional  problem  in  forward  flight  the  orders  of  magnitude  for  the  relevant  quantities 
arc  listed  below: 

w,»  =  v,x  =  $  =  0(e);  X  =  =  Xlc  =  -^  =  =  0(e) 

±  =  |.  =  Pp=0(£):  pi  =  0(1);  e  =  0„  =  0lc=O{£*)  (4) 

!=x,=xa  =  0(e 2>;  C<j0/,  =  0(e3/2);  7  =  -§  =  ^=°(» 

It  is  important  to  note  that  ordering  schemes  arc  not  unique,  a  somewhat  different  combination  of  parameters 

(14)  can  be  also  used  as  the  basis  of  an  ordering  scheme  for  the  same  coupled  flap-lag-torsional  problem  of  a  hmgelcss 
rotor  blade  in  forward  flight.  Furthermore  ordering  schemes  arc  based  on  common  sense  and  experience  with  practical 
blade  configuration  thus  their  application  requires  both  care  and  a  certain  degree  of  flexibility. 

Ordering  schemes  are  also  convenient  when  using  general  purpose  algebraic  manipulators,  such  as  MACSYMA 

( 15)  to  derive  equations  of  motion  in  explicit  form.  Since  algebraic  tasks  are  relegated  to  a  computer,  it  is  fairly  easy 
to  retain  additional  higher  order  terms.  For  example  by  changing  the  basic  assumption  in  the  ordering  scheme  to 

0(1)  +  0(e3)  =  0(1)  * 

one  can  include  the  next  group  of  higher  order  terms.  For  this  case  the  various  elements  of  the  transformation  matrix 
S,,,  Eq  (1),  can  be  rewritten  as  (20) 

S„  =  l-i(v,2  +  w,2);  S|2=v,x;  S]3=w,x 

S2|  =-(v„  +4>w,x  +  -^-v,xw,J):  S77  =  1  -  jv,2-0v..w,.. 

Sa  =<>-^'v,2(5;  S3|  =-(w,x  -<5v,x-|-v,2w,x)  (6) 

S32  =— (0  +  v.xw,x  -  -jv.xW;  S33  =  1-yw,2 
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For  ihis  case  explicit  derivation  of  the  equations  of  motion  is  feasible  only  when  using  symbolic  algebraic 
manipulation  on  a  computer. 

Such  explicit  derivations  of  the  equation  of  motion  have  a  number  of  advantages.  First  it  enables  one  to  write 
out  the  equations  of  motion  with  considerable  detail.  Thus  one  can  inspect  the  equations  and  identify  the  coupling 
terms  introduced  by  retaining  the  geometrically  nonlinear  terms.  I'hysical  intcqirctation  of  such  coupling  terms  facili¬ 
tates  the  understanding  of  equations.  Furthermore  equations  derived  by  various  researchers  can  be  compared  with  each 
other  and  differences  between  various  formulations  can  lie  clarified.  Another  advantage  is  associated  with  the  compu 
tational  efficiency  inherent  in  such  formulations.  Numerical  implementation  of  blade  stability  and  response  calculu 
lions  based  on  such  equations  arc  more  efficient  and  computer  storage  requirements  arc  reduced.  This  represents  an 
advantage  in  structural  optimization  studies  of  rotor  blades  with  ncroclastic  constraints  (16-18)  which  rcqu.re  repetitive 
evaluation  of  the  objective  functions  and  the  constraints. 

A  second,  and  more  recent,  approach  to  generating  rotary-wing  equations  of  motion  is  baseo  on  the  implicit 
approach  (11)  In  this  approach  the  equations  of  motion  are  never  explicitly  written  down  and  they  are  generated 
numerically  by  the  computer  during  the  solution  process.  This  approach  is  particularly  effective  when  combined  with 
the  finite  element  approach,  for  the  spatial  discretization  of  the  equation.  This  approach  allows  the  treatment  of  com 
plicated  configurations,  provides  considerable  flexibility  in  the  representation  of  the  aerodynamic  loads  (19).  and  does 
not  require  the  use  of  ordering  schemes  and  the  inherent  approximations  associated  with  such  ordering  schemes. 

More  recent  studies  use  small  strain  and  large  rotation  type  of  analyses  [21-23)  which  utilize  Rodrigues  param¬ 
eters  to  represent  finite  rotations.  This  approach  when  coupled  with  an  implicit  formulation  completely  eliminates  the 
need  for  using  an  ordering  scheme.  Furthermore  it  also  removes  the  need  to  use  the  torsional  quasi-coordinate  which 
was  present  in  a  previous  formulation  (24)  and  created  complexity  in  the  finite  element  models  based  on  this  formula¬ 
tion  Use  of  Euler  angles  combined  with  an  implicit  formulation  is  equally  effective  for  representing  the  large  dis¬ 
placement  of  both  isotropic  or  composite  beams  [25].  These  models  are  more  consistent  and  mathematically  mor- 
elegant  than  blade  models  based  on  ordering  schemes.  At  the  same  time  the  incorporation  of  such  models  into  a  gen¬ 
eral  analysis  is  more  complicated.  Thus  to  date  the  only  aeroelastic  analysis  capability  based  on  such  general  formula¬ 
tion  is  the  GRASP  program  (26). 

Based  on  the  discussion  presented  above  one  can  now  assess  the  importance  of  third  and  higher  order  terms  in 
rotordynamic  models.  For  the  case  of  aeroelastic  stability  of  hingcless  rotor  blades  in  hover  a  recent  study  [27,28]  has 
investigated  the  influence  of  third  order  nonlinear  terms,  similar  to  those  in  Eq.  (6)  when  compared  to  second  order 
nonlinear  terms,  similar  to  those  in  Eq.  (2).  It  was  concluded  that  for  torsionally  soft  blades  (o>n  =  2.5)  and  high  col¬ 
lective  pitch  settings  (0O  >  0.20)  the  effect  of  these  higher  order  terms  on  the  stability  boundaries  was  fairly  small  (less 
than  10%)  Thus  the  effect  of  third  order,  and  higher  order  nonlinear  terms,  on  blade  aeroelastic  stability  appears  to  be 
limited.  The  retention  of  higher  order  nonlinear  terms  can  be  useful  in  the  study  of  mathematical  properties  of  the 
transformations  which  relate  the  position  vectors  of  the  deformed  and  undeformed  states  of  a  blade  undergoing  non¬ 
linear  deformation  (29). 

Another  interesting  question  which  can  be  asked  is  whether  dynamic  models  of  rotor  blades  which  are  based  on 
nonlinear  beam  kinematics  with  finite  large  rotations  are  superior  to  previous  moderate  rotations  theories,  based  on 
ordering  schemes?  There  is  no  unique  answer  to  this  intriguing  question.  When  one  uses  criteria  based  on  mathemati¬ 
cal  elegance,  consistency  and  accuracy  large  deflection  theories  have  a  slight  advantage  because  the  structural  and  iner¬ 
tia  operators  obtained  when  using  these  theories  are  slightly  more  accurate.  However,  it  is  essential  to  note  that  in 
order  to  complete  the  formulation  of  aeroelastic  problem  the  unsteady  aerodynamic  loads  have  to  be  combined  with  the 
inertia  and  the  structural  parts.  Previous  reviews  of  the  unsteady  aerodynamics  used  in  aeroelastic  problems  [2,5,9- 
11,13]  have  clearly  indicated  that  the  aerodynamic  theories  used  for  this  purpose  are  linear  theories,  except  when  one 
attempts  to  model  dynamic  stall,  or  transonic  effects.  Therefore  the  combination  of  kinematical  model  which  contains 
third  order  nonlinear  terms,  with  a  linear  incompressible  aerodynamic  theory  such  as  Greenberg's  theory  [27,28]  does 
not  produce  necessarily  a  consistent  aeroelastic  model.  This  point  is  well  illustrated  by  a  recent  study  by  Hodges, 
Kwon  and  Shankar  [30]  where  it  was  shown  that  by  combining  three  dimensional  tip  loss  and  unsteady  inflow  effects 
with  a  conventional  moderate  deflection  theory  remarkable  agreement  between  theoretical  and  experimental  results 
were  obtained.  In  Fig.  2,  taken  from  Ref.  30,  the  lead-lag  damping  of  a  stiff-in-planc  hingeless  rotor  blade,  with 
configuration  parameters  chosen  from  the  experimental  model  rotor  with  soft  pitch  flexure,  zero  precone  and  droop,  is 
shown  for  various  collective  pitch  angles.  Since  both  theoretical  predictions  have  the  same  structural  model,  the 
differences  are  due  to  aerodynamics.  Clearly  the  results  with  the  three  dimensional  unsteady  aerodynamics,  denoted  as 
panel  method  in  Fig.  2,  give  much  better  agreement  with  the  experimental  data  than  those  based  on  two  dimensional 
quasi-steady  aerodynamics,  denoted  2-D  theory.  This  example  indicates  that  the  key  to  substantia!  improvements  in 
the  aeroelastic  modeling  capability  of  rotor  blades  is  linked  to  improved  unsteady  aerodynamic  models,  and  not  to  die 
retention  of  third  and  higher  order  geometrically  nonlinear  terms. 

Another  basic  inconsistency  in  the  formulation  of  beam  models  for  moderate  and  large  deflections  is  associated 
with  the  fact  that  frequently  excessive  preoccupation  with  geometrically  nonlinear  terms  leads  to  the  neglect  of  simple 
structural  effects  such  as  those  associated  with  Timoshenko  beams,  namely  shear  and  rotary  inertia. 

Using  a  finite  element  model  capable  of  capturing  Timoshenko  beam  effects  in  an  accurate  manner  [31]  and 
calculating  the  influence  of  Timoshenko  beam  effects  on  the  first  three  flapwisc  and  first  three  in-plane  of  rotation  fre¬ 
quencies  of  graphite -epoxy  beam  having  the  cross-sectional  dimensions  shown  in  Fig.  3,  with  a  length  of  15  feet  and 
rotating  at  £2  =  400  RPM,  yields  the  results  presented  in  Table  1.  It  is  evident  that  for  this  particular  case  the  influence 
of  shear  and  rotary  inertia  on  the  lead-lag  or  in-planc  of  rotation  modes,  and  in  particular  the  frequency  of  the  second 
and  third  mode  can  be  significant.  Thus  incorporation  of  such  effects  could  be  of  equal  importance  to  third  order 
geometrical  nonlincarities.  Obviously  the  incorporation  of  true  anisotropy,  present  in  composite  rotor  blades  is  a  much 
more  important  topic.  This  topic  is  the  subject  matter  of  the  next  section. 
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3.  STRUCTURAL  MODELING,  FREE  VIBRATION  AND  AEROELASTiC  ANALYSIS 
OF  COMPOSITE  ROTOR  BLADES 

Mom  of  (lie  structural  models  developed  to  dale  have  been  restricted  to  isotropic  inatetial  properties.  On  the 
other  hand  modem  helicopter  blades  arc  frequently  built  of  composites.  To  remedy  this  situation  a  substantial  share  of 
the  recent  studies  in  this  field  has  been  aimed  at  the  development  of  models  which  are  suitable  for  the  structural  and 
acroelastic  analysis  of  composite  rotor  blades.  The  important  attributes  of  such  a  structural  model,  require  the  capabil 
tty  to  represent  transverse  shear  deformation,  cross  sectional  warping  and  clastic  coupling,  in  addition  to  an  adequate 
representation  of  geometric  nonlmcanties.  Rotor  blades  arc  typically  modeled  as  a  beam.  In  a  beam  tlicurv,  the  defor 
mations  of  the  cross-section,  both  in  and  out  of  the  plane,  are  assumed  to  be  either  small  or  neglected.  Therefore,  an 
approach  commonly  used  in  the  available  structural  models  for  composite  rotor  blade  analysis  is  to  determine  the  cross 
section  warping  functions,  shear  center  location  and  cross  sectional  properties  based  on  a  linear  theory.  The  linear, 
two  dimensional  analysis  for  (lie  cross-section  is  decoupled  from  the  nonlinear,  one  dimensional  global  analysis  for  the 
beam  and  can  be  done  once  for  each  cross  section  of  a  non-uniform  beam.  This  decoupling  is  usually  assumed  in  the 
literature  without  rigorous  proof.  The  discussion  of  composite  rotor  blade  structural  modeling  can,  therefore,  be 
divided  into  two  categories:  (1)  Modeling  approaches  which  lead  to  the  determination  of  the  stiffness  properties  of 
arbitrary  blade  cross  sections.  Anisotropic  materials  and  the  composite  nature  of  the  blade  are  taken  into  account  in  this 
category.  (2)  Structural  models  which  use  an  one-dimensional  beam  kinematics  suitable  for  composite  rotor  blade 
analysis.  A  typical  structural  model  in  this  category  should  include  geometric  nonlincarttics,  pretwist,  transverse  shear 
deformation  and  cross  section  warping.  Many  of  the  existing  composite  rotor  blade  models  in  category  (1)  were 
reviewed  in  detail  in  a  recent  review  paper  by  Hodges  (32).  Our  objective  here  is  not  to  duplicate  I  lodges'  review,  but 
to  consider  the  subject  from  a  slightly  different  perspective. 

Mansfield  and  Sobey  (33]  initiated  the  first  pioneering  study  of  this  difficult  subject.  They  developed  the 
stiffness  properties  of  a  fiber  composite  tube  subjected  to  coupled  bending,  torsion,  and  extension.  Transverse  shear 
and  warping  of  the  cross  section  was  not  included  in  the  model.  This  model  is  too  primitive  to  be  suitable  for  compo¬ 
site  rotor  blade  aeroelastic  analysis.  However,  the  authors  made  some  attempts  to  explore  the  potential  of  this  model 
for  the  innovative  idea  of  aeroelastic  tailoring. 

Rehficld  (34]  used  a  similar  approach  but  included  out-of-plane  warping  and  transverse  shear  deformation. 
This  was  strictly  a  static  theory  for  a  single  cell,  thin  walled,  closed  cross  section  composite,  with  arbitrary  layup, 
undergoing  small  displacements.  This  mode!  was  suitable  for  preliminary  design,  or  frequency  tailoring  studies,  as  well 
as  linear  free  vibration  analysis  of  nonrotating  composite  beams.  Pretwist,  dynamic  effects  and  geometric  nonlinean- 
ties  which  arc  known  to  play  a  major  role  in  helicopter  rotor  dynamics  are  not  accounted  for.  This  relatively  simple 
theory  was  correlated  by  Nixon  [35]  with  experimental  data  and  by  Hodges,  Nixon  and  Rehficld  (36J  with  a  NAS- 
TRAN  finite  element  analysis  of  a  beam  model  with  a  single  closed  cell.  Womdlc  [37]  developed  a  linear,  two- 
dimensional  finite  element  model  to  calculate  the  cross  section  warping  functions  of  a  composite  beam  under 
transverse  and  torsional  shear.  With  these  warping  functions,  the  shear  center  locations  and  the  stiffness  properties  of 
the  cross  section  could  be  calculated.  The  cross  section  can  have  arbitrary  shape  but  the  material  properties  were  res¬ 
tricted  to  monoclinic. 

A  more  general  model  for  calculating  the  shear  center  and  the  stiffness  properties  of  an  arbitrarily  shaped  com¬ 
posite  cross  section  was  developed  by  Kosmatka  [38],  He  used  a  two-dimensional  isoparametric  eight  node  quadrila¬ 
teral  finite  element  model  to  obtain  the  St.  Venant  solution  of  the  cross-section  warping  functions  of  a  tip  loaded  com¬ 
posite  cantilever  beam  with  an  arbitrary  cross  section.  The  beam  was  assumed  to  be  prismatic  (axially  uniform), 
nonhomogeneous,  and  anisotropic.  The  blade  material  was  generally  orthotropic,  i.e.,  orthotropic  material  whose 
material  principal  axes  arc  oriented  arbitrarily.  Therefore,  the  beam  behaves  in  an  anisotropic  manner.  Subsequently 
this  model  was  combined  with  a  companion  moderate  deflection  beam  theory  suitable  for  the  structural  dynamic 
analysis  of  advanced  prop-fan  blades  and  helicopter  rotors,  which  will  be  discussed  later  in  this  section. 

Giavotto,  et  al.  (39]  also  formulated  a  two-dimensional  finite  element  model  for  determining  cross  section 
warping  functions,  shear  center  location  and  stiffness  properties.  A  special  feature  of  this  formulation  is  that  the  result¬ 
ing  equations  have  both  extremity  solutions  and  central  solutions.  The  central  solutions  correspond  to  the  warping  dis¬ 
placements  due  to  applied  loads  without  considering  end  effects,  while  the  extremity  solutions  correspond  to  the  warp¬ 
ing  displacements  due  to  end  effects.  Subsequently  this  work  was  extended  by  Borri  and  Merlini  (40)  to  include  the 
so-called  geometric  section  stiffness  associated  with  large  displacement  formulations.  Bauchau  (41]  developed  a  beam 
theory  for  anisotropic  materials  based  on  the  assumption  that  the  cross  section  of  the  beam  docs  not  deform  in  its  own 
plane  The  out-of-plane  cross  section  warping  was  expressed  in  terms  of  the  so-called  eigenwarpings.  This  theory  is 
valid  for  thin-walled,  closed,  multi-celled  beams  with  transversely  isotropic  material  properties.  Subsequently  it  was 
extended  by  Bauchau,  Coffenbeny  and  Rehfield  (42)  to  allow  for  general  orthotropic  material  properties. 

All  the  studies  described  in  this  section  up  to  this  point  employed  a  separate  two-dimensional  analysis  to  deter¬ 
mine  the  cross  sectional  warping  functions  and  stiffness  properties.  For  non-uniform  beams,  such  a  two-dimensional 
analysis  is  earned  out  once  for  each  cross  section.  A  new  approach  developed  by  Lee  and  Kim  (43)  and  Stcmple  and 
Lee  (44)  uses  a  finite  element  formulation  which  can  represent  thin- walled  beams  with  arbitrary  cross  sections,  general 
spanwise  taper  and  planform  distributions  and  allows  artiitrary  cross  section  warping.  Tins  was  accomplished  by  dis- 
tnbuting  warping  nodes  over  the  cross  section  situated  al  the  node  of  regular  beam  type  finite  element.  Thus  the  treat¬ 
ment  of  the  cross  section  warping  is  coupled  with  the  treatment  of  the  beam  bending,  torsion  and  extension.  This  for¬ 
mulation  considered  only  the  out-of-plane  warping  and  linear  problems,  or  small  deflection  problems. 

Recently,  Stcmple  and  Lee  (45]  have  extended  this  formulation  to  allow  for  large  deflection  static  and  free 
vibration  analysis  of  rotating  composite  beams.  The  disadvantage  of  this  approach  is  that  the  analysis  is  more  expen¬ 
sive  than  those  whose  cross  section  analysis  is  decoupled  from  the  nonlinear  beam  analysis.  Furthermore  the  numerical 
results  obtained  from  this  theory  were  not  compared  to  other  available  theories,  thus  the  validity  of  the  theory  still 
remains  to  be  determined. 
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The  structural  theories  discussed  so  far  emphasize  the  modeling  approach  associated  with  category  (I),  where 
the  emphasis  is  on  determining  the  shear  center,  warping  and  cross-sectional  properties  of  (he  composite  cross  section 
Tor  category  (2)  structural  modeling,  where  the  emphasis  is  one-dimensional  beam  kinematics  suitable  for  the  analysis 
of  composite  rotor  blades,  two  types  of  theories  are  available  depending  on  the  level  of  geometric  nonlinearity  being 
retained  in  the  one  dimensional  beam  kinematics.  The  first  type  is  based  on  a  moderate  deflection  type  of  theory  while 
the  second  type  is  capable  of  modeling  large  deflections.  Moderate  deflection  theories  usually  rely  on  an  ordering 
scheme  to  limit  the  magnitude  of  blade  displacements  and  rotations.  While  large  deflection  theories  do  not  utilize  an 
ordering  scheme  to  limit  the  magnitude  of  blade  displacements  and  rotations,  for  such  theories  the  only  assumption 
used  to  neglect  higher  order  terms  is  the  assumption  that  the  strains  arc  small. 

In  rotary  wing  aeroclasttcity  moderate  deflection  theories  are  usually  adequate  provided  that  a  consistent  order¬ 
ing  scheme  is  used.  The  first  acroelastic  model  for  a  composite  rotor  blade  in  hover  was  presented  in  a  comprehensive 
study  by  Hong  and  Chopra  [46],  In  this  specialized  model,  the  blade  was  treated  as  a  single-cell,  laminated  box  beam 
composed  of  an  arbitrary  lay-up  of  composite  plies.  The  strain-displacement  relations  for  moderate  deflections  were 
taken  from  Hodges  and  Dowell  (24],  which  does  not  include  the  effect  of  transverse  shear  deformations.  Each  lamina 
of  the  laminate  was  assumed  to  have  onhotropic  material  properties.  The  equauons  of  motion  were  obtained  using 
Hamilton's  principle.  A  finite  element  model  was  used  to  discretize  the  equations  of  motion.  Numerical  results  for  the 
coupled  flap-lag-torsional  behavior  of  hingeless  rotor  blades  clearly  illustrated  the  strong  coupling  effects  introduced 
by  the  composite  nature  of  the  blade.  These  coupling  terms  which  depend  on  fiber  orientation  have  a  strong  influence 
on  blade  stability  boundaries  in  hover.  Subsequently  this  analysis  was  extended  to  the  modeling  of  composite  bearing 
less  rotor  blades  in  hover  [47]  and  a  systematic  study  was  carried  out  to  identify  the  importance  of  the  stiffness  cou¬ 
pling  terms  on  blade  stability  with  fiber  orientation  and  for  different  configurations.  In  this  model  the  composite 
flexbeam  of  the  bearingless  rotor  blade  was  represented  by  an  I-secdon  consisting  of  three  laminates.  Each  laminate  ts 
composed  of  an  arbitrary  lay-up  of  composite  plies.  The  outboard  main  blade  and  the  torque  tube  were  assumed  to  be 
made  of  isotropic  materials.  Titus  this  model  represents  a  somewhat  idealized  model  for  a  composite  bearingless  rotor 
blade  In  addition  to  acroelastic  stability  studies  of  composite  rotor  blades  in  hover.  Panda  and  Chopra  (48)  also  stu¬ 
died  the  aeroclastic  stability  and  response  of  hingeless  composite  rotor  blades  in  forward  flight  using  the  structural 
model  presented  in  Ref.  46.  It  was  found  that  ply  orientation  is  effective  in  reducing  both  blade  response  and  hub 
shears. 


A  more  comprehensive  analysis  for  the  structural  dynamic  modeling  of  composite  advanced  prop-fan  blades, 
which  with  some  modifications,  is  also  suitable  for  the  general  modeling  of  curved,  pretwisted  composite  rotor  blades 
was  developed  by  Kosmatka  [38].  The  cross  section  geometry  of  the  blade  is  arbitrary,  and  the  associated  cross  section 
stiffness  properties  and  shear  center  location  can  be  obtained  from  the  accompanying  linear  two  dimensional  finite  ele¬ 
ment  model  which  has  been  discussed  earlier  in  this  section.  In  the  one-dimensional,  nonlinear  analysis,  the  curved 
pretwisted  blade  was  modeled  by  a  series  of  straight  beam  elements  which  are  aligned  with  the  curved  line  of  shear 
centers  of  the  blade.  Each  beam  finite  element  was  derived  using  Hamilton's  principle  and  the  following  basic 
assumptions:  the  beam  has  an  arbitrary  amount  of  pretwist,  undergoes  moderate  deflections,  is  composed  of  generally 
orthotropic  materials,  has  an  arbitrary  cross  sectional  shape,  and  rotates  about  a  vector  in  space.  Numerical  results  for 
frequencies  and  mode  shapes  obtained  from  this  structural  dynamic  model  were  in  good  agreement  with  modal  tests  on 
conventional  and  advanced  propellers  [49,50],  Bauchau  and  Hong  [51-53]  have  developed  a  scries  of  large  deflection 
composite  beam  models  which  arc  intended  for  rotor  blade  structural  dynamic  and  acroelastic  analysis.  The  first  of 
these  models  (51)  used  a  finite  element  approach  combined  with  a  general  global  coordinate  system.  While  the  model 
was  general  it  also  proved  itself  to  be  computationally  inefficient.  A  second  version  of  this  model  [52]  used  a  curvi¬ 
linear  coordinate  system  and  the  resulting  finite  element  model  was  found  to  be  computationally  more  efficient  than  the 
first  one  However  some  deficiencies  associated  with  the  derivation  of  the  strain  displacement  relation  were  noted  by 
Hong  in  his  dissertation  [54],  After  additional  improvements  in  the  model  a  final  version  of  this  theory  which  is  capa¬ 
ble  of  modeling  naturally  curved  and  twisted  beams  undergoing  large  displacements  and  rotations  and  small  strains 
was  developed.  The  kinematics  of  this_ theory  is  an  extension  of  the  common  approach,  using  the  definition  of  Green 
strains,  to  include  effects  such  as  small  initial  curvature,  transverse  shear  deformations  and  out-of-plane  warpings.  The 
fundamental  assumptions  in  the  kinematics  are  the  indefonmability  of  the  cross-section  in  its  own  plane  and  a  revised 
small  strain  assumption.  In  this  revised  small  strain  assumption,  both  axial  and  shearing  strains  arc  still  neglected  com- 
Pij™  3°  however,  no  assumption  is  made  about  the  relative  magnitude  between  the  axial  and  shearing  strains. 
Therefore,  the  second  order  shear  strain  coupling  terms  in  the  axial  strain  expression  are  retained  under  this  revised 
sma  I  strain  assumption.  The  commonly  used  small  strain  assumption,  which  include  an  additional  assumption  that  the 
axial  and  shearing  strains  arc  of  the  same  order  of  magnitude,  was  often  successfully  used  in  beam  models  with  isotro- 
ci  ^  ^msotn?P‘c  materials.  However,  Bauchau  and  Hong  [53]  showed  that  it  might  not  be  adequate  for  beams 
with  highly  anisotropic  material  by  comparing  the  analytical  and  experimental  results  of  a  thin-walled  kevlar  beam. 

1  he  strain-displacement  relations  were  derived  using  an  implicit  formulation  with  seven  unknown  functions  which 
depend  on  the  space  coordinate,  namely  three  displacement  components,  three  rotation  parameters  (Euler  angles  were 
used),  and  the  amplitude  of  the  torsional  warping.  The  three  rotation  parameters,  which  are  used  to  describe  the  large 
rotations  from  the  undeformed  to  the  deformed  triad  of  unit  vectors,  are  defined  implicitly  in  the  rotation  matrix  instead 
o  appearing  explicitly  in  the  strain-displacement  relations.  An  extension  of  this  model  for  free  vibration  analysis  can 
be  tound  in  Hong  s  dissertation  [54],  Using  this  model  a  number  of  cases  testing  the  static  large  deflection  capability 
and  the  free  vibration  capability  of  the  model  were  computed.  However  an  aeroelastic  analysis  of  a  rotor  blade,  based 
on  this  model,  is  not  available  to  date. 


yvcently  N.tr.guct  and  ^ugunuji  t55,56]  have  developed  a  large  deflection  composite  biade  model  lor  static 
pa]  and  free  vibration  [56]  analysis.  Large  deflections  are  accounted  for  by  using  the  Euler  angles  to  describe  the 
transtormation  between  a  global  and  local  coordinate  system  after  deformation.  However,  transverse  shear  deforma¬ 
tion,  and  cross  section  warping  were  not  incorporated  in  this  model.  Thus  the  model  is  more  suitable  for  the  study  of 
fiat  composite  strips  than  actual  rotor  blades.  1 

i  A*''gan  an<i  *  forges  (57)  have  recently  presented  a  theory  for  nonhomogeneoi  s,  anisotropic  beams  undergoing 
ge  global  rotation,  small  local  rotation  and  small  strain.  They  used  a  perturbation  procedure  to  obtain  a  linear  two- 
dimensional  cross  section  analysis  which  is  decoupled  from  the  nonlinear  one-dimensional  global  analysts.  The  non¬ 
linear  beam  kinematics  was  based  on  a  study  by  Danielson  and  Hodges  [58].  The  nonlinear  beam  kinematics,  which 
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defines  a  conjugate  stress  measure,  describes  a  constitutive  relationship  and  develops  the  momentum  balance  condi 
lions,  was  based  on  Atluri  (59).  TTie  cross  section  warping,  both  near  the  ends  (boundary  layer  solutions)  and  away 
from  the  ends  (St.  Vcnant  solutions)  was  obtained  in  terms  of  force  and  moment  strains  from  a  linear  two  dimensional 
equation.  The  three  force  strains  include  axial  and  both  transverse  shear  strains  at  the  reference  line,  the  three  moment 
strains  are  torsion  and  curvature  quantities  (58).  Since  the  strain  components  from  the  nonlinear  beam  kinematics  are 
expressed  in  terms  of  force  and  moment  strains,  and  cross  seetton  warpings,  it  is  possible  to  obtain  a  one  dimensional 
system  of  equations  in  terms  of  force  and  moment  strains  only  by  substituting  the  solutions  for  the  cross  sectional 
warpings  into  the  strain  expressions.  However,  the  authors  abandoned  this  procedure  as  being  too  tedious  and  unncces 
sary.  Instead,  they  chose  to  follow  the  same  approach  as  Giavotto,  et  al.  (39).  which  solves  (he  cross  section  warping 
functions  in  terms  of  force  and  moment  stress  resultants.  A  geometric  stiffness  matrix  was  also  formulated  in  this 
model.  Tins  matrix  and  the  material  stiffness  matrix  were  intended  to  be  used  in  an  incremental  updated  Lagrangian 
formulation  for  large  deflections.  A  general  mixed  variational  principle  was  used  in  the  formulation  of  the  one- 
dimensional,  geometrically  nonlinear  global  analysis.  No  numerical  results  were  presented  in  this  paper.  This  study 
draws  rather  heavily  on  the  material  presented  in  Refs.  39  and  59  and  its  principal  contribution  seems  to  be  the  dcriva 
tion  of  a  geometric  stiffness  matrix. 

Prom  this  review  of  the  literature  on  available  structural  models  capable  of  represenung  composite  rotor  blades 
it  is  evident  that  this  has  been  an  active  area  of  research  during  the  last  five  years.  Yet  it  is  remarkable  to  note  that 
despite  the  availability  of  such  models  the  only  published  body  of  research  which  actually  contains  acroelastic  stability 
and  response  type  of  results  is  that  published  by  Chopra  and  his  associates  (46-48).  This  is  somewhat  disappointing 
since  the  potential  for  acroelastic  tailoring  for  composite  rotor  blades  has  been  demonstrated  (46).  A  typical  result.  Fig 
4,  taken  from  Ref.  48,  shows  that  for  a  four  bladed  hingeless  soft-in-plane  composite  rotor,  where  the  composite  single 
cell  box  has  symmetric  laminates,  the  ply  orientation  angle,  At,  for  the  vertical  walls  of  the  box  has  a  very  significant 
cfTect  on  flic  peak-to-peak  value  of  nondimensional  vertical  hub  shears.  This  quantity  is  indicative  of  the  vibration  lev¬ 
els  experienced  at  the  hub.  It  is  interesting  to  note  that  both  positive  and  negative  ply  angles  reduce  the  oscillatory  hub 
shear,  however  positive  angles  of  30°  are  more  effective  in  reducing  vibration  levels.  In  view  of  the  positive  nature  of 
coupling  effects  present  in  composite  blades  one  would  expect  to  see  more  studies  of  this  type  which  take  advantage  of 
the  vibratory  load  reduction  potential  present  in  composite  blades. 

Instead  of  the  excessive  preoccupation  with  arbitrarily  large  deflection  theories  for  composite  beams  it  appears 
to  be  more  sensible  to  derive  reahsne  composite  blade  models,  based  on  moderate  deflection  theories,  and  use  them  to 
explore  the  potential  for  designing  composite  rotor  blades  which  have  low  vibration  levels  and  good  acroelastic  stabil 
ity  margins.  Based  on  the  present  state-of-the-art  it  appears  that  a  moderate  deflection  type  of  beam  theory  combined 
with  a  cross-sectional  analysis  of  the  composite  blade  presented  in  Refs.  38  or  39  would  be  both  suitable  and  adequate 
for  this  purpose. 


4.  MODELING  OF  COUPLED  ROTOR/FUSELAGE  AEROMECHANICAL  PROBLEMS 
AND  THEIR  ACTIVE  CONTROL 

The  acromechanica!  instability  of  a  helicopter,  on  the  ground  or  in  flight,  is  caused  by  coupling  between  the 
rotor  and  body  degrees  of  freedom,  this  instability  is  commonly  denoted  air  resonance  when  the  helicopter  is  in  flight 
and  ground  resonance  when  the  helicopter  is  on  the  ground.  The  physical  phenomenon  associated  with  this  instability 
is  quite  complex.  The  rotor  lead-lag  regressing  mode  usually  couples  with  the  body  pitch  or  roll  to  cause  an  instability. 
The  nature  of  the  coupling  which  is  both  aerodynamic  and  inertial  is  introduced  in  the  rotor  by  body  or  support  motion. 
The  importance  of  developing  a  mathematically  consistent  model  capable  of  representing  the  coupled  rotor/fuselage 
dynamic  system  has  already  been  discussed  in  previous  reviews  (5,6,9-11,13).  With  the  use  of  algebraic  symbolic 
manipulative  programs  or  implicit  formulations  (26,60)  sophisticated  mathematical  models  for  coupled  rotor/fusclage 
aeromechanical  problems  can  be  formulated  and  solved. 

The  primary  purpose  of  this  section  is  to  discuss  a  number  of  modeling  aspects  for  this  class  of  problems  which 
have  been  shown  to  be  important  by  recent  research.  Subsequently  research  on  the  use  of  active  controls  to  stabilize 
this  class  of  aeromechanical  problems  is  also  described. 

Many  previous  air  resonance  studies  (6,9- 11,13]  were  limited  by  simplifying  assumptions  such  as  hover,  blades 
which  were  assumed  to  be  torsionally  rigid  or  aerodynamic  loads  based  on  quasisteady  assumptions.  Recent  studies 
(61,62)  have  clarified  the  role  of  unsteady  aerodynamics,  forward  flight  and  torsional  flexibility  on  air  resonance. 
Furthermore  the  mathematical  model  derived  for  this  coupled  rotor  fuselage  system  had  also  a  provision  for  including 
an  active  controller  capable  of  suppressing  air-resonance. 

The  mathematical  model  of  the  rotor/fuselage  system  of  Ref.  61  and  62,  and  its  salient  features  arc  described 
next  The  fuselage  is  represented  as  a  rigid  body  with  live  degrees  of  freedom,  where  three  of  these  are  linear  transla¬ 
tions  and  two  are  angular  positions  of  pitch  and  roll  (Fig.  5).  Yaw  is  ignored  since  its  effect  in  the  air  resonance  prob¬ 
lem  is  known  to  be  small.  A  simple  offset  hinged  spring  restrained  rigid  blade  model  is  used  to  represent  a  hingeless 
rotor  blade  (Fig.  6).  This  assumption  simplifies  the  equations  of  motion,  while  retaining  the  essential  features  of  the  air 
resonance  problem  In  this  model,  the  blade  elasticity  is  concentrated  at  a  single  point  called  the  hinge  offset  point,  and 
torsional  springs  are  used  to  represent  this  flexibility.  The  dynamic  behavior  of  the  rotor  blade  is  represented  by  three 
degrees  of  freedom  for  each  blade,  which  arc  flap,  lag,  and  torsion  motions  The  aerodynamic  loads  of  the  rotor  blades 
arc  based  on  quasi-steady  Greenberg’s  theory,  which  is  a  two  dimensional  potential  flow  strip  theory.  Compressibility 
and  dynamic  stall  effects  arc  neglected,  though  they  could  be  important  at  high  advance  ratios.  Unsteady  aerodynamic 
effects,  which  arc  created  by  the  time  dependent  wake  shed  by  the  airfoil  as  it  undergoes  arbitrary  time  dependent 
motion,  are  accounted  for  by  using  a  dynamic  inflow  model.  This  simple  model  uses  a  third  order  set  of  linear 
differential  equations  driven  by  perturbations  in  the  aerodynamic  thrust,  roll  moment,  and  pitch  moment  at  the  rotor 
hub  The  three  states  of  these  equations  describe  the  behavior  of  perturbations  in  the  induced  inflow  through  the  rotor 
plane.  The  dynamic  inflow  mode!  coefficients  used  are  those  of  Ref,  63. 
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The  equations  of  motion  of  the  coupled  rotor/fuselage  system  are  very  large  and  contain  geometrically  non¬ 
linear  terms  due  to  moderate  blade  deflections  in  the  aerodynamic,  inertial,  and  structural  forces.  Furthermore,  the 
coupled  rotor/fuselage  equations  have  additional  complexity  due  to  the  presence  of  the  fuselage  degrees  of  freedom. 
Tu  reduce  the  equations  to  a  manageable  sire,  an  ordering  scheme  is  used  in  the  derivation  of  the  equations  of  motion 
to  systematically  remove  the  higher  order  nonlinear  terms.  The  ordering  scheme  is  based  on  Eq.  (3)  For  this  class  of 
problems  c  represents  the  slopes  of  the  deflections  of  the  blades,  winch  usually  are  of  an  order  of  magnitude  which  is 
less  than  .15.  The  blade  degrees  of  freedom  arc  assigned  an  order  of  e,  while  the  fuselage  degrees  of  freedom  are  of 
order  e’72,  the  various  other  parameters  have  the  order  of  magnitude  given  in  Eqs.  (4).  A  symbolic  manipulation  pro¬ 
gram  is  then  used  to  generate  the  nonlinear  set  of  equations  of  the  rotor/fuselage  system  using  the  ordering  scheme. 
Five  fuselage  equations  result  of  which  three  enforce  the  fuselage  translational  equilibrium  and  two  enforce  the  roll 
and  pitch  equilibrium  The  .lircc  resulting  rotor  blade  equations  arc  associated  with  flap,  lag,  and  torsional  motions  of 
each  blade.  Also,  the  aerodynamic  thrust  and  roll  moments  at  the  hub  center  arc  determined  for  the  perturbation  aero¬ 
dynamics  in  the  dynamic  inflow  equation.  All  of  these  equations  can  be  found  in  detail  in  Ref.  62. 

An  active  control  device  to  suppress  the  air  resonance  instability  through  a  conventional  swashplate  is  also 
incorporated  in  this  mathematical  model.  The  pitch  of  the  k-th  rotor  blade  is  given  by  the  expression 

OpV  =  (Oo  +  AOo)  +  (0)c  +  AG|C)cos(vji|.)  +  (0|,  +  A0i,)sin(\{ik)  ^ 

The  various  pitch  terms  with  the  symbol  A  arc  small  and  these  represent  the  active  control  inputs,  while  those  without 
A  are  the  inputs  necessary  to  trim  the  vehicle. 

The  stability  of  the  system  is  determined  through  the  linearization  of  the  equations  of  motion  about  a  blade 
equilibrium  solution  and  the  helicopter  trim  solution.  The  helicopter  trim  and  equilibrium  solution  are  extracted  simul¬ 
taneously  using  harmonic  balance  for  a  straight  and  level  flight  condition  (61).  After  linearization,  a  multi-blade  coor¬ 
dinate  transformation  is  applied,  which  transforms  the  set  of  rotating  blade  degrees  of  freedom  to  a  set  of  hub  fixed 
non-rotating  coordinates.  This  transformation  is  introduced  to  take  advantage  of  the  favorable  properties  of  the  non¬ 
rotating  coordinate  representation.  The  original  representation  has  periodic  coefficients  with  a  fundamental  frequency 
of  unity,  however,  the  transformed  system  has  coefficients  with  a  higher  fundamental  frequency.  These  higher  fre¬ 
quency  periodic  terms  have  a  reduced  influence  on  the  behavior  of  the  system  and  can  be  ignored  in  some  analyses  at 
low  advance  ratios  (5).  In  hover,  the  original  system  has  period  coefficients  with  a  frequency  of  unity,  but  the 
transformed  system  has  constant  coefficients. 

Once  the  transformation  is  carried  out,  the  system  is  rewritten  in  first  order  form. 

x  =  A(vp)x  +  B(\p)u  ® 

The  fundamental  frequency  of  the  coefficient  matrices  depends  on  the  number  of  rotor  blades.  For  an  odd  bladed  sys¬ 
tem  the  fundamental  frequency  is  Nb  per  revolution,  while  for  an  even  bladed  system  the  fundamental  frequency  is 
Nb/2  per  revolution.  Stability  can  now  be  determined  using  either  an  eigenvalue  analysis  or  Floquet  theory  for  the 
periodic  problem  in  forward  flight.  An  approximate  stability  analysis  in  forward  flight  is  also  possible  by  performing 
an  eigen-analysis  on  the  constant  coefficient  portion  of  the  system  matrices  in  Eq.  (8). 

The  mathematical  model  was  carer  'y  tested  by  comparing  results  to  other  investigator's  analytical  and  experi¬ 
mental  results.  The  correlation  with  thes.  results  was  good  and  verified  that  the  effects  of  torsion,  unsteady  aero¬ 
dynamics,  and  forward  flight  were  accurat  'y  represented  in  the  model  [61,63). 

This  mathematical  model  was  used  to  analyze  the  behavior  of  a  four  bladed  hingclcss  helicopter  somewhat 
similar  to  the  MBB  105  helicopter.  The  nominal  configuration  differs  from  the  MBB  105  in  that  it  has  an  unstable  air 
resonance  mode.  The  system  has  37  states.  The  five  body  degrees  of  freedom  and  the  twelve  rotor  degrees  of  freedom 
(three  degrees  of  freedom  for  each  blade)  produce  34  position  and  rate  states.  The  dynamic  inflow  model  augments  the 
system  with  three  more  states  giving  a  total  system  order  of  37.  Thus  even  this  relatively  simple  model  requires  a  con¬ 
siderable  number  of  states  (or  degrees  of  freedom). 

As  mentioned  before  the  lead  lag  regressing  mode  is  the  critical  degree  of  freedom  associated  with  air  reso¬ 
nance  therefore  the  essential  features  of  this  instability  arc  described  by  damping  plots  for  this  particular  degree  of 
freedom. 

Figure  7  illustrates  the  influence  of  unsteady  aerodynamics  as  well  as  the  effect  of  periodic  coefficients  (or  for¬ 
ward  flight)  on  tiie  lead-lag  regressing  mode  damping  of  the  open  loop  configuration.  The  two  sets  of  curves  represent 
air  resonance  damping  of  the  configuration  with  quasi-steady  aerodynamics  and  with  dynamic  inflow  at  various 
advance  ratios  Dynamic  inflow  captures  primarily  the  low  frequency  unsteady  aerodynamic  effect  which  is  important 
for  rotor/fuselage  aeromechanical  problems  such  as  air  resonance.  The  stabilizing  effect  in  forward  flight,  which  is  evi¬ 
dent  in  the  figure,  is  consistent  with  behavior  observed  in  previous  studies  [5,1 1).  For  hover,  the  system  has  constant 
coefficients  and  thus  the  constant  coefficient  approximation  and  the  periodic  system  produce  the  same  results,  as  is 
clearly  evident  in  the  figure.  It  is  also  shown  in  the  figure  thatthe  effect  of  periodic  coefficients  is  relatively  minor. 
The  quasisteady  aerodynamic  model  produces  a  more  stable  system  than  the  model  which  includes  the  unsteady  aero¬ 
dynamic  effects  as  represented  by  the  dynamic  inflow  model.  It  is  also  worthwhile  mentioning  that  considerable 
differences  between  the  two  models  exist  particularly  in  low  advance  ratios. 

Figure  8  shows  that  neglecting  the  torsional  degree  of  freedom  on  the  nominal  configuration  increases  the  insta¬ 
bility  of  the  lead-lag  regressing  mode.  The  trend  of  the  two  curves  also  tends  to  diverge  at  high  advance  ratios.  The 
addition  of  torsion  also  tends  to  amplify  the  effect  of  the  periodic  terms.  At  high  values  of  advance  ratio,  the  llap-lag- 
torsion  model  shows  a  much  greater  difference  between  the  constant  and  periodic  stability  analysis  than  does  the  flap- 
lag  analysis. 
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Coupled  rotor  fuselage  analyses  based  on  an  implicit  formulation  developed  by  Done  ct  al.  |60]  have  modeled 
the  cITcct  of  swashplaic  flexibility  and  control  system  stillness  on  this  class  of  problems.  The  limited  number  of 
numerical  results  presented  [60]  do  not  comey  a  clear  picture  of  the  influence  of  these  additional  parameters  on  the  an 
resonance  problem. 

Another  recent  study  by  l.oewy  and  Zotto  (64)  studied  the  effect  of  rotor  shaft  flexibility  and  associated  rotor 
control  coupling  on  the  ground'nir  resonance  of  helicopters.  This  is  of  particular  interest  for  certain  lypcs  of  advanced 
helicopters  scinch  have  a  relatively  flexible  shaft.  Tlic  analysis  was  formulated  using  a  l-agrangian  approach  and  the 
equations  were  derived  using  symbolic  manipulation  based  on  the  MACSYMA  package.  The  equations  are  based  ini¬ 
tially  on  large  Euler  angles  however  subsequently  the  equations  were  simplified  using  an  ordering  scheme.  I  he  blades 
sserc  offset  hinged  spring  restrained  with  flap  and  lag  degrees  of  freedom  for  each  blade,  and  a  four  bladed  rotor  seas 
considered  The  model  included  a  total  of  12  degrees  of  freedom.  Two  cyclic  flap  modes,  tsvo  cyclic  lag  modes, 
fuselage  pitch  and  roll,  fuselage  center  of  mass  translation  in  two  directions,  hub  translation  due  to  shaft  flexibility  in 
two  directions,  and  shaft  bending  slope  in  two  directions.  The  blades  were  assumed  to  be  torsionally  stiff.  Only  hover 
was  considered,  the  aerodynamic  loads  were  based  on  simple  strip  theory,  and  the  periodic  coefficients  weres  elim¬ 
inated  using  multibladc  coordinates.  Numerical  results  were  obtained  for  a  configuration  similar  to  the  OH-58D  hel¬ 
icopter  with  a  four  bladed  articulated  rotor.  The  effect  of  shaft  flexibility  and  associated  coupling  terms,  together  with 
the  influence  of  a  mass  simulating  a  mast  mounted  sight  on  air  and  ground  resonance  stability  boundaries  was  studied 
with  considerable  detail.  It  was  found  that  flcxibility/control  coupling  adds  new  modes  of  instability  for  ground  reso¬ 
nance,  however  this  was  found  to  be  a  weak  instability  which  was  eliminated  by  small  amounts  of  structural  damping. 
Air  resonance  type  of  instabilities  arc  more  susceptible  to  shaft  flexibility/control  coupling  effects  and  small  amount  of 
structural  damping  fails  to  stabilize  the  coupled  rotor/fusclage  system.  Therefore  for  certain  combinations  of  parame¬ 
ters  shaft/blade  pitch  coupling  effects  have  to  be  designed  carefully  to  insure  the  stability  of  the  air  resonance  mode 

The  air  resonance  stability  of  hingelcss  rotors  in  forward  flight  was  also  studied  in  Ref.  76.  The  body  of 
research  available  on  coupled  rotor/fuselage  acromechanical  problems  has  reached  a  remarkable  level  of  maturity  in  a 
fairly  short  period  of  time.  It  is  also  evident  that  reliable  models  for  this  class  of  problem  should  contain  coupled  flap- 
lag-torsional  blades  models  combined  with  a  fuselage  which  has  pitch  and  roll  as  well  as  two  translational  degrees  of 
freedom  For  certain  configuration;  coupling  effects  introduced  by  swashplate  flexibility,  shaft  bending  and  pitch  link 
flexibility  should  be  also  incorporated  in  the  mode!  so  as  to  obtain  reliable  stability  boundaries.  Tlic  aerodynamic 
representation  which  is  most  suitable  for  (his  class  of  problem  is  the  dynamic  inflow  model  (63).  Up  to  advance  ratios 
of  p  =  0  40  the  role  of  periodic  coefficients  is  fairly  small,  and  usually  the  most  unstable  cases  for  air  resonance  occur 
for  the  case  of  hover  For  soft-in  plane  hingelcss  helicopter  configurations  and  articulated  blade  configurations  forward 
flight  usually  introduces  a  stabilizing  effect. 

The  derivation  of  coupled  rotor/fusclage  analyses  or  models  has  become  fairly  routine  when  using  computer 
algebra  or  implicit  formulations.  Therefore  studies  based  on  very  simple  models  which  contain  fewer  degrees  of  free¬ 
dom  than  those  retained  in  Refs.  61-62  and  64  will  rarely  serve  a  useful  purpose,  since  the  validity  of  such  models  is 
questionable. 

Improved  understanding  of  acromechanical  phenomena  such  as  air  and  ground  resonance  has  also  raised  the 
possibility  of  eliminating  or  suppressing  such  instabilities  using  active  controls.  One  possible  means  of  stabilizing  or 
augmenting  stability  of  air  resonance  is  through  an  active  controller  operating  with  a  conventional  swashplate.  This 
approach  is  feasible  from  a  practical  point  of  view  only  if  it  is  simple  to  implement  since  it  must  compete  against  the 
straightforward  mechanical  solution  to  this  problem  based  on  lag  dampers.  Such  an  active  controller  would  need  sens¬ 
ing  and  actuating  devices  leading  to  an  expensive  system.  However,  with  the  inevitable  introduction  of  other  active 
control  devices  such  as  higher  harmonic  control  (HHC)  for  vibration  suppression  [65-67)  this  argument  is  considerable 
weakened.  Vibration  control  requires  sensors  and  actuators  with  bandwidths  well  above  the  1/rev  frequency.  Since. die 
air  resonance  instability  results  in  an  unstable  lead-lag  regressing  mode  (i.c.,  the  mode  associated  with  the  [1  -  coLi) 
frequency)  these  devices  would  also  be  sufficient  for  air  resonance  control.  Thus,  sensing  and  actuator  hardware, 
which  may  be  already  available,  could  be  used  for  additional  purposes  below  the  frequency  range  intended  for  the 
available  vibration  control  objective. 

Research  in  the  active  control  of  air  and  ground  resonance  has  been  limited  to  a  few  studies  [68-70],  where 
various  theoretical  aedve  control  studies  were  presented.  The  helicopter  models  used  in  these  studies  were  quite  lim¬ 
ited  since  important  effects  such  as  torsional  flexibility  of  the  rotor  blades,  forward  flight,  and  unsteady  aerodynamics 
were  all  neglected.  Furthermore,  the  studies  dealing  with  the  active  control  of  air  resonance  did  not  adequately  demon¬ 
strate  the  ability  of  the  control  schemes  to  operate  through  the  wide  range  of  operating  conditions  which  can  normally 
be  encountered. 

A  comprehensive  study  which  demonstrated  the  feasibility  of  designing  a  simple  active  controller  capable  of 
suppressing  air  resonance  throughout  the  flight  envelope  representative  of  a  wide  range  of  operating  conditions  which 
may  be  encountered  by  a  helicopter,  with  hingelcss  blades,  has  been  completed  recently  [61,62,71).  The  coupled 
rotor/fusclage  model  representing  this  system  was  briefly  described  at  the  beginning  of  this  section.  The  equations  of 
motion  describing  this  system  are  represented  by  Eq.  (8),  and  the  pitch  of  the  k-th  blade,  determined  by  the  controller  is 
expressed  by  Eq.  (7). 

The  control  studies  undertaken  consisted  of  two  stages.  In  the  first  stage  (61)  the  coupled  rotor/fuselage  system 
was  combined  with  linear  quadratic  optimal  control  theory  to  design  full  slate  feedback  controllers.  These  controllers 
were  then  used  to  evaluate  the  importance  of  various  modeling  effects  on  the  closed  loop  damping  of  the  unstable  air 
resonance  mode.  It  was  found  that  periodic  terms  in  the  model  seem  to  play  only  a  small  role  at  advance  ratios  less 
than,  p  =  0  40  With  this  in  mind  and  considering  the  cost  of  extracting  periodic  optimal  gains,  it  seems  reasonable  to 
neglect  the  periodic  terms  in  the  initial  stage  of  controller  development.  Knowing  that  the  constant  model  is  a  reason¬ 
able  approximation  also  allows  the  use  of  many  other  control  design  techniques.  The  results  also  indicate  that  unsteady 
aerodynamics  and  blade  torsional  flexibility  seem  to  be  important  modeling  effects  that  should  be  included  in  a  con¬ 
troller  design  model.  Significant  errors  of  between  25  and  50  percent  in  closed  loop  lead-lag  damping  could  result  in 
these  effects  arc  not  included.  The  collective  control  input  scents  to  have  little  influence  in  controlling  air  resonance  at 
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high  advance  ratios,  so  it  was  fell  to  be  unnecessary  for  the  complete  control  task,  finally,  partial  state  feedback  of  the 
body  states  does  not  seem  to  be  a  reasonable  approach  to  controlling  air  resonance.  Poor  lead  lag  damping  results  and 
lead  lag  progressive  mode  excitation  is  a  possible  consequence.  This  particular  conclusion  contradicts  the  behavior 
observed  for  articulated  rotors  in  Ref.  69. 

Since  full  state  feedback  is  not  practical  and  (he  first  stage  of  the  research  showed  that  partial  state  feedback  is 
not  reliable,  the  second  stage  of  this  research  162,71)  was  based  on  more  advanced  control  system  synthesis  techniques 
to  design  a  suitable  controller. 

The  controller  .itned  at  suppressing  air  resonance  in  the  flight  envelope  of  the  helicopter  is  based  on  an  optimal 
state  estimator  in  conjunction  with  optimal  feedback  gains  |72).  A  constant  coefficient  model  is  assumed  since  the 
results  obtained  in  the  first  stage  indicated  that  a  periodic  model  was  unnecessary.  The  objective  in  this  portion  of  the 
study  was  to  design  a  controller  at  an  operating  condition  and  require  it  to  function  adequately  at  off  design  conditions, 
corresponding  to  the  entire  flight  envelop  of  the  helicopter.  Furthermore  in  all  applications,  the  design  model  and  (he 
actual  plant  to  be  controlled  will  base  unavoidable  differences  due  to  the  limitations  associated  with  formulating 
models  of  physical  systems.  To  overcome  these  difficulties,  die  multivariable  frequency  domain  design  methods  of 
Refs  73.74  and  75  were  used.  This  allows  interpretation  of  the  design  process  using  frequency  domain  concepts  and 
accounts  for  the  possibility  of  high  frequency  modeling  error.  All  this  can  be  accomplished  while  retaining  the  struc¬ 
ture  of  the  state  space  approach.  The  technique,  which  is  based  on  transfer  function  singular  values,  proved  to  be  par¬ 
ticularly  effective  in  resolving  problems  that  would  not  be  obvious  if  only  the  covanancc  and  weight  matrices  were 
used  in  the  design  process.  To  select  the  design  loop  shapes,  Loop  Transfer  Recovery  was  used,  which  can  be  inter¬ 
preted  as  an  optimization  balancing  system  performance  requirements  and  the  requirement  of  stability  in  presence  of 
modeling  errors. 

The  controller  used  a  single  roll  rate  measurement  and  both  the  sine  and  cosine  swashplate  inputs.  This 
configuration  is  particularly  simple  since  the  measurement  is  taken  from  a  non-rotating  frame  of  (the  fuselage)  refer¬ 
ence  avoiding  the  need  to  send  signals  across  the  rotor  head.  Using  sine  and  cosine  inputs  is  also  simple  and  can  be 
accomplished  through  a  conventional  swashplate  mechanism.  A  constant  four  mode  design  model  consisting  of  the 
body  roll  and  pitch  modes  and  the  lead-lag  regressing  and  progressing  modes  was  found  to  be  quite  practical  for  con¬ 
trol  design  The  controller  was  shown  to  stabilize  the  system  throughout  a  wide  range  of  loading  conditions  and  for¬ 
ward  flight  speeds  and  it  required  small  inputs  of  the  order  of  three  degrees  or  less. 

Some  typical  results  obtained  in  this  study,  for  the  same  four  bladcd  hingeless  rotor  configuration  for  which 
results  were  presented  in  Figs.  7  and  8,  are  presented  next.  The  open  loop  lead-lag  regressing  mode  damping  of  the 
helicopter  configuration  throughout  its  flight  regime  is  presented  in  Fig.  9.  The  horizontal  axis  is  the  advance  ratio, 
while  the  vertical  axis  is  the  fuselage  mass  Mp  nondimensionalized  by  the  blade  mass  of  52  kg.  A  nondimensional 
fuselage  mass  of  32  plus  four  blades  corresponds  to  the  nominal  total  mass  of  1872  kg.  The  figure  indicates  that  the 
system  experiences  an  air  resonance  instability  throughout  most  of  the  flight  regime.  Marginal  stability  exists  at  an 
advance  ratio  greater  than  (l  =  0.35  and  the  point  of  deepest  instability  is  at  Mp  =  30  and  in  the  vicinity  of  hover.  Fig¬ 
ure  10  shows  the  same  system  after  the  controller  designed  according  to  the  methodology  [71]  discussed  above,  has 
been  applied  on  the  helicopter.  From  the  figure  it  is  clear  that  the  lead  lag  regressing  mode  is  stable  throughout  the 
whole  flight  regime,  and  its  stability  is  lowest  in  the  neighborhood  of  Mp  =  23  and  p  =  0.1 1. 

The  body  of  research  described  above  indicates  that  during  the  last  five  years  considerable  progress  has  been 
made  in  understanding  the  role  of  active  controls  as  a  potential  means  for  stabilizing  air  and  ground  resonance.  It 
should  be  noted  that  air  resonance  is  a  relatively  mild  instability,  while  ground  resonance  is  substantially  stronger  thus 
the  former  is  easier  to  control.  It  is  interesting  to  note  that  partial  state  feedback  of  the  body  states  docs  not  seem  to  be 
a  reasonable  approach  to  controlling  air  resonance  in  four  bladed  hingeless  rotors  (61).  This  is  probably  due  to  the  fact 
that  for  hingeless  rotors  the  coupling  between  the  blade  and  fuselage  degrees  of  freedom  is  much  stronger.  A  relatively 
simple  controller  using  only  a  single  body  roll  rate  measurement  and  two  swashplate  control  inputs  (sine  and  cosine) 
was  shown  to  stabilize  the  air  resonance  throughout  the  whole  flight  envelope  when  using  modem  state  of  the  art  con¬ 
trol  system  design  techniques  [62,71].  It  is  precisely  this  robustness  of  the  modem  control  system  synthesis  techniques 
employed  which  make  them  attractive  when  compared  to  more  classical  control  techniques  employed  by  Ham  and  his 
associates  in  their  research  on  the  use  of  individual-blade-control  for  stability  augmentation  [77]. 

Finally  it  should  be  noted  that  while  the  air  resonance  instability  proved  itself  to  be  fairly  easy  to  control,  the 
active  control  of  more  powerful  instabilities  such  as  the  fuselage  induced  flap-lag  instability,  and  the  coupled  flap-lag- 
torsional  instability  of  stiff-in-plane  hingeless  rotor  blades  in  forward  flight,  which  was  also  studied  in  Ref.  62,  proved 
itself  to  be  much  more  difficult  to  control. 


5.  THE  USE  OF  HIGHER  HARMONIC  CONTROLS  FOR  VIBRATION  REDUCTION 
OF  HELICOPTER  ROTORS  IN  FORWARD  FLIGHT 
i>  use  of  high  frequency  blade  pitch  inputs,  referred  to  as  higher  harmonic  control  (HHC),  to  reduce  hel¬ 
icopter  vibrations  has  been  investigated  in  a  number  of  studies.  Aircraft  flight  tests  [67,78,79],  wind  tunnel  tests  [80- 
83],  and  analytical  simulations  [84-91]  have  shown  that  HHC  is  capable  of  substantial  reduction  in  helicopter  vibiauon 
levels  encountered  in  forward  flight. 

Furthermore  the  literature  in  this  field  has  been  surveyed  in  a  number  of  review  articles  (92,1 1).  Particularly 
noteworthy  is  Johnson’s  review  article  [92],  which  described  the  state  of  the  an  up  to  1982  and  a  recent  paper  by  Shaw 
et  al  [82]  which  provides  an  excellent  perspective  on  this  important  topic.  Currently  it  is  well  understood  that  HHC 
produces  vibration  reduction  by  modifying  the  unsteady  aerodynamic  loads  on  the  rotor  blades. 

it  is  quite  remarkable  that  a  significant  portion  of  the  research  in  this  area  involved  wind  tunnel  testing  [80-83] 
and  flight  testing  [67,78,79].  While  a  number  of  analytical  studies  were  carried  out  [84-91]  they  were  less  comprehen¬ 
sive  than  the  tests  Previous  analytical  studies  have  generally  relied  on  simple  analog  [87]  or  frequency  domain  [87] 
models  of  the  helicopter  response.  Other  studies  (84-86,90,91)  were  based  on  using  a  fairly  old  acroclastic  response 
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code.  the  G400  (93,94)  for  simulation  purposes.  Hie  model  used  in  these  studies  did  not  have  a  consistent  represents 
tion  of  the  geometrical  nonlincantics  due  to  moderate  deflect1  ns  Oilier  shortcomings  of  this  simulation  capacity  were 
the  lack  of  time  domain  aerodynamics  needed  for  captunng  high  frequency  unsteady  aerodynamic  effects  as  well  as  a 
step  by  step  time  integration  method  which  precluded  the  calculation  of  direct  stability  information  in  forward  flight, 
such  as  provided  by  Floquet  theory. 

Recently  a  comprehensive  acroclastic  simulating  capability  has  Isccn  developed  (95-97]  and  used  to  study  a 
number  of  fundamental  issues  in  higher  harmonic  control.  The  principal  topics  studied  using  this  new  acroclastic 
simulation  capability  are  listed  below: 

1  Comparison  of  the  effectiveness  of  deterministic  and  cautious  controllers  based  on  local  and  global  111  1C 
models  in  reducing  vibratory  hub  shears. 

2.  Comparison  of  the  response  of  HHC  of  roughly  equivalent  amculaicd  and  liingcless  rotors. 

3  Evaluation  of  changes  in  hub  moments  when  hub  shears  arc  minimized  using  HI  1C.  Similar  calculations  are 
done  when  HHC  is  used  to  try  to  reduce  both  hub  moments  and  hub  shears  simultaneously.  These  studies  are 
done  for  both  articulated  and  hingeless  blades. 

4  Determination  of  the  influence  of  HHC  on  the  acroclastic  stability  margins  of  the  blade  in  forward  flight. 

5  Comparison  of  the  relative  additional  power  requirements  experienced  when  I II IC  is  applied  to  two  similar 
rotor  configurations,  one  with  articulated  and  one  with  hingcless  blades. 


A  brief  description  of  the  acroclastic  simulation  capability,  the  implementation  of  HHC  and  few  important  con 
elusions  obtained  in  the  course  of  these  studies  [95-97]  are  presented  next. 

The  coupled  flap-lag-torsional  equations  of  motion  which  serve  as  the  basis  of  this  analysis  arc  similar  to  those 
derived  in  Ref.  98.  They  contain  geometrically  nonlinear  terms  due  to  moderate  blade  deflections  as  illustrated  in  Pig. 
1  These  equations  form  the  basis  of  an  implicit  flap-lag-torsional  undergoing  small  strains  and  moderate  deflections. 
Thus  the  equations  contain  geometrically  nonlinear  terms  in  the  structural,  inertia,  and  aerodynamic  operators  associ¬ 
ated  with  this  acroclastic  problem.  The  partial  dependence  in  the  equations  is  eliminated  by  using  a  Galerkm  type 
finite  cletnent  method  (19).  A  modal  coordinate  transformation,  using  six  rotating  coupled  modes,  is  performed  to 
reduce  the  number  of  degrees  of  freedom.  These  modes  are  calculated  at  a  fixed  value  of  collective  pitch  which 
depends  only  on  advance  ratio.  For  the  configurations  considered,  the  six  lowest  modes  arc  usually  the  first  three  dap, 
first  two  lead-lag.  and  the  fundamental  torsional  modes.  The  ordinary  differential  equations  are  solved  using  quasil- 
inearization  in  an  iterative  manner  to  obtain  the  periodic  equilibr  „ni  position  in  forward  flight,  for  a  propulsive  mm 
type  flight  condition. 

The  inertia  loads  arc  determined  by  using  D’Alembert's  principle.  An  implicit  formulation  for  the  aerodynamic 
loads  is  used  At  each  iteration  an  approximation  to  the  blade  response  is  produced.  This  response  is  then  used  to  gen¬ 
erate  numerical  values  of  the  modeling  quantities  needed  in  expressions  for  the  aerodynamic  loads  to  be  used  in  die 
next  iteration.  The  equations  are  linearized  by  writing  perturbation  equations  about  the  nonlinear  equilibrium  position. 
Stability  is  determined  by  using  Flcquet  theory. 

The  unsteady  aerodynamics  arc  finite-state  time-domain  aerodynamics  presented  in  Ref.  99  which  arc  an 
improved  version  of  the  theory  developed  by  Dinyavari  and  Friedmann  [100],  Stall  and  compressibility  effects  are 
neglected  The  acroclastic  model  is  combined  with  a  new  trim  procedure  (95,96)  which  provides  fully  coupled  simul¬ 
taneous  solution  to  both  the  acroclastic  response  and  the  trim  problem. 

The  rotor  dynamic  model  has  provisions  for  inclusion  of  HHC  pitch  changes  in  the  structural,  inertial,  and  aero¬ 
dynamic  portions  of  the  model  to  allow  HHC  to  be  input  to  the  model  as  4/rev.  sine  and  cosine  components  in  each  of 
the  collective,  longitudinal,  and  lateral  control  degrees  of  freedom.  The  general  HHC  input  may  be  expressed  as: 

%!  =(0ossino>HHV+0oc««WHHY] 

+  (0Cs  sin  “mt  V + ©cccos  WiikvIcos  \j / 

+  (0sssinG)HHV+0sccos<i)Hny]sin\j/  ^ 

where  0oc.  0os_.  ®cs.  0cc.  0ss>  a,,d  Ssc  are  independent  of  y.  The  hub  shears  to  be  minimized  arc  calculated  by 
integrating  the  inertia  and  aerodynamic  loads  due  to  a  given  blade  response  over  the  blade.  The  loads  due  to  the  four 
blades  are  then  combined  and  the  total  loads  are  transformed  to  the  non-rotating  coordinate  system.  Fourier  series 
represer'ations  of  these  loads  are  then  found  and  the  4/rcv.  components  extracted.  Alternately  to  calculate  blade  root 
loads  the  Fourier  analysis  is  done  in  the  rotating  system  at  the  blade  root. 

Operating  the  electro-hydraulic  actuators  needed  to  implement  HHC  will  require  power  from  the  helicopter 
powcrplant.  In  addition,  the  helicopter  rotor  may  require  more  or  less  power  than  at  the  baseline  condition  because  of 
the  additional  aerodynamic  loads  which  arc  imposed  on  it  by  the  HHC  inputs.  The  total  power  used  to  operate  the 
rotor  consists  of  the  power  needed  to  drive  the  rotor  and  the  power  needed  to  input  control  pitch  at  the  blade  root. 

The  various  contributions  to  the  power  required  to  actuate  the  blade  were  calculated  with  a  level  of  care  and 
detail  (96,97)  which  was  not  available  in  previous  studies. 

The  vast  majority  of  all  HHC  investigations  to  date  have  used  linear  optimal  control  sohitions  based  on  a  qua¬ 
dratic  cost  functional.  Therefore  this  approach  was  also  used  in  Refs.  95-97.  Minimum  variance  control  is  based  on 
the  minimization  of  a  cost  functional  which  is  the  expected  value  of  a  weighted  sum  of  the  mean  squares  of  the  control 
and  vibration  variables. 
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Minimum  variance  controllers  arc  obtained  by  minimization  of  the  cost  functional: 


j  =  n 


!zT(i)\v/z<i)  + 


0T(i)Wo0(i)  4  A0T(i)WA0AG(i) 


(10) 


Typically  Z,  0,  and  AO  consist  of  the  sine  and  cosine  components  of  the  N/rev.  vibrations  and  HI  1C  inputs.  The 
weightings  of  each  of  these  parameters  may  be  changed  to  make  it  more  or  less  important  than  the  other  components. 


The  minimum  variance  controllers  are  obtained  by  taking  the  partial  derivative  of  J  with  respect  to  G(i)  and  set 
ting  this  equal  to  zero: 


JLL 

<10(i) 


=  0 


(11) 


The  resulting  set  of  equations  may  lie  solved  for  the  optimal  I II  1C  input  *Q(i). 


The  form  of  the  resulting  algorithm  will  depend  on  whether  the  global  or  local  system  model  is  used  and  on 
whether  a  deterministic  or  cautious  controller  is  desired. 


The  global  model  of  the  helicopter  response  to  HHC  assumes  linearity  over  the  entire  range  of  control  applica 

non. 

Z(i  +  1)  =  Zo  +  T0(i)  (l2) 

The  vibration  vector  Z  at  step  i  +  I  is  equal  to  the  baseline  vibration  Zo  plus  the  product  of  the  transfer  matrix  T  and  the 
control  vector  0  at  step  i.  This  implies  that  T,  the  transfer  matrix  relating  HHC  inputs  to  vibration  outputs,  is  indepen 
dent  of  G(i). 


The  local  model  of  the  helicopter  response  to  HHC  is  a  linearization  of  the  response  about  the  response  to  the 
current  value  of  (he  controller  vector: 


Z(i  +  1 )  =  Z(i)  +  T(G(i  +  1 )  -  G(i)] 


(13) 


or: 


AZ(i  +  1)  =  TAQ(i+  I) 


(14) 


The  vibration  vector  Z  at  step  i  +  1  is  equal  to  the  vibration  vector  at  step  i  plus  the  product  of  the  transfer  matrix  and 
the  difference  in  the  control  vector  from  step  i  to  step  i  +  1.  Tins  allows  for  variation  of  the  transfer  matnx  T  with  input 


A  deterministic  and  cautious  minimum  variance  controller  can  be  programmed  into  two  algorithms,  one  for  the 
local  and  one  for  the  global  HHC  model  (90,92). 

A  few  selected  results  together  with  a  summary  of  the  most  important  conclusions  obtained  in  the  course  of  this 
study  (95,97)  arc  presented  below. 

An  interesting  test  of  the  ability  of  the  controllers  to  adapt  to  changing  flight  conditions  was  performed  by 
introducing  a  step  change  in  advance  ratio  from  n  =  0.30  to  p  =  0.35.  This  was  done  by  starting  with  a  converged 
optimal  solution  and  its  response  at  p  =  0.30,  changing  the  propulsive  trim  values  to  those  for  p  =  0.35,  and  proceeding 
with  iterative  control  calculations  and  quasilinearization  solution.  Results  for  a  four-bladed  soft-in-plane  hingeless 
rotor  are  shown  in  Figs.  1 1  and  12. 

The  iteration  history  of  4/rcv.  hub  shears  for  the  deterministic  and  cautious  local  control  arc  shown  in  Fig.  1 1 
and  the  iteration  history  for  the  global  control  arc  shown  in  Fig.  12.  Application  of  HHC  eliminated  essentially  all  the 
4/rev.  hub  shears  within  five  iterations.  With  the  local  controller,  the  vertical  shears  rise  slowly  after  the  second  itera¬ 
tion  indicating  that  the  transfer  matrix  identification  has  not  been  ideal.  For  lack  of  space  the  control  inputs  are  not 
shown  here.  However  examination  of  these  inputs  (95)  reveals  that  when  the  step  change  in  advance  ratio  was  applied 
to  the  local  controller  there  were  large  oscillations  in  the  calculated  control  inputs  and  the  resultant  hub  shears.  On  the 
other  hand  when  a  step  change  in  p  was  applied  to  the  global  controller,  the  oscillatory  behavior  in  the  required  control 
input  was  reduced  substantially  and  the  controller  moved  fairly  smoothly  from  an  initial  large  increase  in  shears  toward 
a  minimum.  A  comparison  of  the  three  shear  components  and  their  baseline  values  for  the  local  and  global  controllers 
is  given  in  Fig.  13.  It  is  evident  that  the  global  controller  has  been  more  successful  in  reducing  shears. 

Another  interesting  result  is  associated  with  the  instantaneous  control  power  as  a  function  of  azimuth  which 
was  evaluated  for  both  the  hingeless  rotor  and  a  comparable  articulated  rotor,  at  an  advance  ratio  of  p  =  0.30.  Figures 
14  and  15  show  the  variation  in  instantaneous  trim  power,  HHC  power,  and  total  control  power  respectively  for  die 
articulated  and  hingeless  blades. 


For  the  bngeless  blade  die  HHC  power  contributes  relatively  much  more  to  the  total  control  power  as  can  be 
observed  from  the  nature  of  diesc  curves  which  tend  to  follow  one  another  closely  in  Fig.  15.  As  is  evident  from  com¬ 
paring  the  vertical  scales,  the  maximum  power  excursions  arc  on  the  order  of  five  times  larger  for  the  hingeless  blade 
than  for  the  articulated  blade  and  in  addition  are  much  sharper.  The  peaks  arc  biased  toward  positive  power  values  so 
that  the  total  control  power  over  one  revolution  is  much  higher  for  the  hingeless  than  die  articulated  blade. 


The  most  important  conclusions  obtained  from  these  studies  were: 


Overall  blade  response  was  compared  for  a  baseline  no  HHC  condition  and  the  optima)  reduced  vibration  con¬ 
dition  Principal  differences  were  in  the  torsional  and  flap  response.  Overall  response  magnitudes  changed  lit¬ 
tle  but  large  4/rcv.  components  were  introduced  to  modify  the  airloads  and  cancel  out  vibrations. 
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2  A  global  controller  was  used  in  comparison  of  the  effects  cf  I UIC  on  roughly  equivalent  articulated  and  lunge 
less  rotors  Shears  were  successfully  reduced  for  both  rotors  by  a  Hi  1C  algorithm  winch  tinninu/ed  just  huh 
shears.  As  shears  were  reduced  there  were  large  increases  in  hub  moments  for  the  hingeles  rotor  but  only 
moderate  increases  for  the  articulated  rotor.  Much  larger  I II  1C  angles  were  required  to  reduce  shears  for  the 
hingelcss  rotor. 

1  Attempts  to  reduce  both  hub  shears  and  moments  were  not  very  successful  for  either  blade  With  the  ar'icu- 
latcd  blade  moments  were  already  low  and  were  only  decreased  slightly  while  shears  increased  slightly.  I  or 
the  hingeless  blade  the  large  moments  decreased  greatly  but  only  at  Ibc  expense  of  poor  shear  reduction  'I  iicse 
results  indicate  that  vibration  reductions  with  I U  1C  may  be  more  diflictilt  in  hingeless  rotor  configurations  than 
for  articulated  configurations. 

4  Application  of  HHC  to  the  lungeless  rotor  lead  to  an  increase  in  required  power  of  1.44%  while  for  the  articu¬ 
lated  rotor  this  increase  was  only  0.18%.  The  required  power  increase  for  the  lungeless  rotor  was  somewhat 
mitigated  by  a  0.2%  increase  in  rotor  thrust. 

5  Overall  blade  acroclastic  stability  margins  were  not  significantly  degraded  by  application  of  HHC  for  culler  the 
articulated  or  hingelcss  blade. 


It  is  interesting  to  note  that  two,  somewhat  similar  studies,  on  the  application  of  higher  harmonic  control  to 
hingelcss  rotor  systems  were  also  recently  completed  by  Nguyen  and  Chopra  [101,102).  These  studies  were  based  on 
an  advanced  HHC  acroclastic  simulation  capability.  The  analysis  utilities  finite  element  in  both  space  and  time.  A 
nonlinear  time  domain  unsteady  aerodynamic  model  is  used  to  obtain  the  air  loads,  and  the  rotor  induced  tnllow  is  cal¬ 
culated  using  a  free  wake  model.  The  vehicle  trim  is  also  obtained  from  a  fully  coupled  trim/rotor  acroclastic  analysis. 
Thus  this  simulation  capability  has  many  common  features  with  Refs.  95-97.  except  that  the  aerodynamic  model  is 
more  refined  when  compared  to  that  used  by  Robinson  and  Friedmann  [95-97].  The  higher  harmonic  control  algorithm 
is  also  very  similar  to  that  employed  in  Refs.  95-97.  This  acroclastic  simulation  capability  was  validated  by  comparing 
it  with  the  wind  tunnel  tests  conducted  on  a  one-sixth  dynamically  sealed  thrce-bladcd  articulated  rotor  model  tested  by 
Boeing  Helicopter  Co.  [82]  up  to  high  advance  ratios,  p  =  0.40,  and  fairly  good  correlation  with  the  expenmeiiMl  data 
was  obtained. 

The  conclusions  obtained  from  this  study  arc  in  agreement  with  the  maj'ority  of  the  conclusions  noted  in  Refs. 
95-97  Thus  the  performance  of  die  global  controller  was  superior  to  the  local  controller.  When  applying  HHC  on 
hmgeless  rotors  larger  control  angles  were  required  than  for  articulated  rotor  and  a  substantial  increase  in  the  HHC 
power  requirements  was  also  noted.  Furthermore  HHC  effects  on  rotor  performance  were  found  to  be  small.  An 
interesting  effect  noted  was  the  increase  in  stall  areas  over  the  rotor  disk,  when  the  local  HHC  controller  was  at  the 
edge  of  the  flight  boundary. 

Fairly  extensive  tests  on  a  four  bladed  hingeless  rotor,  with  a  minimum  variance  controller  w-cre  also  described 
in  detail  in  a  recent  paper  by  Lehmann  and  Kube  [103].  An  adaptive  local  controller  was  found  to  perform  quite  well 
over  the  whole  envelope  tested.  Considerable  detail  on  the  actual  hardware,  and  digital  control  implementation  of 
HHC  for  the  lest  configuration  is  presented  in  the  paper.  Effective  vibration  reduction  for  the  lungeless  rotor  was 
obtained  with  a  gain  adjustment  (adaptation)  algorithm.  The  power  requirement  for  HHC  during  the  test  was  not  dis¬ 
cussed  in  this  study. 

A  very  comprehensive  wind  tunnel  lest  program  on  a  three  bladed  articulated  model  rotor  (CH-47D)  was  also 
conducted  by  Shaw  et  al.  [82],  It  was  found  dial  HHC  is  highly  effective  on  articulated  three  bladed  rotors.  Harmonic 
hub  load  response  to  HHC  was  found  to  be  linear  up  to  3  degrees  and  it  was  insensitive  to  flight  condition.  Further¬ 
more  a  fixed-parameter  control  law  was  found  to  be  fully  effective  for  vibration  reduction.  Figure  16,  taken  from  Ref. 
82,  shows  that  closed  loop  HHC  was  extremely  effective  in  suppressing  vibratory  hub  forces.  The  fixed-gain  con¬ 
troller,  when  configured  to  regulate  3/rev  vertical  force  and  2/rev  and  4/rev  inplane  shears,  suppressed  all  three  of  these 
components  simultaneously  by  90  percent  in  almost  all  of  the  trimmed  and  quasi-stcady  maneuvering  test  envelope. 

From  the  discussion  presented  in  this  section  it  is  evident  that  comprehensive  aeroelastic  simulations  of  HHC 
have  a  very  useful  role  in  this  field.  In  addition  to  being  much  more  cost  effective  than  wind  tunnel  tests  and  flight 
tests,  they  provide  one  with  an  ideal  basis  for  planning  such  tests.  Furthermore  experience  with  flight  tests  [67]  indi¬ 
cates  that  scale  model  tests  do  not  always  correlate  well  with  wind  tunnel  tests.  Thus  correlation  between  aeroelastic 
simulation  programs  of  HHC  which  correlate  well  with  wind  tunnel  tests  provide  one  with  a  tool  suitable  for  simulat¬ 
ing  flight  tests. 

The  evidence  available  from  recent  studies  seems  to  imply  that  the  practical  implementation  of  HHC  on  hinge¬ 
less  rotor  systems  could  be  more  difficult,  and  less  effective,  than  similar  implementation  of  HHC  on  articulated  rotors. 
Control  angles  needed  for  HHC  on  hingeless  rotors  are  significantly  larger  and  power  requirements  could  be  also 
between  2-5  times  larger.  Simultaneous  reduction  of  both  hub  shears  and  moments  on  hingelcss  rotors  is  also  more 
difficult,  and  rotating  blade  root  loads  can  be  also  substantially  larger.  One  possible  explanation  for  these  observations 
is  the  strong  physical  coupling  between  the  bending  and  torsional  degrees  of  freedom  which  exists  in  most  hingelcss 
rotor  blades  This  can  also  cause  the  increase  in  the  power  required  for  implementing  HHC  on  such  rotors.  Since  die 
higher  harmonic  movement  of  the  complete  blade  during  its  HHC,  introduces  more  complex  motion  due  to  the  cou¬ 
pling  effect  and  thus  requires  more  power. 

Finally,  it  appears  that  adaptive  controller  may  not  be  always  required  in  order  to  implement  an  effct.ive  HHC 
i  vibration  reduction. 


6.  CONCLUDING  REMARKS 

This  paper  provides  a  detailed  discussion  of  four  important  topics  in  helicopter  rotor  dynamics  and  acroclasu- 
city  Hopefully  this  discussion  will  provide  an  improved  fundamental  understanding  of  the  current  state  of  the  an. 
Tf'us  future  research  on  these  topics  can  be  focused  on  problems  which  remain  to  be  solved  instead  of  producing  mar¬ 
ginal  improvements  on  problems  which  are  already  well  understood. 
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Tor  acroclastic  and  acromcchamcal  stability  problems  incorporation  of  geometrically  nonlinear  terms  due  to 
moderate  dellcctions  has  become  fairly  routine.  Both  moderate  deflection  beam  theories  and  large  deflection  theories 
are  available,  and  both  can  serve  as  a  basis  for  acroclastic  stability  and  response  analysis.  It  is  noted  that  improvement 
of  the  unsteady  aerodynamics  in  rotary  wing  aeroelastic  analyses  is  much  more  important  than  the  incorporation  of 
higher  order  geometrically  nonlinear  terms.  Excessive  preoccupation  with  such  higher  order  terms  provides  only 
diminishing  benefits. 

A  number  of  recent  composite  beam  theories,  suitable  for  modeling  composite  blades  have  been  developed 
however  these  have  yet  to  be  implemented  in  a  comprehensive  acroclastic  stability  or  response  analysts.  The 
significant  potential  for  acroclastic  tailoring  inherent  m  composite  rotor  blades  remains  to  be  exploited. 

The  body  of  research  available  on  coupled  rotor/fuselagc  acromechanical  problems  has  reached  a  remarkable 
level  of  maturity  in  a  fairly  short  period  of  lime.  Reliable  models  for  this  class  of  problems  should  contain  coupled 
flap  lag  torsional  blade  models  combined  with  a  fuselage  which  has  pitch  and  roll  as  well  as  two  translational  degrees 
of  freedom  Simpler  models  can  easily  lead  to  inaccurate  results.  Air  resonance  in  hingcless  rotors  can  be  actively 
controlled  throughout  the  entire  flight  envelope  using  a  simple  control  system  designed  using  modern  control  system 
synthesis  techniques. 

The  state  of  the  art  in  applying  HHC  to  articulated  rotors  appears  to  be  quite  promising.  The  application  of 
HHC  to  hingelcss  .otors  could  have  practical  implementation  problem.  Aeroelastic  simulations  of  HHC  have  a  very 
useful  role  in  providing  an  improved  fundamental  understanding  on  the  implementation  of  HHC  on  different  types  of 
rotor  systems.  A  substantial  amount  of  additional  research  in  this  topic  is  required  so  as  to  guarantee  its  implementa¬ 
tion  on  the  next  generation  of  rotorcraft. 
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Table  I 


Frequency  Increase  due  to  ft 

Frequency  Decrease  due  to 

Timoshenko  Beam  Effects 

Inplane  Modes 

Inplane  Modes 

Mode  1 

Mode  2 

Mode  3 

Mode  1 

Mode  2 

Mode  3 

3.5% 

2.6% 

1.0% 

-1.6% 

-9.7% 

-19.2% 

Out  of  Plane  Modes 

Out  of  Plane  Modes 

71.2% 

12.8% 

4.6% 

’-0.4% 

-2.5% 

-5.5% 

Figure  la:  Geometry  of  (lie  Blade  Elastic  Axis  Before  and  Alter  Deformation. 


-17- 


rOSITlON  Of  8t*0C  CROSS  SCCTION 
eUOfU  OtfOflMATIOM 


Figure  lb:  Blade  Cross  Sccnon.il  Gcomciry. 
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COLLECTIVE  PITCH  ANGLE  (OEC) 


Figure  2:  Lead-lag  Damping  Versus  Blade  Pilch  Angle,  Soil 
Pilch  Flexure,  Zero  Droop.  Zero  Preconc  (Ref.  30). 
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Figure  3:  Box  Beam  for  Study  of  Timoshenko  Beam  Effects. 


Figure  5:  Coupled  Rotor/Fusclagc  Model 


Figure  4:  Pcak-to-Pcak  Vertical  Hub  Shear.  Ilingcless  Com¬ 
posite  Blade  (Ref.  48). 
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Figure  6:  Offset  Hinged  Spring  Restrained  Model  for  llingelcss  Blade 
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Figure  7:  Open  Loop  Lead-lag  Regressing  Damping  of  the 
Nominal  Configuration  With  and  Without  Dynamic  Inflow. 
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Figure  8:  Open  Loop  Lead-lag  Regressing  Damping  of  die 
Nominal  Configuration  With  and  Without  Blade  Torsional 
Flexibility. 


Figure  9:  Open  Loop  Lead-Lag  Regressing  Mode  Damping  at  Figure  10:  Closed  Loop  Lead-Lag  Regressing  Mode  Damping 
Various  Fuselage  Weights  and  Advance  Ratios.  Using  die  Active  Controller.  ' 
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Figure  11.  Iteration  History  of  Hub  Shears  for  Deterministic 
versus  Cautious  Local  Control  Optimizauon  at  p  =  0.30.  then 
Step  change  from  n  =  0.30  to  ft  =  0.35. 
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Figure  12:  Iteration  History  of  Hub  Shears  for  Dcicrmimsuc 
versus  Cautious  Global  Control  Optimization  at  p  =  0.30.  then 
Step  change  from  p  =  0.30  to  p  =  0.3S. 
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Figure  13:  Baseline  Shears  for  p  =  0.35 
lions  After  Step  Change  from  Controllers. 


and  p  =  0.30  to  p  =  0.35,  Local 


and  Global  Shears  Five  Itcra- 


power  (<rocl»on  of  lofol) 


Figure  14:  Vanation  of  Trim  Power.  HHC  Power  and  Total 
Control  Power  with  Azimuth  Articulated  Blade,  Advance  Ratio 
P  =  0.30;  5%  Blade  Root  Offset,  Cautious  Global  Control. 
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Figure  15:  Vanation  of  Trim  Power.  HHC  Power  and  Total 
Control  Power  with  Azimuth,  llingclcss  Blade,  Advance  Ratio 
p  =  0.30, 5%  Blade  Root  Offset,  Cautious  Global  Control. 
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Figure  16:  Measured  HHC  Effectiveness  in  Trimmed  Level  Right,  (Ref.  82). 
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PHYSICAL  ASPECTS  OF  RO'J'OH-BODY  COUPLING 


H .  C.  Curtiss,  Jr. 

Princeton  University 

The  influence  of  various  modelling  assumptions  on  the 
dynamic  stability  and  response  characteristics  of  helicopters 
will  be  discussed  using  a  linearized  theory  with  body  and  rotor 
degrees  of  freedom.  Emphasis  will  be  placed  on  physical 
interpretation  of  the  primary  sources  of  rotor-body  coupling  and 
their  influence  on  the  response  and  stability  characteristics. 
Prediction  of  the  time  and  frequency  domain  response 
characteristic  of  holicopters  to  control  inputs  will  be  discussed 
and  the  effect  of  various  approximations  on  the  correlation  of 
theory  with  experiment  described. 


Analysis  of  Sectional  Properties  of  Curved, 
Twisted,  Nonhomogeneous,  Anisotropic  Beams 

Marco  Borri 

Dipartimento  di  Ingegneria  Aerospaziale 
Politecnico  di  Milano 
Milano,  Italy 

Dewey  H.  Hodges,  Ali  R.  Atilgan,  and  Victor  L.  Berdichevsky 
School  of  Aerospace  Engineering 
Georgia  Institute  of  Technology 
Atlanta,  Georgia 

The  nonlinear  theories  of  Borri  and  Mantegazza  (1985)  and  Hodges  (1989)  express  the  strain  field 
in  terms  of  so-called  intrinsic  strain  measures  which  have  the  unique  property  that  the  derivatives 
of  the  strain  energy  with  respect  to  these  measures  are  the  actual  section  forces  and  moments  in  the 
deformed  beam  basis.  Therefore,  the  constitutive  relations  can  be  expressed  in  terms  of  material 
and  geometric  stiffness  constants  for  a  cross  section.  An  analysis  of  the  sectional  properties  of 
nonhomogeneous,  anisotropic  beams  is  presented  which  takes  initial  curvature  and  twist  into 
account.  All  possible  types  of  couplings  between  extension^,  shear,  bending,  and  torsional 
deformation  are  included.  The  analysis  is  an  extension  of  the  earner  work  by  Giavotto  et  al.  (1983) 
and  leads  to  a  two-dimensional  finite  element  formulation  over  the  section  under  consideration.  The 
method  will  be  compared  with  other  methods  in  the  literature  including  the  work  of  Soviet 
mechanicians. 


Nonlinear  Analysis  of  Curved,  Twisted, 
Nonhomogeneous,  Anisotropic  Beams 


Dewey  H.  Hodges  and  Ali  R.  Atiigan 
School  of  Aerospace  Engineering 
Georgia  Institute  of  Technology 
Atlanta,  Georgia 

Tremendous  progress  has  been  made  in  the  last  few  years  toward  simplifying  the  analysis  of 
beams  which  are  undergoing  large  deformation.  In  previous  years,  analysts  became  accustomed  to 
writing  beam  equations  explicitly  in  displacement  form  and  to  tedious  ordering  schemes  which 
could  not  be  applied  rigorously.  For  some  applications  the  size  and  complexity  of  the  equations 
become  overwhelming,  leading  to  the  use  of  symbolic  manipulation  just  to  derive  the  equations.  In 
what  will  be  presented,  the  exact  equations  for  a  beam  undergoing  arbitrarily  large  deformation  will 
be  written  in  only  a  few  lines.  l.'Ic  ‘cover,  these  equations  lead  to  a  finite  element  method  which  is 
not  only  elegant,  but  very  simple  to  program.  The  analysis  methodology  leads  to  a  natural 
breakdown  of  the  analysis  into  th  ee  parts:  (1)  exact  intrinsic  equations  of  equilibrium  in  weakest 
possible  form;  (2)  exact  kinematical  equations  in  weakest  possible  form;  (3)  approximate 
constitutive  law  in  terms  of  sectional  elastic  constants.  An  overview  of  the  process  is  shown  below 
in  Figure  1. 


Figure  1:  Overview  of  Anisotropic  Beam  Analysis 


Investigation  of  the  Use  of  Extension-Twist  Coupling  in 
Composite  Rotor  Blades  for  Application  to  Tiltrotor  Aircraft 


Renee  C.  Lake 
U.S.  Army  Aerostructures  Directorate 
NASA  Langley  Research  Center 
Hampton  VA  23665 

Abstract 


Tiltrotor  aircraft  are  designed  to  operate  in  both 
helicopter  and  airplane  modes  of  flight.  This  operational 
flexibility  results  in  several  conflicting  design 
requirements.  One  such  design  requirement,  which  has 
significant  effects  on  aerodynamic  performance,  is  blade 
twist.  Typically,  the  twist  employed  is  a  compromise 
between  that  required  in  the  two  different  modes  of  flight. 
Performance  could  be  improved  if  it  were  possible  to  vary 
the  blade  twist  between  the  two  flight  modes.  Tiltrotor 
aircraft  typically  vary  rotor  speed  by  about  20  percent 
between  the  two  flight  modes.  This  change  in  rotor  speed 
induces  a  rather  substantial  change  in  the  centrifugal  force 
along  the  blade,  which  could  be  used  to  passively  change 
the  twist  of  an  extension-twist-couplec-  r,.omposite  blade. 
An  in-house  investigation  was  initiated  several  years  ago 
within  the  Army  Aerostructures  Directorate  at  NASA 
Langley  Research  Center  with  the  overall  objective  of 
demonstrating  the  increased  aerodynamic  performance 
which  can  be  realized  in  tiltrotor  blades  that  have  been 
elastically-tailored  to  exhibit  extension-twist  structural 
coupling.  Analytical  and  experimental  studies  have  been 
designed  to  demonstrate  the  improvements  in  tiltrotor 
blade  performance,  and  ultimately  include  whirl  tests  and 
wind-tunnel  tests  conducted  on  a  model-scale  blade  in 
addition  to  advanced  finite-element  analyses.  However,  the 
investigation  of  research  areas  which  are  fundamental  to 
the  development  of  elastically-tailored  composite  tiltrotor 
blades  must  first  be  addressed.  For  example,  the  feasibility 
of  producing  the  desired  twist  deformation  within  the 
operational  and  material  design  limits  for  a  tiltrotor  blade 
is  a  principal  driver  in  the  effective  development  and 


evolution  of  extension-twist-coupled  tiltrotor  blades.  In 
addition,  the  structural  dynamic  characteristics  of 
extension-twist-coupled  rotor  blades  must  be  assessed,  as 
the  elastic  couplings  can  have  a  significant  influence  on  the 
dynamic  characteristics,  such  as  frequency  and  mode  shape. 
Also,  the  bending  and  twisting  deformations  of  the  primary 
load  carrying  members  of  a  rotor  blade  are  usually 
accompanied  by  additional  deformations  associated  with 
cross-sectional  warping.  Therefore,  the  significance  of  such 
warping  effects  for  extension-twist-coupled  rotor  blades 
must  also  be  determined.  While  each  of  these  research 
areas  serves  to  drive  the  overall  design  of  a  rotor  blade  and 
is  being  addressed,  this  paper  focuses  on  the  investigation 
of  structural  dynamic  behavior  of  extension-twist-coupled 
composite  rotor  blades. 

A  study  was  initiated  in  which  dynamic  results  from 
shell-finite-element  models  were  correlated  with 
experimental  results  obtained  from  a  vibration  test  of 
extension-twist  coupled  composite  tubular  specimens.  The 
specimens  all  had  non-circular  cross  sections  and  were 
therefore  subject  to  warping  deformations.  A  potentially 
significant  effect  on  the  dynamic  behavior  of  the  specimens 
could  occur  because  the  presence  of  extension-twist 
coupling  in  a  laminate  acts  to  "magnify"  the  warping  effect. 
Such  "warping-prone"  designs  are  more  closely 
representative  of  true  rotor  blade  structures  and  therefore 
require  a  thorough  investigation  of  the  dynamic 
characteristics.  Three  thin-walled  cross-sectional 
geometries  were  studied:  square,  D-shape,  and  flattened 
ellipse  (see  figure  1).  The  test  specimens  were  fabricated 
utilizing  an  expandable  rubber  core  and  a  segmented 
aluminum  female  tool  configured  to  the  desired  final 
specimen  geometry.  The  overall  size  and  dimensions  of 
each  specimen  type  are  also  indicated  in  figure  1.  A  free- 
free  boundary  condition  was  simulated  in  the  test  by 
suspending  the  specimens  from  a  test  support  structure 
with  rubber  tubing.  A  vibration  survey  was  conducted  on 
each  specimen  to  determine  structural  modes  and 
frequencies.  A  calibrated  impact  hammer  provided  the 
input  excitation  for  each  specimen  such  that  a  "roving 
hammer"  concept  was  employed.  Three  single-axis 
accelerometer  transducers  were  block-mounted  at  an  offset 
on  one  end  of  the  specimens  in  a  tri-axial  configuration. 


2 


Data  acquisition  was  performed  using  an  HP-3565S 
computer  system  using  HP-VISTA  software.  Subsequent 
curvefitting  and  parameter  extraction  was  facilitated  with 
the  use  of  SDRC-Tdas  software.  These  results  were 
compared  with  those  obtained  from  the  MSC-NASTRAN 
shell-finite-element  analyses  for  each  specimen  type.  In 
addition  to  the  global  mode  shapes  of  the  shell  structures, 
numerous  shell  breathing  modes  were  obtained.  However, 
these  breathing  modes  are  not  of  current  interest  in  this 
study  and  will  not  be  addressed  in  the  discussion  of  results. 

The  preliminary  comparison  of  experimental  and 
analytical  data  has  been  conducted  and  is  shown  in  figures 
2  through  4.  The  data  correlation  for  the  elliptical 
specimens  (figure  2)  showed  the  best  overall  agreement 
with  experiment,  yielding  results  within  5%  for  the  first 
five  global  modes  (first  vertical  bending,  first  lateral 
bending,  second  vertical  bending,  second  lateral  bending, 
first  torsion),  and  within  13%  for  the  sixth  mode  (third 
vertical  bending).  The  comparison  of  results  for  the  D- 
shape  specimen  (figure  3)  produced  agreement  generally 
within  6.5%  for  the  five  identified  modes,  with  the  third 
global  mode  (second  vertical  bending)  producing  a  "worst 
case"  agreement,  dropping  off  to  14%.  Frequency  and  mode 
shape  results  for  the  square  cross-section  specimen  (see 
figure  4)  showed  agreement  within  10%  for  the  first  four 
global  modes.  In  general,  curvefitting  of  data  was 
performed  in  the  0-2000  Hz  range,  as  frequencies  above 
this  range  generally  showed  diminished  phase  and 
coherence  characteristics. 

In  parallel  with  the  aforementioned  studies  on  the 
tubular  specimens,  studies  have  also  been  initiated  on  an 
extension-twist-coupled  model  rotor  blade.  The  model 
blade,  which  was  manufactured  from  existing  blade  molds 
for  an  isotropic  rotor  blade,  was  designed  with  a  view 
towards  establishing  a  preliminary  "proof  of  concept"  for 
extension-twist-coupled  rotor  blades.  Whirl  tests  of  the 
model  blade  will  provide  the  basic  objective  of  measuring 
blade  twist  as  a  function  of  rotor  speed.  This  paper  will, 
however,  focus  on  the  investigation  of  the  structural 
dynamic  characteristics  of  the  model  blade.  The  resultant 
shell-finite-element  model  of  the  model  rotor  blade  is 
shown  in  figure  5.  Results  from  a  normal  modes  analysis  of 
the  rotor  blade  shell-finite-element  model  have  been 
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compared  with  those  obtained  from  a  ground  vibration  test 
conducted  on  a  fabricated  model  blade,  for  both  free-free 
and  cantilevered  boundary  conditions,  and  are  shown  in 
figure  6.  In  general,  the  mode  shape  and  frequency 
correlation  shows  agreement  within  10%  for  the  first  five 
free-free  modes.  Similar  trends  are  observed  for  the 
cantilevered  blade,  with  all  identified  frequencies  generally 
within  this  10%  agreement  range,  with  the  exception  of  the 
first  torsion  and  fourth  flap  bending  modes,  which  were 
found  to  agree  within  approximately  13%.  This  comparison 
of  experimental  and  analytical  results  serves  to  establish  a 
foundation  for  accurately  modeling  and  predicting  the 
structural  dynamic  behavior  of  extension-twist-coupled 
rotor  blades. 
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Abstract 

In  the  aeroelastic  modeling  and  analysis  of  rotorcraft,  there  are 
several  key  steps  that  must  be  accomplished.  First,  the  various  equations 
for  the  pieces  of  the  model  (i.e.,  the  structure,  lift,  induced  flow,  and 
control  system)  must  be  assembled  into  a  set  of  equations  that  can  be  solved 
by  some  solution  strategy.  Second,  one  must  find  the  trim  settings  that  will 
result  in  a  periodic  solution  at  a  given  flight  condition.  Third,  one  must 
perturb  this  flight  condition  to  obtain  perturbation  equations.  Fourth, 
these  perturbation  equations  must  be  solved  in  order  to  find  frequencies, 
modes,  and  stability  boundaries;  and  then  control  systems  can  be  designed. 
As  it  turns  out,  the  determination  of  trim  and  the  development  of 
perturbation  equations  can  be  greatly  expedited  if  the  assembled  rotorcraft 
model  can  be  written  in  terms  of  a  finite  number  of  state  variables. 
Unfortunately,  most  lift  and  inflow  models  do  not  lend  themselves  to  state 
representations  and,  thus,  are  difficult  to  use  in  aeroelasticity  analysis 
beyond  simple  time-marching. 


In  this  paper,  we  will  review  the  ARO-sponsored  work  in  which  we  are 
attempting  to  model  the  entire  rotor  system  with  state-variable  models.  The 
structure  is  modeled  in  terms  of  generalized  Rayleigh-Ritz  coordinates,  the 
lift  is  modeled  in  terms  of  the  modified  ONERA  dynamic-stall  state  variables 
(circulations  and  moment  coefficients),  the  wake  is  modeled  in  terms  of  a 
finite-state  induced  flow  model,  and  trim  is  found  from  a  state-space 
auto-pilot.  This  combination  is  also  interesting  from  the  point  of  view  that 
an  unsteady  wake  model  implicitly  gives  some  unsteady  aerodynamic  effects 
which  must  then  be  expunged  from  the  lift  model  so  as  not  to  be  counted 
twice.  In  our  case,  this  is  quite  simple  to.  do  because  the  ONERA  model  is 
already  segregated  into  the  part  due  to  linear,  unsteady  aerodynamics  and  the 
part  due  to  stall.  Thus,  we  simply  eliminate  the  linear  state  variables  from 
the  lift  model.  Further,  we  treat  the  drag  in  a  quasi-steady  manner.  This 
leaves  four  aerodynamic  states  per  aerodynamic  control  point. 

Thus,  for  a  typical  computation,  we  would  have  6  structural  modes  (3 
flap,  2  lag,  1  torsion)  which  is  12  states,  5  aerodynamic  modes  which  is  20 
states,  four  harmonics  of  inflow  which  is  15  states,  and  three  controllers 
which  is  6  states.  Therefore,  there  are  53  total  states  for  time-marching 
and  47  states  for  stability  (the  trim  controller  is  turned  off  for 
perturbation  dynamics).  This  is  a  reasonable  number  of  states  with  which  to 
do  meaningful  Floquet  theory  and  control  system  design.  We  are  correlating 
our  results  with  the  SA349  Gazelle  measurements  and  will  show  some  of  these 
at  the  workshop. 
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A  variety  of  applications  for  spatial  and  temporal  finite  elements  is  described.  A  weak  Hamiltonian 
formulation  is  first  presented  for  the  dynamics  of  rigid  bodies.  With  this  method,  one  can  either 
time  march  or  solve  periodic  two-point-boundary-value  problems.  All  boundary  conditions  are  of 
the  “natural”  type,  and  one  can  use  the  crudest  possible  shape  functions  and  still  achieve 
remarkable  accuracy.  The  method  is  extended  to  include  formulation  and  solution  of  optimal 
control  problems  in  terms  of  states,  co-states,  and  control  variables.  Results  are  presented  for  a 
simple  trajectory  optimization  problem.  The  method  is  further  extended  for  determining  the 
nonlinear  static  response  of  beams  in  terms  of  displacement  variables,  generalized  strain,  and 
section  stress  resultants.  A  separate  two-dimensional  finite  element  code  leads  to  characterization  of 
the  constitutive  law  in  terms  of  21  elastic  constants  for  nonhomogeneous,  anisotropic  beams.  For 
the  global  deformation  of  a  beam,  the  entire  formulation  is  so  compact  that  it  can  be  written  in  a 
few  lines  -  even  for  arbitrarily  large  displacement  and  rotation.  Results  are  presented  that  show  the 
convergence  properties  of  this  mixed  method.  The  formulation  has  also  been  developed  for 
nonlinear  statics  and  dynamics  of  moving  beams  although  no  result.,  are  yet  available.  Future  work 
includes  space-time  response  of  moving  beams,  deployment  problems,  and  disturbance 
propagation.  Development  of  a  similar  approach  for  plate  problems  appears  to  be  feasible.  The 
potential  impact  of  this  methodology  on  rotorcraft  dynamics  problems  is  also  discussed. 


The  Rotorcraft  Center  of  Excellence  at  Rensselaer 


R.  G.  Loewy 


The  U.  S.  Army  Research  Office  Center  of  Excellence  at  Rensselaer, 
established  in  1982  and  supported  under  two  five-year  contracts,  is  in 
its  eighth  year.  This  ARO-COE,  known  as  the  Rensselaer  Rotorcraft 
Technology  Center  (RRTC) ,  continues  under  its  originally  defined 
three-part  concept:  first,  to  provide  advanced  and  comprehensive 
educational  programs  which  will  attract  some  of  the  best  and  brightest 
young  minds  to  the  study  of  rotorcraft  technology  and  to  up-date 
professionals  in  the  field,  in  programs  specialized  to  continuing 
education;  second,  to  perform  research  which  is  at  the  cutting  edge  of 
advances  in  knowledge  and  understanding  of  the  phenomena  pacing  new 
rotorcraft  development  and  to  transfer  these  advances  expeditiously  to  U. 
S.  developers;  and  third,  to  develop  new  experimental  and  computational 
resources  in  rotorcraft  technology  for  the  use  and  benefit  of  all 
appropriate  components  of  the  U.S.  rotorcraft  community. 

At  this  writing,  seven  faculty/staff  and  ten  graduate  students  are 
sponsored  in  their  pursuit  of  the  Center's  goals  under  ARO  auspices. 
Another  eight  faculty  members  and  twelve  graduate  students  are  active ,  to 
varying  degrees,  in  the  activities  of  the  Center.  Some  of  the  latter 
group  are  conducting  research  under  industry  contracts,  others  under 
contracts  or  grants  won  from  government  agencies  in  addition  to  the 
ARO-COE  funding,  and  still  others  with  Rensselaer  sponsorship. 

Areas  of  research  concentration  with  ARO-COE  support  include  innovative 
composite  structural  design,  analysis  and  fabrication;  composite  drive 
shaft  system  design;  non-linear  structural  dynamics  and  aeroe] asticity  - 
particularly  aimed  at  violent  maneuvers;  blade -vortex  interaction 
phenomena;  rotor  shaft  flexibility  influences  on  ground/air  resonance 
instabilities;  and  improving  means  of  mode  identification  and 
modification  in  fuselage  ground  shake  tests.  Similar  areas  of 
concentration  with  other  than  ARO-COE  funding  include  continuing 
transonic  wind  tunnel  tests  of  passive  wave  drag  reduction  mechanisms  in 
unsteady  flow,  developing  advanced  CAD/CAM-FEM  linkages  using  expert 
systems  techniques,  conducting  shake  tests  and  correlative  analyses  for 
the  natural  modes  and  frequencies  of  full  scale  helicopter  fuselage 
components,  completing  development  of  a  model  rotor  impedance  testing 
rig,  establishing  a  major  laboratory  for  symbollic  computation,  and 
developing  techniques  for  conducting  ground/air  resonance  stability 
flight  testing  using  radio-controlled  scale  models. 

Clearly,  attempts  to  make  the  Center  self-supporting,  ie  independent  of 
ARO-COE  continuation  beyond  July  1992,  the  end  of  the  present  contract, 
are  receiving  determined  attention.  Judging  the  RRTC  by  such  standard 
measures  as  (a)  the  quality  and  numbers  of  graduate  students  who  have  and 
are  entering  the  program,  (b)  the  number  and  ranking  of  positions 
graduates  of  the  program  have  and  are  to  continuing  to  enter  in  the 
rotorcraft  community,  and  (c)  publications  in  professional  journals  and 
presentations  at  national  and  international  meetings,  we  believe  the 
ARO-COE  program  can  be  pronounced  a  success  of  significant  proportions, 
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ABSTRACT 


Rotorcraft  designs  are  evolving  rapidly  towards  the  high-speed  end.  These  de¬ 
signs  challenge  the  conventional  barriers  of  rotorcraft  capability.  The  analytical 
support  and  evaluation  of  these  concepts  poses  new  challenges  to  the  dynamic  ana¬ 
lyst.  A  number  of  assumptions  in  traditional  analysis  are  held  inviolate,  and  analyst 
must  define  new  approaches  in  order  to  evaluate  designs. 

At  McDonnell  Douglas  Helicopter  Company,  the  high-speed  rotorcraft  challenge 
is  being  addressed  by  the  “TrailRotor  Convertiplane  (TRC)”  Concept.  The  concept 
calls  for  a  high  speed  aircraft  with  hovering  capability.  The  design  considered 
is  a  twin-rotor  configuration  (such  as  in  a  tiltrotor  configuration)  that  meets  the 
hovering  and  low-speed  forward  flight  requirement,  and  jet  engines  to  provide  the 
high-speed  forward  flight  requirements.  In  the  transition  from  helicopter  to  cruise 
conditions,  the  rotor  pylons  tilt  aft  through  ninety  (90)  degrees,  simultaneously 
slowing  RPM  to  near  zero,  and  eventually  stop  and  fold  (about  the  flapping  hinge). 

The  analytical  problem  at  hand  is  to  (i)  examine  the  feasability  of  the  concept, 
(ii)  examine  the  requirements  of  the  rotor  and  fuselage  designs,  (iii)  suggest  designs 
to  promote  the  doability  of  the  approach. 

Existing  analytical  tools  were  used  in  performing  the  analyses  presented  in  this 
paper.  The  key  issue  to  be  examined  is  the  feasability  of  the  conversion  from 
helicopter  mode  to  aircraft  (or  cruise)  mode,  and  the  loads  generated  during  the 
conversion.  Linearized  stability  and  gust  response  during  conversion  are  additional 
items  to  be  examined. 


The  conversion  requirements  are  steady  forward  flight  conditions  at  135  Kts, 
1G  load  condition,  and  nominally  10  seconds  for  complete  conversion.  Attendant 
requirements  are  small  excursions  in  hub  and  ship  pitching,  rolling  moments,  sim¬ 
plicity  of  conversion  requirements  fro 


mm  a  design  and  pilot  workload  standpoint,  and 


“conversion"  includes  the  shaft 


minimal  drag  forces  due  to  conversion.  The  term 
rotation  through  90  degrees,  as  well  as  the  RPM  bleed  from  100  percent  to  near 
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QO  =  HOVER  POWER 
WO  =  NOMINAL  RPM 

Figure  1:  Quui-iteady  RPM  and  Torqus  Schedules  During  Conversion 

zero.  Several  approaches  to  conversion  are  possible.  Two  specific  approaches  are 
considered  in  this  paper.  The  analysis  was  performed  using  the  Comprehensive 
Analytical  Model  for  Rotorcraft  Aerodynamics  and  Dynamics  (CAMRAD/JA). 

A  quasi-steady  equilibrium  approach  is  taken  in  performing  the  analysis.  By 
employing  the  rotor  in  a  windmill  brake  state,  it  is  possible  to  simplify  the  require¬ 
ments  of  the  transmission  system.  In  a  windmill  brake  state,  specific  requirements 
on  the  rotor  torque  can  be  defined  during  both  RPM  bleed  as  well  as  speed-up. 
Given  these  requirements,  it  is  possible  to  find  trim  angles  that  will  control  the 
torque  and  RPM  of  the  rotor.  The  assumptions  made  in  performing  the  analysis 
and  the  results  obtained  are  outlined  in  the  paper. 

Consider  the  RPM  and  the  corresponding  torque  schedules  versus  shaft  angle 
shown  in  figure  1.  The  RPM  schedule  is  assumed;  and  corresponding  to  it,  is  the 
torque  requirements  that  can  be  fairly  easily  derived.  A  trim  analysis  may  now  be 
performed  at  each  intermediate  shaft  angle,  at  the  desired  RPM,  and  the  desired 
torque  requirement.  The  trim  control  angles  so  obtained  determine  the  necessary 
schedule  in  order  to  require  the  appropriate  rotor  torque  and  therefore  result  in  an 
RPM  bleed.  The  effects  of  the  shaft  tilt  dynamics  are  not  included  in  the  model.  As 
the  rotor  tilts  back,  in  near  auto-rotation,  and  a  torque  limit  is  specified  (derived 
from  a  desired  RPM  rate),  there  is  a  non-zero  thrust  generated  by  the  rotor.  In 
response  to  this  thrust,  the  blades  respond  through  flapping.  The  flapping  generated 
must  be  minimal,  in  order  to  keep  the  rotor  controllable  and  minimize  blade  and 
hub  loads.  Since  the  rotor  tilts  aft,  the  thrust  manifests  itself  as  an  apparent  drag 
force.  Additional  requirements  on  the  conversion  process  are  therefore  to  maintain 
small  thrust  as  well  as  small  flapping.  The  net  torque  on  the  rotor  can  be  controlled 
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Figure  2:  Quasi-steady  Control  Requirement*  During  Conversion 

by  prescribing  the  collective,  while  the  flapping  can  be  controlled  by  prescribing  the 
cyclic.  Since  the  cyclic  flapping  is  controlled  to  a  small  amount  (2  degrees),  the 
hub  moments  are  also  small.  For  the  RPM  and  torque  schedules  shown  in  figure 
1,  the  required  control  angles  are  presented  in  figure  2.  The  drag  force  penalty  is 
presented  in  figure  3.  The  only  parameters  that  can  reduce  the  drag  penalty  are 
twist,  RPM  and  allowable  power.  A  parametric  variation  of  the  effects  of  twist 
variation  is  obvious  from  figure  3. 

The  quasi-Bteady  trim  conditions  obtained  at  each  of  the  intermediate  shaft 
angles  may  be  used  to  perform  a  stability  and  linearized  gust  response  analysis. 
These,  and  other  analyses  performed  to  characterize  possible  approaches  to  conver¬ 
sion,  the  limitations  imposed  on  these  analyses  due  to  analytical  assumptions  and 
requirements  of  future  analytical  tools  are  presented  in  the  paper. 
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1  OBJECTIVE 


The  effects  of  refined  aerodynamic  modeling  on  aeroelastic  blade  stability  are  inves¬ 
tigated.  An  advanced  non-linear  aerodynamic  model  based  on  an  indicial  response  method 
is  included  into  a  comprehensive  coupled  rotorcraft  analysis.  The  aerodynamic  analysis  ac¬ 
curately  predicts  blade  loadings  in  attached  and  separated  flow,  as  well  as  dynamic  stall. 
Compressibility  effects  are  implicitly  included  in  the  formulation.  In  addition,  advanced  wake 
models,  both  prescribed  and  free  wake  analyses,  are  coupled  into  the  trimmed,  rotor  response 
and  stability  solution.  The  rotor  analysis  is  formulated  using  a  finite  clement  technique  in 
space  and  time.  Rotor  equations  of  motion  and  trim  equilibrium  equations  are  solved  as  one 
coupled  solution.  Aeroelastic  stability  is  determined  utilizing  Floquet  theory  for  a  linearized 
system.  Blade  stability  for  the  complete  non-linear  system  is  determined  using  a  modified 
Sparse  Time  Domain  technique.  The  rotor  analysis  is  validated  in  a  correlation  study  with 
Aerospatiale  Gazelle  flight  test  data.  The  refined  aerodynamics  significantly  affect  the  blade 
loads,  and  at  extreme  flight  conditions,  small  perturbations  in  the  trim  state,  substantially 
alter  these  loadings.  Consequently,  the  refined  aerodynamics  play  an  important  role  in  de¬ 
termination  of  blade  stability.  The  effects  of  flow  separation,  dynamic  stall,  compressibility 
and  refined  wake  models,  as  well  as  other  key  parameters,  on  rotor  blade  stability  are  deter¬ 
mined.  Results  are  validated  in  a  correlation  study  with  model  hingeless  rotor  experimental 
data. 


2  INTRODUCTION 


Determination  of  the  aeroelastic  stability  of  a  rotor  in  forward  flight  presents  a  clial- 
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lenging  problem  to  researchers.  This  problem  is  magnified  when  considering  high  speed  flight 
conditions  near  the  extremes  of  a  vehicles’  flight  envelope.  Under  such  conditions,  the  aero¬ 
dynamic  environment  of  the  rotor  is  highly  non-linear,  exhibiting  regions  of  dynamic  stall, 
flow  separation  and  transonic  conditions.  At  low  speed,  these  non-linear  effects  are  reduced, 
but  the  complexities  of  the  wake  structure  also  present  challenges  to  the  accurate  prediction 
of  rotor  blade  stability. 

Due  to  the  complexities  of  this  problem,  many  simplifications  are  made  to  make  the 
analysis  more  manageable.  In  recent  work  [1,2],  simplifications  include  rigid  blades,  quasi¬ 
steady  aerodynamics,  and  uniform  or  Drees  inflow  models  [3].  Unsteady  effects  are  often 
approximated  with  a  dynamic  inflow  model  [4].  This  method  is  limited  to  low  frequency 
effects,  and  thus  cannot  predict  detailed  flow  phenomena.  Additional  investigations  are 
cited  in  Reference  [5].  These  approximations  yield  satisfactory  results  under  moderate  flight 
conditions.  At  more  extreme  flight  conditions,  however,  non-linear  aerodynamic  effects  sig¬ 
nificantly  affect  blade  loadings,  and  thus  simpler  models  are  no  longer  applicable.  At  low 
speeds,  the  linear  Drees  inflow  model  breaks  down,  and  more  advanced  wake  models  are 
required.  The  present  study  utilizes  a  refined  non-linear  aerodynamic  model,  as  well  as  pre¬ 
scribed  and  free  wake  models,  to  predict  rotor  blade  stability  for  several  flight  conditions. 
The  analysis  is  validated  in  a  correlation  with  model  hingeless  rotor  test  data. 


3  ROTOR  ANALYSIS 

3.1  AERODYNAMIC  MODEL 


In  an  effort  to  capture  the  complete  range  of  flow  characteristics,  and  still  maintain 
simplicity  and  efficiency,  a  recent  aerodynamic  model,  developed  by  Beddoes  and  Leishman 
[6,7]  is  adopted.  This  model  accurately  predicts  attached  and  separated  flow,  as  well  as 
dynamic  stall.  An  important  feature  of  the  model  is  that  compressibility  effects  are  implicitly 
included.  The  model  requires  only  a  few  empirical  factors,  most  of  which  are  derived  from 
static  airfoil  data.  The  parameters  derived  from  the  unsteady  data,  however,  are  insensitive 
to  airfoil  shape,  and  thus  need  not  be  recalculated  for  different  a'.rfoils.  The  model  has 
been  extensively  validated  for  a  series  of  different  airfoil  sections  over  a  wide  range  of  flight 
conditions.  This  model  has  been  recently  incorporated  into  a  comprehensive  rotor  aeroelastic 
analysis  [8].  This  work  carried  out  a  parametric  study  to  determine  the  significance  of  the 
refined  aerodynamics  on  the  dynamic  response  of  rotors  in  forward  flight.  Results  showed 
that  at  high-speed  conditions,  the  non-linear  aerodynamics,  or  separated  and  stalled  flow, 
dominate  the  outboard  blade  section  forces,  and  significantly  affect  blade  response  and  blade 
root  loads.  These  effects  were  amplified  at  increased  thrust  levels.  Compressibility  effects 
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were  also  found  to  substantially  affect  the  extent  of  separated  flow  on  the  rotor  blade. 


3.2  WAKE  MODELS 


Refined  wake  modeling  has  proven  to  be  very  important  under  certain  flight  condi¬ 
tions.  Recent  literature  has  emphasized  the  need  for  further  refinements  to  improve  cor¬ 
relations  of  experimental  inflow  measurements  with  predicted  values.  Current  technology 
provides  several  models,  from  classical  skewed  helices  to  free  wake  calculations,  and  compar¬ 
isons  have  been  made  on  their  respective  ability  to  predict  rotor  inflow  over  a  range  of  flight 
conditions  [9,10,11].  The  present  analysis  incorporates  both  the  prescribed  and  free  wake 
models  from  CAMRAD  [12], 


3.3  COUPLED  TRIM  AND  RESPONSE 


The  trim  solution  determines  rotor  pitch  controls  and  vehicle  orientation  for  a  propul¬ 
sive  flight  condition.  The  steady  response  determines  the  deflected  position  of  the  blade  for 
one  complete  azimuthal  cycle.  The  present  investigation  utilizes  an  analysis  which  solves  the 
coupled,  non-linear  periodic  equations  of  blade  motion  and  six  trim  equilibrium  equations. 
The  analysis  is  based  on  a  finite  element  approximation  in  both  space  and  time.  Each  blade 
is  assumed  to  be  an  elastic  beam  undergoing  flap  bending,  lag  bending,  elastic  twist  and 
axial  deflection.  The  blades  are  discretized  into  a  number  of  beam  elements,  each  having 
fifteen  degrees  of  freedom.  To  reduce  computational  time,  a  large  number  of  finite  element 
equations  are  transformed  to  the  modal  space  as  a  few  normal  mode  equations  (6-9  modes 
are  normally  used).  The  periodic,  non-linear,  response  solution  is  determined  by  using  a 
temporal  finite  element  discretization.  Fourth  order  Lagrangian  shape  functions  are  utilized 
within  each  trie  element.  The  finite  element  in  time  formulation  is  based  on  Hamilton’s 
weak  principle.  Resulting  normal  mode  equations  are  ultimately  solved  as  a  set  of  non-linear 
algebraic  equations.  The  solution  is  obtained  utilizing  an  iterative,  modified  Newton  method 
[1,13]. 


Concurrent  with  the  solution  of  the  blade  equations  of  motion,  is  the  solution  of 
the  vehicle  trim  equations.  The  vehicle  is  trimmed  on  all  axes  for  both  force  and  moment 
equilibrium  (6  equations).  The  modified  Newton  solution  method  allows  the  two  solutions 
to  be  coupled  to  produce  a  unique  solution  that  satisfies  both  blade  and  vehicle  equilibrium 


conations  The  force  summation  method  is  used  to  determine  hub  foices  as  a  sum  of  inertial 


and  aerodynamic  loadings  [14]. 
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3.4  AEROELASTIC  STABILITY 


Linearized  blade  stability  is  determined  using  Floquet  transition  matrix  theory  [15]. 
A  normal  mode  approach  is  used,  utilizing  the  coupled  blade  modes  determined  about  the 
trimmed  mean  deflections.  The  perturbation  equations  of  motion  are  numerically  integrated 
to  determine  the  transition  matrix.  A  complex  eigenanalysis  is  performed  on  this  matrix 
to  assess  the  stability  of  each  blade  mode.  Modal  decay  rates  are  derived  from  resulting 
eigenvalues.  A  positive  decay  rate  denotes  an  unstable  mode. 

Rotor  blade  stability  for  a  non-linear  response  is  determined  using  a  modified  Sparse 
Time  Domain  (STD)  technique,  recently  developed  [16].  The  STD  technique  is  a  multi¬ 
degree  of  freedom  technique  for  estimating  the  modal  damping  parameters  from  transient 
response  data.  Solutions  of  system  equations  are  performed  utilizing  a  windowing  technique. 
This  method  ultimately  determines  the  damping  characteristics  for  each  mode,  which  in  turn 
determine  blade  stability. 


4  RESULTS  AND  CONCLUSIONS 


The  effects  of  refined  aerodynamic  analyses  have  been  studied  in  depth  in  a  correlation 
study  with  flight  test  data  [17].  The  data  used  is  from  a  modified,  experimental  SA349/2 
Modified  Gazelle  helicopter  [18].  The  SA349/2  is  outfitted  with  an  Aerospatiale  advanced 
geometry  rotor.  This  articulated  rotor  includes  a  damped  lag  articulation  hub,  and  three 
high  speed  (Grande  Vitesse)  research  blades.  The  rotor  blade  airfoil  section  is  an  OA209, 
9  percent  thickness  airfoil,  and  the  blade  has  a  non-linear  twist  distribution. 

The  correlation  included  validation  of  the  analysis  for  three  flight  conditions,  varying 
both  advance  ratio  and  thrust  coefficient;  low  speed-low  thrust,  high  speed-low  thrust,  and 
high  speed-high  thrust  [17].  Aerodynamic  blade  loads  (normal  force  and  pitching  moment) 
and  blade  bending  moments  (flap  and  lag)  were  investigated.  For  a  moderate  flight  condition, 
high  speed-low  thrust,  blade  section  normal  force  was  well  predicted,  Figure  1.  Non-linear 
aerodynamic  effects  were  small  for  this  flight  condition.  Flap  bending  moments  were  also 
satisfactorily  predicted,  Figure  2.  Use  of  a  refined  wake  analysis  proved  necessary  to  predict 
the  harmonics  of  the  bending  moments.  At  high  speed,  prescribed  and  free  wake  models 
gave  similar  results.  At  the  low  speed  flight  condition,  the  free  wake  model  substantially 
improved  the  correlation  of  blade  section  normal  force  on  the  retreating  side  of  the  disk, 
Figure  3.  At  a  high  speed-high  thrust  flight  condition,  the  non-linear  aerodynamic  effects 
were  quite  significant.  Pitching  moment  correlations  were  improved  by  including  non-linear 
aerodynamic  effects,  Figure  4.  The  effect  of  small  perturbations  in  the  rotor  thrust  were 
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also  studied.  A  small  increase  in  thrust  significantly  altered  the  blade  section  pitching 
moment  as  dynamic  stall  occurred  on  the  blade  giving  rise  to  a  large  spike  in  the  resulting 
moment,  Figure  5.  The  sensitivity  of  blade  loads  to  such  perturbations  are  quite  important  in 
determining  blade  stability.  Additional  results  depicted  good  correlation  of  trimmed  vehicle 
controls  and  vehicle  orientation,  and  a  substantial  influence  of  blade  twist  and  compressibility 
effects  on  the  extent  of  separated  flow  on  the  rotor. 

The  validation  of  the  analysis  yields  confidence  in  the  predictive  capabilities  of  the 
refined  aerodynamics.  The  current  research  effort  extends  these  advanced  analyses  into  the 
determination  of  rotor  blade  stability.  Results  are  being  correlated  with  model  hingeless 
rotor  test  data  to  validate  the  analyses  and  assess  the  effects  of  the  refined  aerodynamics. 
The  test  configuration  was  an  isolated,  soft  in-plane,  rotor  with  discrete  flap  and  lead-lag 
flexures  and  relatively  rigid  blades.  Data  was  acquired  for  both  hover  and  forward  flight  [19]. 
Preliminary  results  show  that  a  linearized  Floquet  analysis  yields  accurate  results  under 
moderate  flight  conditions,  Figures  6-7.  At  more  extreme  conditions,  namely  increasing 
advance  ratios,  a  discrepancy  arises  in  the  correlation.  Further  investigation,  utilizing  a 
non-linear  modified  STD  technique,  will  determine  if  these  discrepancies  are  attributable 
to  non-linear  aerodynamic  effects.  Ultimately,  the  effects  of  flow  separation,  dynamic  stall, 
compressibility,  refined  wake  modeling  and  other  key  parameters  on  rotor  blade  stability  will 
be  investigated  over  a  range  of  flight  conditions. 
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Abstract 


Rotary  wing  aeroelasticity  has  been  an  area  of  considerable  activity  for  the  last  one  dec¬ 
ade.  A  lot  of  work  has  gone  into  developing  an  accurate  and  efficient  model  of  rotor 
blades  for  analysis  purposes,  most  of  which  have  been  referred  to  in  the  review  by 
Friedmann  [1],  One  of  the  prime  areas  of  research  has  been  structural  modelling,  with  a 
continued  interest  shown  in  geometrically  non-linear  models  for  hingeless  and 
bearingless  rotor  blades.  Inspite  of  this  spate  of  research  .there  is  a  distinct  lack  of  liter¬ 
ature  on  the  analysis  of  unsymmetrically  laminated  composite  rotor  blades. 

A  simple  and  efficient  element  has  been  developed  by  Kapania  and  coworkers  [Refs.  2 
and  3]  for  non-linear  analysis  of  unsymmetrically  laminated  beams  and  for  aeroelastic 
response  of  wing  boxes  with  laminated  skins.  The  ongoing  research  Intends  to  include 
study  of  the  aeroelastic  response  of  rotor  blades.  The  element  can  be  used  to  analyse 
built-up  box  structures  reinforced  by  stringers.  The  formulation  takes  into  account  the 
effects  of  transverse  shear,  the  bending-twist  and  the  bending-stretching  coupling  inherent 
in  the  blades  with  unsymmetric  skins  and  unsymmetric  cross-sections.  An  important 
feature  of  this  element  is  that  even  though  it  is  a  one-dimenslonal  element,  it  is  capable 
of  simulating  the  structural  behaviour  of  box  type  structures  wiih  the  same  accuracy  as 
that  obtained  using  WIDOWAC. 

The  analysis  accounts  for,(i)  flap-lag  coupling  due  to  the  unsymmetric  nature  of  the 
cross-section  of  the  rotor  blades;  (ii)  axial-  bending  coupling,  and  the  axial-torsion  coupl¬ 
ing  due  to  the  either  the  anisotropic  nature  of  the  skin  or  the  presence  of  sweep  as  in  the 
case  of  rotor  blades  with  swept  tips. 

The  effect  of  axial  forces  is  considered  in  the  Incremental  stiffness  matrix  [4]  and  the 
non-linear  load-displacement  path  is  predicted  by  a  linearized  mid-point  tangent  incre¬ 
mental  procedure  [Refs.  5  and  6].  The  generalised  aerodynamic  forces  is  obtained  using 
the  principle  of  virtual  work  in  conjunction  with  the  2-dimensional  quasi-steady  aerodyna¬ 
mic  theory.  The  results  obtained  will  be  compared  with  other  works  to  validate  the 
method. 
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EXTENDED  ABSTRACT 
Introduction  and  Problem  Statement 

The  dynamic  and  aeroelastic  problems  of  hingeless  rotor  helicopters  have 
been  the  subject  of  a  considerable  amount  of  research  in  the  past  twenty 
years.  Comprehensive  reviews  of  such  work  have  been  published  by  Fried¬ 
mann  [1],  Johnson  [2],  Ormiston  [3],  and  Ormiston  et  al.  [4].  One  of  the 
advantages  of  hingeless  rotor  helicopters  over  articulated  rotor  helicopters 
is  greater  agility.  Therefore  it  is  rather  surprising  that  little  or  no  infor¬ 
mation  has  been  published  on  the  aeroelastic  stability  of  hingeless  rotors  in 
maneuvering  flight. 

The  analytical  and  experime.  tal  studies  on  hingeless  rotor  aeroelastic 
stability  have  generally  focused  on  steady,  straight,  1-g  level  flight  conditions. 
Maneuvering  flight  conditions  affect  the  rotor  dynamics  mainly  by  changing 
the  load  factor,  by  changing  the  orientation  of  the  resultant  rotor  forces  and 
moments,  and  by  introducing  pitch  and  roll  angular  velocities.  Secondary 
effects  are  due  to  the  changing  of  pitch  control  settings,  and  the  generation 
of  sideslip  angles.  The  effect  of  these  factors  on  the  aeroelastic  stability  of 
the  rotor  has  been  largely  ignored  in  the  published  literature. 

The  objectives  of  the  proposed  paper  are  the  following: 

1.  to  present  the  results  of  a  numerical  study  of  the  aeroelastic  stability 
and  response  of  a  trimmed  hingeless  rotor  in  a  steady  coordinated  turn, 

2.  to  study  the  effects  of  steady  stall  and  compressibility  airfoil  charac¬ 
teristics  on  the  trim  conditions  of  the  aircraft,  and  on  the  aeroelastic 
stability  of  the  rotor. 
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Mathematical  model 


The  rotor  blades  are  modeled  as  isotropic,  Bernoulli-Euler  beams  under¬ 
going  moderate  deflections  [5].  The  nonlinear,  partial  differential  equations 
of  motion  of  the  blade  are  discretized  using  a  finite  element  formulation  based 
on  Galerkin  method  of  weighted  residuals  [6]. 

The  formulation  of  the  rotor  aeroelastic  problem  is  based  on  the  ” im¬ 
plicit”  approach  introduced  in  Ref.[7)  for  the  aerodynamic  operator  of  the 
equations  of  motion,  and  extended  in  Ref.[8]  to  the  inertia  operator.  Using 
this  approach,  the  mathematical  expressions  for  the  aerodynamic  and  inertia 
loads  are  not  expanded  symbolically.  Instead,  they  are  built  numerically  as 
part  of  the  solution  process.  This  approach  all  but  eliminates  the  tedious  and 
error  prone  task  of  deriving  long,  explicitly  expanded  equations  describing 
the  complex  nonlinear  dynamics  of  the  elastic  blades  attached  to  a  rotating 
and  translating  hub. 

The  equations  describing  the  trim  state  of  the  aircraft  in  the  coordinated 
turn  are  based  on  those  presented  by  Chen  [9]: 

1.  Force  equilibrium  along  the  x  body  axis: 

X  —  mg  sin  8  +  m(qw  —  rv) 

2.  Force  equilibrium  along  the  y  body  axis: 

Y  =  —mg  cos  0  sin  <f>  +  m(ru  —  pw) 

3.  Force  equilibrium  along  the  z  body  axis: 

Z  =  —mg  cos  6  cos  <j>  -f  m(pv  —  qu ) 

4.  Roll  moment  equilibrium: 

L  =  ~Iyz{q2  ~r2)~  Ixzpq  +  Ixypr  -  { Iy  -  Iz)qr 

5.  Pitch  moment  equilibrium: 


M  —  !Xz{t  p )  ^xpqT  4"  lyzpq  (I*  ^x)yp 
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The  airspeed  components  u,  v,  and  w  in  the  aircraft  body  axes  system  x, 
y ,  z  are  given  by: 


V  cos  a  cos p 

(1) 

V  sin  ft 

(2) 

V  sin  a  cos  p 

(3) 

The  13  unknowns  are  the  collective  pitch  of  the  main  rotor  90  and  of  the 
tail  rotor  6t,  the  longitudinal  and  lateral  cyclic  pitch  settings  9U  and  9ic,  the 
steady  state  roll,  pitch,  and  yaw  rates  p ,  q ,  and  r,  the  Euler  pitch  angle  9, 
the  aerodynamic  angle  of  attack  o,  the  sideslip  angle  /?,  and  the  constant 
portion  of  the  inflow  for  the  main  rotor  A  and  the  tail  rotor  Xt. 

The  turn  rate  ^  (>  0  for  a  right  turn),  the  magnitude  V  of  the  aircraft 
velocity  vector,  and  the  climb  angle  7  (>  0  for  an  ascending  turn)  define  the 
turning  flight  condition,  and  are  provided  as  input. 

Detailed  expressions  for  the  force  components  X,  Y,  and  Z,  and  for  the 
moment  components  L,  M ,  and  AT,  are  not  presented  in  this  abstract  for 
brevity,  but  will  be  provided  in  the  complete  version  of  the  paper. 

Steady  stall  and  compressibility  effects  are  taken  into  account  through  the 
aerodynamic  coefficients  of  the  airfoil.  Semi-empirical  analytical  expressions 
for  the  normal  force  coefficient  c„,  the  drag  coefficient  q,  and  the  moment- 
coefficient  cm,  as  a  function  of  angle  of  attack  a  and  Mach  number  M  have 
been  presented  by  Beddoes  for  the  NACA  0012  airfoil  [10].  The  values  of  cn, 
cm,  and  Cd  generated  using  these  expressions  are  plotted  in  Figures  1,  2,  and 
3  respectively,  for  a  range  of  values  of  a  and  M.  These  expressions  are  used 
in  the  solution  of  both  the  coupled  trim  problem  and  the  aeroelastic  stability 
and  response  problem. 

Three-dimensional  tip  effects  are  taken  into  account  in  an  approximate 
way,  by  using  Goldstein’s  circulation  correction  factor  k  [11,  pp. 103-1 10], 
and  multiplying  the  normal  force  coefficient  cn  and  the  moment  coefficient 
cm  at  a  given  blade  section  by  k. 

Solution  Technique 

The  solution  process  is  divided  into  two  phases,  namely: 

1.  determination  of  the  trim  state  of  the  helicopter  in  the  turn,  and 
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2.  calculation  of  the  aeroelastic  response  and  stability  of  (he  rotor. 

Phase  I — Trim  in  coordinated  turn 

The  determination  of  the  trim  state  of  the  helicopter  is  coupled  with  the 
calculation  of  the  steady  state  equilibrium  position  of  the  rotor.  For  the  trim 
calculations  the  system  of  nonlinear  ordinary  differential  equations  (ODE)  of 
motion  of  the  blade  is  transformed  into  a  set  of  nonlinear  algebraic  equations 
using  a  classical,  or  global  Galerkin  method,  as  described  in  Ref.  [12].  Before 
Galerkin  method  is  applied,  a  modal  coordinate  transformation  is  performed. 
If  m  modes  are  used  in  the  coordinate  transformation,  and  each  generalized 
coordinate  is  approximated  with  a  Fourier  series  truncated  at  the  n-th  har¬ 
monic,  then  the  application  of  Galerkin  method  results  in  the  generation 
of  m(2n  +  1)  nonlinear  algebraic  equations.  The  set  of  13  trim  equations, 
Eq.  (1)  through  (13)  is  appended  to  those  equations,  and  the  combined  sys¬ 
tem  is  solved  simultaneously. 

Phase  II — Rotor  aeroelastic  stability  and  response 

For  the  solution  of  the  rotor  aeroelastic  problem,  the  rotor  is  assumed  to 
be  attached  to  a  hub  of  infinite  mass,  that  rolls,  pitches,  and  yaws  at  the 
rates  calculated  in  the  trim  process.  The  problem  of  solving  the  system  of 
nonlinear  ODE  of  motion  of  the  blades  is  transformed  into  that  of  solving  a 
sequence  of  linearized  ODE,  using  quasilinearization.  Because  the  equations 
have  periodic  coefficients,  the  stability  of  the  linearized  system  of  ODE  is 
calculated  using  Floquet  theory.  The  details  of  the  solution  procedure,  for 
the  case  of  implicitly  formulated  aerodynamic  and  inertia  loads,  are  presented 
in  Refs.  [7]  and  [8]. 


Results 

The  results  presented  in  this  section  refer  to  a  soft-in-plane  blade  config¬ 
uration,  with  fundamental,  rotating,  coupled  natural  frequencies  of  0.73/rev, 
1.12/rev,  and  3.17/rev  in  lag,  flap,  and  torsion  respectively.  The  blades  have 
uniform  mass  and  stiffness  properties  along  the  span.  In  steady,  straight  and 
level  flight  the  rotor  is  operating  at  a  Gr  —  0.005.  The  solidity  is  o  =  0.07. 
Blade  precone,  built-in  twist,  and  cross-sectional  offsets  from  the  elastic  axis 
are  equal  to  zero. 

The  first  set  of  results  was  obtained  with  a  linear,  incompressible  airfoil 
model,  and  neglecting  three  dimensional  aerodynamic  effects.  Results  were 
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derived  for  values  of  the  advance  ratio  /x  =  0.1,  0.2,  and  0.3,  and  for  turn  rates 
=  0,  0.1, 0.2, 0.4,  and  0.6  rad/sec,  corresponding  to  straight  flight,  and  to 
one  complete  360  degree  turn  in  63,  31,  16,  and  11  seconds  respectively.  The 
load  factors  corresponding  to  the  various  combinations  of  speed  and  turn  rate 
are  summarized  in  Figure  4.  Three  coupled  modes  were  used  to  calculate  the 
trim  condition.  The  modes  were  the  first  flap,  the  first  lag,  and  the  first 
torsion  coupled  modes.  Each  generalized  modal  coordinate  was  represented 
as  a  two  harmonics  Fourier  series.  Therefore  m  =  3  and  n  =  2,  for  a  total 
of  15  equations  associated  with  the  rotor  elastic  deformations.  Including 
the  13  equations  associated  with  the  fuselage  equations,  the  trim  state  of  the 
helicopter  is  defined  by  a  system  of  28  algebraic  equations.  The  quasi-Newton 
solver  HYBRD  [13]  was  used  to  solve  the  system.  The  solution  required 
between  110  and  130  seconds  of  CPU  time  on  an  IBM  3081D  mainframe. 

The  aeroelastic  stability  calculations  were  conducted  using  six  coupled 
modes.  Each  generalized  coordinate  was  represented  as  a  ten  harmonics 
Fourier  series.  The  calculations  were  conducted  at  a  fixed  advance  ratio, 
for  increasing  turn  rate.  The  initial  approximation  to  the  response  at  a 
given  turn  rate  was  the  solution  at  the  previous  turn  rate;  the  initial  ap¬ 
proximation  for  the  straight  flight  cases  was  a  rigid  blade.  Two  iterations  of 
quasilinearization  were  sufficient  for  all  but  one  of  the  turning  flight  condi¬ 
tions;  three  iterations  were  required  for  the  straight  flight  conditions.  Each 
iteration  required  180-200  seconds  of  CPU  time. 

Figures  5  and  6  show  the  real  part  of  the  characteristic  exponent  of  the 
first  and  second  lag  mode  respectively,  as  a  function  of  the  advance  ratio 
and  turn  rate.  The  real  parts  of  the  characteristic  exponents  are  indicative 
of  the  damping,  for  a  system  with  periodic  coefficients.  Figures  5  and  6 
show  that  a  right  turn  has  a  stabilizing  effect  on  these  two  lowly  damped 
modes,  and  that  this  effect  becomes  larger  with  increasing  advance  ratio. 
The  first  flap  mode  and  the  first  torsion  mode  are  essentially  unaffected  by 
the  turn  rate,  therefore  the  appropriate  results  are  not  presented  here.  The 
second  and  third  flap  modes  are  somewhat  destabilized  by  a  right  turn.  The 
characteristic  exponent  of  the  second  flap  mode  is  shown  in  Figure  7.  The 
damping  of  the  second  flap  mode  is  so  high  that  the  mode  remains  well 
damped  even  at  high  load  factors. 

Additional  results  were  obtained  for  the  same  blade  configuration,  us¬ 
ing  the  nonlinear  airfoil  characteristics  mentioned  previously.  Figures  8 


through  10  show  the  real  parts  of  the  characteristic  exponents  of  the  first 
lag,  second  lag,  and  first  torsion  mode  respectively.  The  advance  ratio  is 
[i  =  0.2,  the  turn  rate  is  ip  =  0.25  rad/sec,  and  the  thrust  coefficient  is 
Ct/<t  w  0.13. 

Stall  and  compressibility  affect  directly  the  stability  of  the  blade  by  mod¬ 
ifying  the  aerodynamic  load  distribution.  Because  rotary-wing  aeroelastic 
stability  problems  are  inherently  nonlinear,  stall  and  compressibility  also 
affect  the  stability  level  of  the  blade  indiitctly  by  modifying  the  trim  condi¬ 
tions  of  the  helicopter,  including  the  pitch  control  settings.  In  an  attempt  to 
separate  the  two  mechanisms,  Figures  8  through  10  each  show  three  plots. 
The  first  refers  to  results  obtained  by  neglecting  stall  and  compressibility  in 
both  the  trim  and  the  aeroelastic  stability  calculations;  in  the  second,  stall 
and  compressibility  were  modeled  in  the  aeroelastic  stability  and  response 
calculations  but  not  in  the  trim  calculations;  in  the  third,  stall  and  com¬ 
pressibility  are  modeled  in  both  the  trim  and  the  aeroelastic  calculations. 
Each  of  the  figures  shows  the  characteristic  exponents  computed  at  each  it¬ 
eration  of  quasilinearization.  The  number  of  iterations  required  to  achieve 
convergence,  and  the  size  of  the  changes  of  the  characteristic  exponents  as 
the  iterations  proceed,  give  an  indication  on  the  strength  of  the  nonlinearities 
in  the  mathematical  model. 

Figures  8  through  10  indicate  that  stall  and  compressibility  have  a  stabi¬ 
lizing  effect  on  the  damping  of  the  first  and  second  lag  modes,  and  a  desta¬ 
bilizing  effect  on  the  first  torsion  mode.  Stall  and  compressibility  increase 
the  nonlinearity  of  the  problem.  This  is  indicated  by  the  need  for  an  addi¬ 
tional  iteration  of  quasilinearization  when  these  effects  are  present.  From  a 
numerical  standpoint,  quasilinearization  converged  without  problems  when 
stall  and  compressibility  were  included,  but  the  computer  time  required  to 
complete  one  iteration  increased  to  450-500  seconds  of  CPU  time.  This  is 
probably  due  to  the  smaller  time  steps  that  the  variable-step,  variable-order 
ODE  solver  [7]  selected  to  deal  with  faster  varying  aerodynamic  loads. 

Additional  results  are  currently  being  derived,  to  gain  a  better  under¬ 
standing  of  the  physical  mechanisms  that  contribute  to  modify  the  aeroelas¬ 
tic  stability  of  the  rotor.  In  particular,  results  are  being  derived  to  assess 
separately  the  effects  of  stall  and  of  compressibility.  These  results  will  be 
presented  in  the  final  version  of  the  paper. 


7 


Conclusions 


The  study  outlined  in  this  abstract  represents  a  first  step  toward  a  basic 
understanding  of  an  important,  yet  largely  ignored  practical  problem,  namely 
the  aeroelastic  behavior  of  a  hingeless  rotor  in  maneuvering  flight. 

The  preliminary  results  presented  in  this  abstract  indicate  the  following 
trends: 

1.  When  stall  and  compressibility  effects  are  neglected,  the  lag  modes 
appear  to  be  strongly  stabilized  by  a  turn.  The  stabilizing  effect  in¬ 
creases  with  increasing  turn  rate  and  advance  ratio.  Some  of  the  higher 
modes  may  be  destabilized  by  a  turn  at  very  high  load  factor,  but  their 
damping  level  in  rectilinear  flight  is  high,  and  they  remain  stable  in  a 
turn. 

2.  Preliminary  results  on  the  efTect  of  stall  and  compressibility  suggest 
that  the  damping  of  the  lag  modes  increases  when  these  effects  are 
taken  into  account.  All  the  remaining  modes  are  destabilized,  but  not 
to  the  extent  that  they  become  unstable.  Stall  and  compressibility 
act  directly  by  modifying  the  aerodynamic  loads,  and  indirectly  by 
modifying  the  helicopter  trim  condition. 

3.  Stall  and  compressibility  increase  the  nonlinearity  of  the  aeroelastic 
problem.  Quasilinearization  still  converges  reliably,  but  the  computer 
time  required  to  complete  an  iteration  increases  considerably. 
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The  bearingless  rotorcraft  offers  reduced  weight,  structural 
simplicity,  and  superior  flying  qualities.  Bearingless  rotor 
technology  has  been  successfully  applied  to  the  tail  rotors  of  the 
Blackhawk  and  S-76  helicopters  by  Sikorsky.  Boeing  Vertol  built  the 
first  successful  bearingless  main  rotor  and  flew  it  in  1978.  Most  of 
the  next  generation  rotorcraft  are  likely  to  be  designed  with 
bearingless  hubs. 

All  practical  designs  of  bearingless  rotors  have  multiple -load-path 
hubs.  The  bearingless  rotor  tested  by  Boeing  Vertol  has  three  load 
paths:  two  fiberglass  flexbeams  and  a  filament  wound  torque  tube. 

Because  of  the  multiple  load  paths,  the  dynamic  behavior  of 
bearingless  rotors  is  significantly  different  than  conventional 
rotors.  It  is  important,  therefore,  to  have  the  capability  to 
accurately  analyze  bearingless  rotor3  and  how  they  affect  the 
stability  of  the  entire  helicopter. 

CAMRAD  (  Comprehensive  Analytical  Model  for  Rotorcraft  Aerodynamics 
and  Dynamics)  is  a  comprehensive  rotorcraft  analysis  program  developed 
at  NASA-Ames.  The  program  utilizes  new  technology  to  analyze  a  wide 
range  of  rotorcraft  problems.  This  includes  the  trim  solution, 
performance,  loads  and  noise,  stability  derivatives  and  handling 
qualities,  aeroelastic  stability,  and  vibration  and  gust  response. 

CAMRAD  is  limited,  however,  to  single  load  path  modelling.  The 
axial  degree  of  freedom  for  blade  deformation,  which  is  important  for 
multiple  load  paths,  is  also  not  included. 
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The  CAMRAD  approach  to  solving  the  elastic  blade  equations  i3 
termed  the  modified  Galerkin's  method.  Exact,  nonrotating  beam 
solutions  are  used  in  a  Galerkin's  procedure  to  solve  the  rotating  beam 
problem  approximately. 

Procedures  similar  to  the  ones  that  exist  in  CAMRAD  are  being 
developed  for  multiple-load-path  blades.  Once  the  procedures  are 
developed  and  validated,  it  will  be  possible  to  extend  the  CAMRAD 
procedures  to  include  multiple-load-path  bearingless  blades. 

Application  of  the  modified  Galerkin'3  method  to  multiple-load-path 
blades  has  been  completed  and  the  results  are  shown  in  the  next  five 
pages  for  a  model  rotor  tested  by  Dawson.  Work  is  in  progress  for 
obtaining  the  trim  solution. 
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ABSTRACT 


Many  papers  have  been  published  which  analyze  the  behavior  of  cracked 
rotating  shafts  [1].  The  present  investigation  is  more  general  than  previous  work 
in  several  respects.  The  most  important  are  the  inclusion  of  (a)  shear  deformation 
and  rotatory  inertia  (Timoshenko  model),  (b)  longitudinal,  bending,  and  torsional 
motions,  and  (c)  three  types  of  damping. 

Except  for  a  single  transverse  crack,  the  shaft  is  assumed  to  be  uniform 
and  to  have  simply  supported  ends.  Its  axis  need  not  be  horizontal,  and  the 
centroid  of  its  cross-sectional  area  need  not  coincide  with  its  center  of  mass.  The 
material  is  assumed  to  be  viscoelastic,  which  models  internal  damping.  External 
viscous  damping  is  also  included.  In  addition,  when  the  crack  is  closed,  dry 
friction  acts  at  the  crack  face. 

The  crack  is  assumed  to  be  open  when  the  axial  strain  at  the  crack  location 
is  positive,  and  closed  when  the  axial  strain  there  is  negative.  Following  an  idea 
developed  by  Petroski  [2],  the  discontinuities  in  stiffness  and  damping  at  the 
crack  are  represented  by  concentrated  forces.  The  magnitudes  of  these  forces 
are  determined  from  the  compliance  matrix  associated  with  the  crack,  which  was 
obtained  by  the  use  of  fracture  mechanics  in  Reference  [3]. 

The  governing  equations  of  motion  consist  of  six  coupled,  partial 
differential  equations  involving  three  deflections  and  three  rotations,  which  are 
functions  of  time  and  of  position  along  the  shaft  [4].  The  equations  are  linear 
during  each  period  of  time  in  which  the  axial  strain  has  the  same  sign;  when  it 
changes  sign  (i.e.,  when  the  crack  opens  or  closes),  some  of  the  stiffness  and 
damping  coefficients  change.  Hence  the  problem  is  nonlinear. 
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Galerkin’s  method  is  applied  with  two  terms  for  each  of  the  deflection  and 
rotation  functions.  This  leads  to  a  system  of  12  ordinary  differential  equations  for 
the  time-varying  amplitudes  of  the  assumed  global  shape  functions.  The  system 
is  solved  with  the  use  of  numerical  integration.  Time  histories  of  the  deflections 
and  rotations  are  computed  and  their  frequency  spectra  are  obtained. 

Attention  is  focused  on  the  steady-state  response  of  the  shaft.  The  effects 
of  the  crack  depth  and  the  angular  velocity  are  investigated.  Exchanges  of 
energy  between  the  longitudinal,  bending,  and  torsional  motions  are  studied. 
Both  free  and  forced  vibrations  are  examined. 

In  conclusion,  then,  this  paper  generalizes  previous  work  on  cracked 
rotating  shafts  in  several  respects.  The  results  demonstrate  the  effect  of  cracks 
on  shaft  response.  They  improve  our  understanding  of  the  behavior  of  cracked 
shafts,  and  can  be  utilized  in  the  detection  of  cracks  before  the  occurrence  of 
failure. 
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ABSTRACT 

The  dynamic  behavior  of  a  long  slender  beam  is  of  Interest  in  connection  with 
helicopter  rotor  blades,  lengthy  robot  manipulators,  spacecraft  antennae,  flexible 
satellites,  wind  turbine  blades,  and  other  systems  that  perform  large  and/or  complex 
motions. 

In  this  paper,  Newton's  second  law  is  used  to  develop  the  nonlinear  equations 
of  motion  and  boundary  conditions  describing  the  flexural-flexural-torsional  vibration 
of  composite  beams.  The  equations  contain  structural  coupling  terms  and  cubic 
nonlinearities  due  to  curvature  and  inertia.  Three  consecutive  Euler  angles  are  used 
to  evaluate  the  transformation  matrix  between  the  deformed  and  undeformed  states. 
Because  the  twisting-related  Euler  angle  </>  is  not  an  independent  Lagrangian  coordi¬ 
nate,  the  torsional  equation  is  used  to  define  the  real  twist  angle. 

The  nonplanar  responses  of  a  cantilevered,  symmetrically  laminated 
graphite-epoxy  composite  beam  subject  to  lateral  harmonic  base  excitations  (see  Fig. 
1)  are  investigated  using  the  derived  nonlinear  coupled  integro-differential  equations 
of  motion.  A  powerful  technique  based  on  the  state  space  concept  and  the  funda¬ 
mental  matrix  method  is  used  to  solve  for  the  linear  mode  shapes.  Then  a  combi¬ 
nation  of  the  Galerkin  procedure  and  the  method  of  multiple  scales  is  used  to 
construct  a  first-order  uniform  expansion  for  the  case  of  a  primary  resonance  and  a 
one-to-one  internal  resonance  in  which  the  linear  natural  frequency  of  flexural  vi¬ 
brations  in  one  plane  is  approximately  equal  to  the  linear  natural  frequency  of  flexural 
vibrations  in  the  other  piane.  We  determine  the  ranges  of  the  driving  frequency 
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which  result  In  unstable  planar  motions  and  describe  quantitatively  the  ensuing 
nonplanar  motions.  We  study  some  of  the  main  characteristics  of  the  nonplanar 
motions  and  investigate  the  possible  bifurcations  of  the  response  amplitudes  when 
the  excitation  frequency  is  varied.  We  also  demonstrate  the  existence  of  Hopf 
bifurcations  and  the  nonexistence  of  periodic  motions  for  certain  excitation  parame¬ 
ters.  For  different  range  of  parameters,  the  nonplanar  motions  can  be  steady  whirling 
motions,  whirling  motions  of  the  beating  type,  or  chaotic  motions. 

To  the  first-order  approximation  the  response  of  the  beam  is  given  by 
v(s,t)  =  vjsja^t)  cos(£it  -  y,)  +  ...  (1) 

w(s,t)  =  wn(s)a2(t)  cos(£2t  -  y2)  +  ...  (2) 

where  vm(s)  and  w„(s)  are  the  linear  mode  shapes  of  the  beam,  £2  is  the  excitation 


frequency,  and 

2oj1ma1'  +  [Rt  +  R2a2  sin  2 (y2  -  yjfla,  =  0  (3) 

Oimyi'  ~  R3  —  ^4ai  +  Rsa2  +  R2a2  cos  2(y1  —  y2)]a1  =  0  (4) 

2<w2na2'  +  [E1  +  E2ai  sin  2^  -  y2)]a2  -  /?0fou2n  sin  y2  =  0  (5) 

[2<y2ny2'  -  E3  -  E4a2  -  E5a*  -  E2a*  cos  2(y,  -  y2)]a2  -  fit fw2n  cos  y2  =  0  (6) 


Here,  o)lm  and  o)2n  are  the  linear  undamped  frequencies  of  the  interacting  modes,  f  is 
a  measure  of  the  excitation  amplitude,  an£j  the  /?„  R„  and  E,  are  constants  that  are 
given  in  quadrature. 

Periodic  solutions  of  the  beam  correspond  to  the  fixed  points  (i.e.,  constant 
solutions)  of  (3)-(6).  which  in  turn  correspond  to  a/  =  a/  =  y,'  =  y2'  =  0.  Thus,  (3)-(6) 
reduce  to  a  set  of  algebraic  equations  that  can  be  solved  numerically. 
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Figures  2  and  3  show  representative  frequency-response  curves  of  the  first 
and  second  modes.  The  results  show  that  planar  motions  in  the  z-direction  may  be 
destabilized  by  disturbances  in  the  y-direction,  resulting  In  vibrations  In  both  planes. 
For  certain  beam  and  load  parameters,  nonplanar  responses  can  be  either  (a)  un¬ 
steady  whirling  motions  (see  Fig.  4(a)-(c))  (lengths  and  directions  of  the  axes  of  the 
elliptical  route  continuously  change  and  the  motion  also  shows  a  low-frequency 
twisting  motion),  or  (b)  steady  whirling  motions  (see  Fig.  4(d))  (the  principal  axes  of 
the  elliptical  route  are  not  parallel  to  the  y-  or  z-axes),  or  (c)  chaotic  motions.  The 
unsteady  whirling  motions  may  be  periodic  if  the  period  of  the  amplitudes  is 
commensurable  with  that  of  the  base  motion.  But  this  periodicity  is  not  obvious  be¬ 
cause  the  twisting  frequency  is  too  low.  Figure  5  shows  the  frequency-response 
curves  of  a  beam  with  different  modes  in  two  principal  directions. 
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W(0,t)  =  F  cos  Qt 


Figure  1.  Coordinate  systems  :  x-y-z  =  inertial  reference  frame;  £  —  7  —  £  =  principal 
axes  of  the  beam's  cross  section  at  position  s,  which  are  fixed  on  the  cross  section. 
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Figure  2.  Response  curves  of  the  first  mode  for  a  beam  with  an  aspect  ratio 
b/h=s1.0  :  mode  (1,1)  ,  w„  =  (y2„  50  =  O.O  ,  <5*  =  -0.01  ,  ^  =  0.04,  /?, f  =  0.002, 
/?,  =  0.6489;  a?  =  planar  response  amplitude;  aln,aJn=  nonptanar  response  ampli¬ 
tudes;  ( _ )  stable,  ( — )  unstable  with  at  least  one  eigenvalue  being  positive. 
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Figure  3.  Response  curves  of  the  second  mode  for  a  beam  with  an  aspect  ratio 
b/h=t  1.0  :  mode  (2,2)  ,  cotl  =  coM,  3 „  =  0.0  ,  32  =  0.002  ,  /z  =  0.04,  /?9f  =  0.00006, 
/?,  =  0.6489;  a?=  planar  response  amplitudes;  aIn,aJn=  nonplanar  response  ampli¬ 
tudes;  ( _ )  stable,  ( — )  unstable  with  at  least  one  eigenvalue  being  positive,  (....) 

unstable  with  the  real  part  of  a  complex  conjugate  pair  of  eigenvalues  being  positive; 
(a)  response  curves  of  w(s,t),(b)  response  curves  of  v(s,t). 
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Abstract 


Rotating  frame  turbulence,  or  RFT,  refers  to  the  actual 
turbulence  experienced  by  the  blades,  and  requires  noneulerian 
description  and  sampling  of  measurements.  Its  characteristics  are 
strikingly  different  from  those  of  conventional  models  based  on 
eulerian  or  space-fixed  description,  and  these  differences  increase 
with  decreasing  advance  ratio  and  altitude.  Therefore,  RFT  effects  on 

gust  sensitivity  of  rotorcraft  and  tilt  rotors  are  appreciable  during 

low-altitude  and  low-speed  operations.  Typical  examples  include  NOE 
maneuvers,  light-house  and  offshore  oil-platform  missions  and  ship 
landing  of  rotorcraft.  The  tilt  rotor  "chugging"  problem,  which  is 
basically  a  fore-to-aft  low-frequency  acceleration  of  the  rigid  body 
mode  coupling  with  the  rotor  torque  mode,  is  another  example.  That 
chugging  occurs  at  relatively  low  speeds  and  aggravates  during 
descent  shows  that  RFT  effects  should  be  included  in  the  treatment 
of  this  problem. 

In  axial  flight,  though  the  RFT  and  space-fixed  turbulence  are 
stationary,  the  RFT  spectral  density  shows  large  peaks  at  integer 

multiples  of  the  rotor  speed  (IP,  2P,  etc.)  and  transfer  of  energy 

essentially  from  the  low-frequency  region  (<  IP)  to  the  high 
frequency  (>  IP).  Figure  1  shows  the  spectral  density  S(f)  of  the 
vertical  turbulence  velocity  at  0.7R  blade  station  on  the  basis  of  the 
Dryden  model.  The  results  refer  to  L/R=4,  very  low  altitude 
conditions  at  which  turbulence  scale  length  L  is  comparable  to  the 
rotor  radius  R.  They  are  presented  as  f*S(f)  versus  In  f,  which  gives  a 
better  picture  of  the  spectral  density  variation  or  transfer  of  energy 
in  the  low-frequency  region.  It  is  noteworthy  that  the  area  under  the 
curve,  ^f*S(f)d(ln  f),  gives  the  variance,  as  in  the  conventional  S(f) 
versus  f  presentation  shown  in  the  upper-left  inset  figure.  The 
upper-right  inset  figure  refers  to  the  autocovariance  function  R(x). 
Since  the  area  under  the  spectral  density  curve  is  the  same  in  both 
the  space-fixed  and  rotating  frame  description,  the  transfer  of 
energy  from  the  low-frequency  region  to  the  high-frequency  region 
is  balanced  by  the  occurrence  of  peaks.  The  RFT  effects  of  transfer  of 
energy  and  the  occurrence  of  peaks  cannot  be  predicted  in  the  space- 
fixed  description  of  turbulence  and  severely  affect  the  helicopter 
gust  response  statistics  of  standard  deviation,  threshold  crossing  and 


peak  distribution.  In  forward  flight,  the  differences  in  characteristics 


are  much  more  fundamental  in  that  the  RFT  becomes  periodically 


nonstationary  or  cyclostationary  although  the  space-fixed  turbulence 


remains  stationary.  Concerning  the  RFT's  stochastic  structure  and  its 


effects  on  helicopter  body  response,  this  paper  addresses  the 
following  topics:  1)  Instantaneous  or  frequency-time  spectrum  of 
RFT,  2)  Finite  impulse  response  (FIR)  filter  for  RFT,  and  3) 
Application  of  FIR  filter  to  generate  the  helicopter  body  response  to 
RFT,  including  comparisons  with  a  conventional  approach  in  which 
response  statistics  are  computed  directly. 

Instantaneous  Spectrum 

The  frequency-time  or  instantaneous  spectrum  Sw(f,t)  provides 
a  means  of  describing  both  the  energy  transfer  with  respect  to 
frequency  f  and  the  periodically  varying  nonstationarity 
(cyclostationarity)  with  respect  to  t.  Here  the  subscript  w  refers  to 
the  vertical  turbulence  velocity,  the  only  component  included  in  this 
investigation.  Moreover,  the  instantaneous  spectrum  Sw(f,t)  shows 
the  distribution  of  energy  in  the  (f,t)  plane  and  is  therefore  real.  By 
comparison,  the  conventional  double-frequency  spectrum  is  a 
complex  quantity  involving  a  real  part  (co-spectrum)  and  an 
imaginary  part  (quad  spectrum)  and  has  a  much  diminishing 
physical  meaning.  For  atypical  low-speed  (advance  ratio  p=0.05)  and 
low-altitude  (L/R=l)  condition,  Figures  2  and  3,  respectively,  show 
the  covariance  function  R(t,t)  and  the  instantaneous  spectral  density 
function  Sw(f,t).  It  is  highly  instructive  that  both  of  these  figures 
clearly  show  the  velocities.  It  is  this  cyclostationarity  that  has 
significant  bearing  on  the  development  of  shaping  filters  (driven  by 
white  noise)  and  the  treatment  of  helicopter  body  sensitivity  to 
turbulence. 


Finite  Impulse  Response  (FIR)  Filters 

In  axial  flight  the  RFT  is  stationary.  Therefore,  the  design  of 
shaping  filter  is  straightforward,  but  algebraically  cumbersome  due 
to  peaks  at  integer  multiples  of  the  rotor  frequency  (see  Figure  1). 
However,  in  forward  flight,  the  RFT  is  nonstationary  although  its 
nonstationarity  varies  periodically,  as  shown  in  Figures  2  and  3. 
Therefore,  the  shaping  filter  system  has  periodically  varying 
coefficients,  and  its  design  is  by  no  means  trivial.  The  FIR  filter  is  an 
intrinsically  stable  digital  filter  that  exploits  the  linearity  of  the 
problem  as  well  as  the  cyclostationarity  of  RFT.  Calculation  of  the  FIR 
filter  coefficients  involves  the  following  steps: 


a)  select  the  spectral  density  function  of  the  turbulence  for  which  a 
sample  function  is  to  be  created; 

b)  obtain  the  frequency  response  of  the  filter  using  the  spectral 
density  function  of  step  (a); 

c)  transform  the  frequency  response  into  impulse  response  by  a 
direct  use  of  the  inverse  Discrete  Fourier  Transforms; 

d)  window  the  impulse  response  to  the  desired  filter  length; 

e)  use  the  impulse  sequence  obtained  in  step  (d)  as  the  filter 
coefficients; 

f)  calculate  the  inner  product  of  the  filter  coefficients  obtained  in 
step  (e)  and  a  sequence  of  random  numbers  to  arrive  at  the  sample 
function  of  the  turbulence. 

Basically,  the  coefficients  of  the  FIR  filter  are  calculated  at  each 
instant  in  time  with  respect  to  the  instantaneous  power  spectrum 
Sw(f,t).  Because  the  RFT  is  periodically  nonstationary,  an  FIR  filter 
designed  using  this  approximation  will  have  periodic  coefficients  as 
well.  Conceptually,  this  amounts  to  approximating  the  nonstationary 
behavior  of  the  RFT  with  small  time  intervals  of  stationary  behavior. 

Body  Response  to  RFT 

An  attractive  feature  of  the  filter  approach  is  the  feasibility  of 
generating  the  response  sample  functions  to  RFT  excitations  and  the 
subsequent  evaluation  of  response  statistics  by  ’ensemble'  averaging. 
Therefore,  this  approach  can  be  directly  implemented  on  current 
flight  simulation  programs  such  as  GENHEL.  Figure  4  shows  a 
response  sample  function  of  the  UH-60A  Black  Hawk  helicopter  body 
normal  acceleration  response  to  turbulence  in  hover  using  the 
GENHEL  flight  simulation  program.  The  dominant  RFT  effects  are 
clearly  seen,  and,  as  expected,  these  effects  increase  with  decreasing 
values  of  turbulence  scale  length. 

Conclusions 

Major  findings  of  this  investigation: 

1.  The  instantaneous  or  frequency-time  spectrum  of  RFT 
provides  a  feasible  approach  to  modeling  RFT  to  predict  both 
the  cyclostationarity  (periodically  varying  nonstationarity)  and 
transfer  of  energy  from  the  low-frequency  (<1P)  to  the  high- 
frequency  (>  1 P)  region. 


2.  The  FIR  filter  approach  modified  to  treat  the  cyclostationarity 
RFT  can  be  directly  implemented  in  conventional  flight 
simulation  programs.  Since  it  generates  response  sample 
functions  that  are  subsequently  'ens<?  Me'  averaged,  this 
approach  can  be  used  to  verify  analytic...  _  computed  response 
statistics.  Preliminary  results  offer  considerable  promise. 

3.  RFT  effects  are  appreciable  on  helicopter  body  response.  This 
finding  has  a  significant  bearing  on  ride  quality  investigations 
and  future  development  of  ride  quality  specifications  in  the 
presence  of  turbulence. 
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Figure  1.  RFT  Effects  on  Space-fixed  Dryden  Model  for  the 
Vertical  turbulence  at  0.7R  blade  station. 
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Effect  of  Turbulence  on  the  Black  Hawk  Helicopter 
Normal  Acceleration  Response  for  LVR  =•  1. 
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Figure  k .  Effect  of  Turbulence  on  the  Black  Hawk  Helicopter 
Normal  Acceleration  Response  for  LVR  o  10. 
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ABSTRACT 

The  dynamics  of  a  rotor  blade  in  a  maneuvering  flight  condition  is  addressed  in  this  paper.  For 
simplicity,  a  uniform  rigid  rotor  blade,  with  small  chord  and  a  simple  flap  and  lag  spring  connection 
at  an  offset  hinge  is  considered.  The  aerodynamic  forces  are  modelled  by  quasi— steady  aerodynamic 
theory .  The  equations  were  derived  with  the  aid  of  the  symbolic  manipulation  package  MACSYMA  and 
a  number  of  fortran  programs  were  developed.  The  analysis  of  the  infinitesimally  small  perturbations 
about  a  trimmed  steady  state  solution  in  a  curved  maneuvering  flight  is  investigated  by  Floquet  theory. 
The  characteristic  multipliers  associated  with  the  perturbed  motion,  and  stability  boundaries  in  the 
two-parameter  space  of  flap  and  lag  spring  stiffness,  are  found  for  several  flight  parameters. 

INTRODUCTION 

To  date,  the  analysis  of  the  stability  of  the  motion  of  helicopter  rotor  blades  has  literally  been 
restricted  to  two  flight  conditions,  namely  hover  and  constant  velocity  forward  flight,  e.g.,  [1-5].  A  set 
of  differential  equations  of  motion  governing  the  dynamics  of  a  rotor  blade,  including  hub  acceleration, 
was  recently  presented  in  [6].  However,  to  the  authors’  knowledge,  no  comprehensive  analysis  dealing 
with  the  effect  of  maneuvers  on  the  blade’s  response  is  available.  This  work  addresses  this  void  in 
the  literature.  To  assess  the  influence  of  maneuvers  of  a  helicopter  on  the  stability  of  the  rotorblade 
motion,  the  simplest  possible  case  is  considered  first.  This  is  a  level  turn  in  steady  flight  with  constant 
speed  along  the  hub’s  trajectory,  resulting  in  a  constant  acceleration  of  the  hub.  For  steady  flight, 
the  homogeneous  differential  equatk  s  of  the  infinitesimal  motion  about  a  steady  state  solution  have 
periodic  coefficients  and,  thus,  stability  of  the  motion  can  be  assessed  by  making  use  of  well  known 
methods  such  as  Floquet  theory. 


EQUATIONS  OF  MOTION 

The  model  considered  in  the  analysis  consists  of  a  rigid  rotor  blade  with  length  R  and  hinge  offset 
e.  Fig.  1  shows  the  generalized  coordinates  0  and  (  and  feathering  pitch  angle  0  used  in  the  analysis. 
The  hub  is  assumed  to  be  moving  with  constant  speed  along  a  horizontal  circular  trajectory  of  radius 
Ro,  with  bank  angle  $  and  pitch  angle  0.  For  simplicity,  the  rotor  disk  force  vector  is  assumed  normal 
to  the  hub.  A  single  lead-lag  linear  spring  with  stiffness  K and  a  single  flap  linear  spring  with  stiffness 
I(p  is  considered.  Adding  up  all  the  relative  linear  and  angular  velocities  and  transforming  to  blade 
coordinate  axes,  the  absolute  velocity  of  a  point  on  the  blade  and  thus,  the  kinetic  energy  of  the  blade 
is  found.  Similarly,  from  coordinate  transformations,  the  angle  between  blade  chord  and  relative  air 
velocity  is  found,  from  which  the  generalized  aerodynamic  forces  Qp  and  are  obtained. 
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Figure  1.  Rotor  blade  coordinate  systems 


The  differential  equations  of  motion  are  obtained  from  Lagrange’s  equations  as 

If  -  I  +  K'V-M  =  *' 

If  '  f  +  *c(C-<w)  =  <?c 

producing  two  non-linear  equations  which  are  written  in  vector  form  as 

z(Ma  (.«,#)=  a  (2) 


and  the  constraint  equation 


9  (ip)  =  do  +  9ccos  (ip)  +  9ssin  (ip) 


(3) 
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The  rotordisk  forces,  obtained  from  aerodynamic  forces  on  the  blade,  are  written  in  normalized 
form  as 


2x 

°T  =  'k 

0 

2tt 

=  (4) 

0 

2?r 

Cy  =  ^7r/  C,C,0)# 

0 

Equations  (4)  and  the  expressions  for  the  functions  /r,  /;/  and  /y  were  obtained  with  MACSYMA[7]. 
They  are  too  lengthy  and,  due  to  lack  of  space,  cannot  be  reproduced  here. 

A  constant  helicopter  weight  is  assumed,  and  the  bank  angle  $  is  determined  from  radial  flight 
path  equilibrium.  To  eliminate  the  additional  parameters  of  the  position  of  the  center  of  gravity  of  the 
fuselage  in  this  study,  the  moment  trim  equations  are  disregarded,  and  the  pitch  angle  0  is  determined 
from  equilibrium  with  an  assumed  fuselage  drag  coefficient.  The  required  thrust  coefficient  Ct  is  then 
defined.  Trim  conditions  are  then 


Ct  :  given  ;  C//  =  0  ;  Cy  =  0  (5) 

These  replace  the  constraint  equation  (3)  so  that  do  ,  dc  and  ds  become  three  quantities  dependent  on 
Ct,  Ch  and  Cy. 


SOLUTION  OF  THE  EQUATIONS  AND  STABILITY 

Defining  a  vector  q(rp)  as  q(ip)  =  [P(ip),C(ip)]T,  the  solution  to  the  equations  of  motion  is  written 
as  q(ip)  =  £p(V’)+£a(V’)*  where  qp(ip)  is  the  particular  (trim)  solution,  and  qj(ip)  is  a  small  perturbation. 
The  particular  solution  is  approximated  by  a  finite  Fourier  series 


?p  = 


A)  +  Pucos  (ip)  +  f3usin  (ip)  + 
Co  +  Ciccos  (ip)  +  Cl  ssin  (ip)  + 


+  Pn/2<z>s  {Nip 1 2) 
+  Cjv/2Cos  (Nip/ 2) 


(6) 


The  solution  is  obtained  numerically  by  introducing  the  Fourier  series  for  q  into  equations  (2)  and 
(4).  The  resulting  2N+3  non-linear  algebraic  equations  for  the  2N+3  coefficients  fio,  ,  /?n/2,  Co,  . 
•  •  >  Cn/2,  t'o.  0c,  and  9S  are  solved  numerically  by  Newton-Raphson’s  method  using  exact  analytical 
expressions  for  derivatives  in  order  to  ensure  good  convergence. 

The  differential  equations  of  motion  for  the  infinitesimally  small  disturbance  q  (ip)  are  obtained 
by  linearizing  the  non-linear  equation  (2)  about  its  particular  solution  q  (ip).  The  linearized  equations 
are  of  the  form 


M  (gp)  9s +C  (V)  rjs+K  y  qs  =  0 


(7) 
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where  M  =  3*,  C  =  -g*  and  K  =  .  This  is  a  set  of  two  linear  differential  equations  with  periodic 

coefficients  determined  by  the  particular  solution  qp(i>)  and  by  the  periodic  pilot’s  input  0(^>). 

Converting  the  above  equation  into  a  system  of  four  equations  written  in  state  variable  form  as 
x  =  A(qp)x,  where  x  =  (/?,/?,  (,C]T,  Floquet  theory  is  then  used  to  obtain  the  characteristic  multipliers 
associated  with  the  linearized  equations,  and  to  determine  the  influence  of  a  number  of  system  parameters 
in  the  stability  of  the  response.  An  eigenvalue  routine  is  used  to  determine  the  four  complex  characteristic 
multipliers  zj,  (i=l, .  .  4),  from  which  four  complex  characteristic  exponents  s,  =  log can 

be  generated.  To  find  the  stability  boundary  in  the  parameter  plane,  these  are  regarded  as  functions  of 
the  parameters  (Kp,I(^).  Taking  the  smallest  u  =  Real(logzt)  for  a  chosen  point  in  the  parameter 
plane  (Jfy,  Ify),  the  particular  solution  xp(if>)  and  the  scalar  quantity  u(Kp,  K^)  are  calculated.  To  find 
the  stability  boundary  u(Kp,  K^)  =  0  in  the  parameter  plane,  we  start  with  the  solution  at  some  point, 

then  construct  the  gradient  grad(u)  =  by  computing  the  solution  at  small  increments 

dKp ,  dK(.  A  finite  step  AK  =  A  +  A is  then  used  to  “march”  in  the  parameter  plane  along 
the  gradient  line  to  find  one  point  on  the  boundary.  From  this  point  on  we  proceeded  in  chosen  steps 
AI<  along  the  tangent  direction,  normal  to  the  gradient  direction  (either  to  the  right  or  left  of  the  gradient 
line),  of  the  boundary  u  (Kp,  K^)  =  0.  A  final  Newton-Raphson  iteration  at  the  end  of  each  step  is 
used  to  obtain  an  accurate  solution  along  the  gradient  line. 

Due  to  lack  of  space  here,  results  will  be  shown  at  the  presentation  of  this  work. 
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Control  and  vibration  reduction  through  the  use  of  Intelligent  actuator  elements  has 
been  shown  to  be  effective  on  several  types  of  aerospace  structures.  Applications  include 
vibration  suppression  In  space  trusses  and  dynamic  control  of  camber  and  twist  for  gust 
alleviation  and  flutter  suppression  on  fixed  wing  surfaces  (1-3).  Vibration  suppression  In 
rotorcraft  could  be  enhanced  through  the  use  of  Intelligent  actuators  because  current 
methods  of  vibration  reduction  In  rotorcraft  do  not  address  some  of  the  vibration 
inducing  phenomena  that  occur  In  actual  helicopters  Including  differences  In  individual 
blade  tracking  and  magnitude  and  locations  of  dynamic  stall  (4-6).  Because  the  unsteady 
bending  moments  In  the  rotor  blade  are  several  orders  of  magnitude  greater  than  present 
intelligent  actuators  can  Impart,  the  direct  manipulation  of  the  rotor  blade  in  bending  is 
currently  not  feasible.  However,  through  blade  twist  manipulation  many  types  of 
vibration  reduction  methods  could  be  employed  including  suppression  of  blade 
vibrational  modes,  in-flight  tracking  and  dynamic  stall  reduction  through  small 
amplitude  pitch  oscillations  (7,8). 

Implementation  of  dynamic  blade  twist  requires  the  introduction  of  torsional 
forcing  to  quasi-orthotroplc  blade  structures  like  uncoupled  composite  or  aluminum 
blades.  Since  nearly  all  composite  rotor  blades  In  use  today  have  such  characteristics,  a 
method  of  torsional  forcing  must  be  developed.  Many  types  of  intelligent  actuation 
devices  (including  piezoelectric  actuators)  in  production  today  are  incapable  of  imparting 
this  torsional  forcing  due  to  their  quasi-isotropic  nature.  Accordingly,  this  paper  deals 
with  the  analytical  and  experimental  development  of  torsional  forcing  methodology  as 
it  relates  to  the  manipulation  of  Intelligent  structures. 

To  introduce  torsional  forcing  from  a  quasi-isotropic  material,  a  method  of 
directional  attachment  was  developed.  The  purpose  of  the  directional  attachment  is  to 
maximize  the  equivalent  longitudinal  stiffness  of  the  piezoelectric  crystal  while 
minimizing  the  effective  transverse  stiffness  (as  experienced  by  the  laminate).  One  type 
of  directional  attachment  is  shown  In  Figure  1. 


Figure  1  Directional  Attachment  of  Piezoelectric  Crystal  to  Substrate 

Two  types  of  test  specimens  were  fabricated  to  demonstrate  this  principle.  A  bending 
specimen  with  the  arrangement  of  Figure  2  demonstrated  that  significant  reductions  of 
transverse  stiffness  could  be  achieved  without  corresponding  reductions  In  longitudinal 
stiffness.  The  crystals  and  the  attachment  direction  were  oriented  longitudinally  and 
transversely  for  experimentation. 


Figure  2  Directional  Attachment  Specimen  Configuration 

The  other  plate  test  article  had  crystals  at  ±45°  angles  from  the  longitudinal  axis. 
The  crystals  produced  torsional  deflections  of  the  aluminum  substrate  which  further 
demonstrated  that  anisotropic  behavior  could  be  achieved  by  directionally  attaching  two 
Isotropic  materials.  Piezoelectric  crystals  were  then  Imbedded  in  the  skin  of  a  Froude 
scaled  rotor  blade.  From  the  experiments  that  followed,  the  feasibility  of  using 
directional  attachment  to  introduce  significant  torsional  deflections  of  aerodynamic 
surfaces  was  demonstrated. 


The  analytical  treatment  of  piezoelectric  crystals  in  a  laminate  takes  the  form  of 
strain-energy  relationships.  For  the  static  actuation  case,  the  strain  energy  in  a  given 
cross-section  is  set  equal  to  zero  and  the  laminate  strains  are  solved  as  can  be  seen  as 
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Where  the  subscripts  "a';  and  "1"  denote  the  piezoelectric  actuator  and  laminate  and 
the  integration  variable  "x"  is  in  the  longitudinal  direction  of  the  laminate.  Expanding 
equation  1  to  show  the  components  of  the  stiffness  and  actuation  matrices  gives  greater 
insight  to  the  problem. 
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For  piezoelectric  crystals,  the  number  of  actuation  strains  reduce  to  only  ex  and  £y  as 
the  other  types  of  strain  actuation  including  shear  and  bending  are  negligible.  It  should 
be  noted  that  one  inherent  characteristic  of  piezoelectric  crystals  is  that  the  amount  of 
transverse  strain  is  usually  lower  than  the  amount  of  longitudinal  strain  (9).  Assuming 
completely  isotropic  behavior  of  the  substrate  and  the  actuator  as  well  as  the  same 
Poisson's  ratios,  the  solution  of  the  beam-bending  case  is  seen  for  the  configuration  of 
Figure  2. 
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Predictions  of  the  beam  performance  followed  closely  to  the  predicted  performance  of 
equation  3.  The  modeling  of  the  directional  attachment  requires  closer  attention  to  the 
boundary  conditions  on  each  segment  of  the  crystal  and  finite  bond  thickness  also 
affects  the  amount  of  directionality  that  the  crystal  can  impart  to  the  structure.  With 
greater  crystal  aspect  ratios,  the  directionality  increases  due  to  the  transverse  shear  lag 
that  the  attachment  area  encounters. 

Future  work  on  intelligent  structures  for  rotorcraft  will  incorporate  more  involved 
lamination  schemes  for  the  integration  of  the  actuators  into  the  rotor  blade  structure. 
The  Froude  scaled  rotor  blade  that  was  constructed  for  this  investigation  demonstrated 
significant  torsional  deflections  for  the  purposes  of  some  types  of  vibration  control. 
However,  the  Froude  scaled  blade  used  non-conducting  composite  plys  around  the 
crystals  (fiberglass  and  epoxy).  This  facilitated  greater  transfer  of  the  crystal  actuation 
shear  strains  due  to  the  extremely  thin  bond  line  (1-2  mil  between  crystal  surface  and 
fibers).  Since  many  types  of  composites  are  good  conductors  (graphite  and  boron),  an 
insulation  layer  must  be  present  between  the  crystal  and  the  laminate  to  prevent 
shorting  of  the  crystal  [10).  For  directionally  attached  crystals,  the  insulation  layer  may 
induce  unacceptably  high  amounts  of  shear  lag  and  cause  the  structures  to  have  low 
amounts  of  deflection.  Further  investigations  of  directionally  attached  crystals  in 
Froude  scaled  rotor  blades  will  determine  these  characteristics. 
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The  aeromechanical  instability  of  a  helicopter,  on  the  ground  or 
in  flight,  is  caused  by  coupling  between  rotor  and  body  degrees 
of  freedom.  This  instability  is  commonly  denoted  air  resonance 
when  the  helicopter  is  in  flight  and  ground  resonance  when  the 
helicopter  is  on  ground.  The  phenomenon  is  one  of  the  classical 
problems  of  helicopters  known  since  the  early  days  of  helicopter 
history,  but  still  today  it  is  a  very  actual  one,  especially  for 
helicopters  with  hingeless  or  bearingless  rotors  which  are  sus¬ 
ceptible  to  such  instabilities. 

There  is  a  good  basic  understanding  of  the  problem  and  how  to 
avoid  the  instability  in  general.  But  for  bearingless  rotors  it 
is  often  very  difficult  to  provide  the  required  damping.  Therefo¬ 
re  the  possibilities  of  active  control  in  the  sense  of  artificial 
stabilization  of  ground  or  air  resonance  find  increasingly  inte¬ 
rest.  Such  an  active  controller  would  need  sensing  and  actuating 
devices,  which  can  operate  with  a  conventional  swashplate,  and 
the  hardware  realization  should  be  no  basic  problem. 

The  aim  of  the  paper  is  to  discuss  the  physical  understanding  of 
the  phenomenon,  for  which  some  modeling  aspects  seem  to  be  very 
important.  Special  consideration  will  be  given  to  possible  con¬ 
trol  schemes  working  through  the  full  range  of  operating  condi¬ 
tions  of  the  helicopter  which  can  normally  be  encountered.  Espe¬ 
cially  stabilization  of  the  air  resonance  instability  has  to  be 
seen  in  direct  relation  to  stability  and  control  angmentation  of 
the  whole  helicopter.  Adverse  effects  on  the  flight  mechanical 
modes  and  thereby  on  the  handling  qualities  have  to  be  avoided. 
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Abstract 

A  direct  analysis  of  the  response  of  a  rotor  blade  is  presented  without  having  to  replace  the 
full  nonlinear  equations  of  motion  of  the  system  by  a  set  of  equations  expanded  into  polynomial 
nonlinearities.  After  the  equilibrium  solution  of  the  full  nonlinear  equations  is  determined, 
the  perturbed  motion  about  the  equilibrium  is  then  analyzed.  Infinitesimally  small  motions 
are  considered,  and  both  the  eigenfunctions  and  the  eigenvalues  associated  with  linearized 
motions  are  determined  by  a  direct  method  which  is  based  on  the  properties  of  the  transition 
matrix  associated  with  such  motion. 


Introduction 

The  investigation  of  the  response  of  a  rotor  blade  involves  the  determination  of  a  par¬ 
ticular  solution  to  the  differential  equations  of  motion  and  the  eigenvalues  associated  with 
infinitesimally  small  motions  about  such  solution.  Because  of  the  complexity  of  the  governing 
differential  equations  of  motion,  it  is  common  practice,  in  the  rotorcraft  dynamics  literature, 
to  first  expand  such  equations,  to  a  pre-determined  degree,  about  the  undeformed  state  of 
the  blade.  The  resulting  equations,  which  contain  polynomial  nonlinearities  truncated  to  the 
desired  order,  are  then  used  to  analyze  the  aeroelastic  response  of  the  blade.  One  must  be 
careful  when  working  with  expanded  equations  in  order  to  obtain  results  that  are  mathemati¬ 
cally  consistent  with  the  approximating  assumptions  and  which  do  not  produce  results  that  are 
physically  invalid  such  as  those  presented  in  [1]  with  quadratic  equations.  As  shown  in  [2-4], 
the  inclusion  of  cubic  nonlineariries  is  necessary  in  order  to  circumvent  such  inconsistencies. 

An  analysis  of  the  response  of  a  rotor  blade  in  hover,  including  the  effects  of  cubic  nonlin¬ 
earities,  was  presented  in  [5,6]  making  use  of  Galerkin’s  method  with  a  set  of  eigenfunctions 
of  a  non-rotating  beam.  The  number  of  Galerkin  coefficients  is  dependent  on  the  number,  N, 
of  eigenfunctions  and  on  the  order  of  truncation  of  the  differential  equations  of  motion.  As 
shown  in  [6],  such  number  is  greatly  increased  when  one  increases  the  truncation  order  from 
quadratic  to  cubic.  Also,  as  shown  in  [7]  for  a  cantilever  with  a  large  tip  mass,  even  the  cubic 
appxoximation  for  the  equations  may  yield,  in  some  cases,  very  innacurate  eigenvalues  for  the 
linearized  motion. 

In  the  present  work,  the  problem  of  determining  the  equilibrium  solution  and  the  eigenval¬ 
ues  associated  with  the  perturbed,  infinitesimally  small,  flap-lag-torsional  aeroelastic  response 
of  a  rotor  blade  in  hover  is  reformulated.  The  equilibrium  state  is  determined  directly  by  solv¬ 
ing  the  two-point  boundary  value  problem  governed  by  the  original,  full-nonlinear,  equations 


of  motion.  A  transfer  matrix  approach,  based  on  the  properties  of  the  transition  matrix  of 
the  system,  is  used  for  determining  the  eigenfunctions  and  the  eigenvalues  associated  with  the 
perturbed  motion  of  the  system. 


Analysis 

The  differential  equations  governing  the  flexural-flexural-torsional  motion  of  inexten- 
sional  beams,  taking  into  account  all  the  geometric  nonlinearities  in  the  system,  were  for¬ 
mulated  in  [2],  and  in  [3]  for  both  extensional  and  inextensional  beams.  The  system  con¬ 
sidered  here  consists  of  an  initially  straight  and  untwisted  cantilevered  rotor  blade  of  un¬ 
deformed  length  R,  pre-cone  angle  0  and  collective  pitch  #o  [5].  The  differential  equations 
of  motion  for  the  blade  were  formulated  using  the  approach  followed  in  [3],  with  the  aero¬ 
dynamic  forces  modelled  as  in  [5].  For  convenience,  they  are  written  in  terms  of  orien¬ 
tation  tingles  (&z,Oy,0x),  which  are  related  to  the  elastic  deformations  (normalized  by  R) 
w(x,  t ),  u(x,  t), u(x,  t )  shown  in  Fig.l,  and  the  elongation  eo  =  \/((l  +  u')2  +  v12  +  to'2)  —  1  as 
v1  =  (1  +  eo)cos0j,  sin#-,u>'  =  (1  +  e0)sin^,  1  +  u1  =  (1  +  e0)cos#j,  cos 6Z.  The  blade’s  angle 
of  twist,  <f>(x,  f),  is  given  as  <j>(x ,  t)  =  9X  +  fQ  0z(jj,  t )  sin  #5,(77,  i)  drj.  Here,  primes  denote  partial 
differentiation  with  respect  to  the  spatial  variable  2,  which  is  the  distance  along  the  unde¬ 
formed  blade’s  reference  axis,  normalized  by  R,  and  dots  denote  partial  differentiation  with 
respect  to  time,  normalized  by  the  rotational  angular  speed  f2.  Considering,  for  simplicity,  a 
homogeneous  blade,  the  normalized  differential  equations  of  motion  are 

- —  +  tan#z  f  fudx]  =  v  -f  2ri  cos  0  -  2w  sin  0  -  v  -  Qv 

\(1  +  e0 )  cos  cos  #.  J\  ) 

G'w  =  (  Gv  sin  9,  tan  8y  +  —— — ^ - —  +  tan  9y  cos  6Z  [  fu  dx 

\  y  (1  +  e0)cos#y  y 

—w  +  2v  sin  0  +  (x  +  u)(sin  20)/2-w  sin2  0  —  Qw 
Aex  =Qez 


where  fu  =  il  -  2v  cos  0  +  w(sin20)/2  -  (2  +  u)cos2  0  -  Qu  -,Qu,Qv,Qw  and  Q6x  are  the 
generalized  aerodynamic  forces,  whose  virtual  work  is  QuSu  +  QvSv  +  Qw9w  +  Qex89x.  The 
expressions  for  Agt,Ae  and  Aex  are  obtained  as  shown  in  [2,3]. 

To  analyze  the  motion,  Eqs.  (1)  are  written,  for  convenience,  in  state  variable  form  as 
y'(x,t)  =  f(y,y,ysi/,y)-  The  state  equations  admit  the  equilibrium  solution  y  =  y  (x). 
The  static  equilibrium  solution  y£(x)  was  determined  numerically  by  solving  the  two-point 
boundary- value  problem  governed  by  the  resulting  ordinary  differential  equation  y'  (x)  = 
/(ye(x)).  By  perturbing  the  equilibrium  solution  as  y  =  ye  +  ey^(x,t),  expanding  the  state 
equations  about  e  =  0,  and  truncating  the  expansion  to  0(e) ,  a  set  of  linearized  differen¬ 
tial  equations  are  obtained.  With  y^x,*)  =  F(x)exp(rt),  these  equations  are  of  the  form 

F'(x)  =  A(r)Y_(x).  The  eigenfunctions  F(x)  and  complex  eigenvalues  r  =  rr  +  ^  —  l)r,  were 
determined  numerically  by  a  transfer  matrix  method. 
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<£,(  S=1)  wA(s-1) 


Fig.  1  (a  and  b)  show  plots  of  the  static  equilibrium  solution  at  the  blade’s  tip,  ve(l),  tue(l) 
and  <^e(l),  and  one  of  the  first  three  eigenvalues  associated  with  the  perturbed  motion,  ver¬ 
sus  u*w  —  the  blade’s  non-dimensional  rotating  flap-frequency  in  vacuum  for  #0=0  [8]. 
The  results  shown  in  these  figures  are  for  the  four-bladed  rotor  considered  in  [6],  with  non- 
dimensional  rotating  lead-lag  and  torsional  frequencies  in  vacuum, equal  to  u;*  =  1.5  and 
=  2.5,  respectively.  The  values  chosen  for  other  parameters  of  the  problem  are  those 
used  in  [6].  The  equilibrium  solution  and  the  eigenvalues,  for  #0  =  0.5,  determined  in  [6]  by 
a  Galerkin’s  procedure  with  five  non-rotating  beam  modes  are  also  shown  in  these  figures 
(points  marked  x )  for  comparison  with  the  more  accurate  results  determined  by  the  method 
presented  here.  Comparison  of  the  results  of  the  two  analyses  disclose  the  effect  of  order 
truncation  when  the  full  nonlinear  differential  equations  of  motion  of  the  system  are  expanded 
about  a  state  that  is  not  an  equilibrium,  and  replaced  by  the  resulting  set  of  approximate 
equations  with  polynomial  nonlinearities.  These  conclusions  are  consistent  with  the  results 
presented  in  [7]  for  a  cantilever  with  a  large  tip  mass. 


Fig.  1.  Static  equilibrium  (ve,we,<f>t),  at  blade’s  tip,  and  real  (cr^)  and  imaginary  (u^)  parts  of 
an  eigenvalue  versus 
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ABSTRACT 

Many  problems  arising  in  f.he  dynamics  of  rotating  systems  such  as  the 
helicopter  blades,  the  rotor  bearing  systems  and  the  dynamic  stability  of 
structures  lead  to  a  set  of  differential  equations  with  periodic  coefficients. 
In  general,  the  equations  are  nonlinear,  however,  the  stability  analyses  of  all 
such  systems  require  investigation  of  linear  equations  only.  In  spite  of  the 
availability  of  many  '-echniques,  Floquet  analysis  has  served  as  the  main  tool 
in  various  applications.  Floquet  analysis  is  a  powerful  technique  which  is 
easily  implemented  on  a  computer  and  holds  most  promise  in  the  analysis  of  large- 
scale  periodic  systems.  Evaluatj  \  of  'Floquet  Transition  Matrix' (FTM)  and 
eigenannlysis  are  the  two  major  c  imputational  problems  encountered  in  the 
analysis  of  periodic  systems.  Although  many  techniques  have  been  suggested  in 
the  literature,  the  computation  of  FTM  for  large-scale  system  is  still  a 
challenging  task.  At  the  present  time,  the  application  of  a  higher  order  initial - 
value  codes(such  as  the  fourth-order  Runge-Kutta  or  Hamming's  Predictor- 
Corrector)  in  a  'single  pass'  scheme  is  regarded  as  the  most  efficient  algorithm 
in  the  evaluation  of  FTM.  The  'single  pass'  scheme  was  suggested  about  fifteen 
years  ago,  but  since  then  not  much  progress  has  been  made  toward  the  development 
of  new  computational  techniques. 

In  this  study,  a  direct  numerical  evaluation  of  Floquet  Transition  Matrix  is 
proposed.  The  solution  of  periodic  systems  is  represented  in  terms  of  shifted 
Chebyshev  polynomials  of  first  or  second  kind.  The  main  feature  of  this  technique 
is  that  it  reduces  the  original  differential  problem  to  a  system  of  linear 
algebraic  equations. 

A  multidimensional  linear  periodic  system  may  be  represented  in  the  state- 
space  form  as 

y(t)  -  (A(t)J  y(t)  ;  (A(t+T)]  -  (A(t)]  ,  y(0)  -  y0  (1) 

where  y(t)  is  n  vectors  and  (A(t)j  is  an  n  x  n  periodic  matrix  with  period  T. 
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For  a  set  of  Chebyshev  polynomials,  S(t)  defined  over  the  period  [0,T],  the 
vectors  and  matrices  appearing  in  equation  (1)  may  be  expanded  as 


m-1 

y(t)  -  2  CjSj(t)  (2) 

i-0 

m-1 

A(t)  -  2  DjSjCt) 

i-0 

and 

y(0)  -  y0S0( t )  ,  (3) 


where  Cj  are  the  unknown  vectors,  yg  is  the  initial  vector,  and  Dj  are  known 
vectors . 

Integrating  equation  (1),  utilizing  relations  in  equations  (2)  and  (3),  and 
equating  the  coefficients  of  Sj  leads  to  a  set  of  linear  algebraic  equations 
in  Cj ,  which  can  be  solved  by  a  standard  technique. 

Once  Cj's  are  known,  the  solution  y(t)  can  be  obtained  from  equation  (2),  and 
the  FTM  can  be  directly  evaluated. 

Several  illustrative  examples  including  the  stability  of  Mathieu's  equation 
and  column  under  periodic  load  are  discussed.  Potential  applications  to 
rotorcraft  stability  problem  are  also  indicated. 
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There  are  several  probabilistic  methods  which  have  been  presented 
in  the  literature  for  investigating  the  stability  of  stochastic  dynamical 
systems.  Most  investigations  have  been  concerned  with  stability  in 
moments  and  almost  sure  or  sample  stability.  In  practice,  it  is  generally 
most  desirable  to  examine  the  almost-sure  sample  stability  and  the  results 
thus  obtained  hold  true  with  probability  one.  It  is  clear  that  for  first  order 
Ito  equations,  the  region  of  stability  for  higher  integer  moments  are 
included  in  the  regions  of  stability  of  lower  integer  moments  and 
furthermore,  all  integer  moment  stability  regions  are  included  in  the 
region  of  sample  stability.  Recently,  an  exact  relation  between  conditions 
for  sample  stability  and  conditions  of  moment  stability  were  established  in 
Kozin  and  Sugimoto  [1]  and  extended  by  Arnold  [2].  It  has  been  shown  by 
Kozin  and  Sugimoto  [1]  that  sample  stability  region  in  some  parameter 
space  is  a  limit  of  the  regions  of  the  pth  moment  stability  as  p  i  0.  However, 
it  is  difficult  to  obtain  sample  stability  in  such  a  fashion,  since  the 
conditions  for  the  pth  moment  stability  is  more  difficult  to  obtain  for 
arbitrary  p,  in  multidimensional  multiplicative  systems.  For  this  reason, 
sample  stability  in  this  paper  is  obtained  using  a  very  well  developed 
concept  of  Lyapunov  characteristic  exponents  (see,  e.g.,  Khasminskii  [3] 
and  Arnold  and  Klieman  [4]).  Furthermore,  it  is  well  known  that  in 
regions  where  the  system  is  stable  in  probability  one,  the  second  moments 
may  grow  exponentially.  It  is,  therefore,  important  from  the  design  point  of 
view  that  the  higher  order  moments  decay  in  some  parameter  region.  To 
this  end,  mean  square  stability  conditions  are  also  obtained  in  this  work. 
The  theory  developed  here  is  applicable  to  many  engineering  systems, 


however,  in  this  paper  attention  is  focussed  on  helicopter  dynamics. 


When  a  helicopter  is  in  forward  flight,  even  under  ideal  smooth  flow 


assumption,  the  blade  encounters  different  flow  regions  in  a  period  of 


revolution  and  the  aerodynamic  loads  along  the  blade  are  distributed 
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nonuniformly.  By  virtue  of  this  complex  but  periodic  aerodynamic  loading, 
a  helicopter  rotor  blade  is  a  time  variant  system  and  is  susceptible  to  severe 
vibrations.  Noteworthy  investigations  of  rotor  dynamics  in  forward  flight 
have  been  made  by  Stammers  [5],  Peters  [6],  Friedman  [7],  et  al.  Stammers 
paper  discussed  the  nature  of  flapping  torsion  flutter  and  dealt  with  the 
presence  of  periodic  coefficients  in  the  equations  of  motion.  It  was  shown 
that  forward  flight  can  have  a  significant  stabilizing  influence  on  flutter 
and  that,  as  tip  speed  increases,  flutter  occurs  predominantly  at  half¬ 
integer  frequencies.  Peters  [6]  considered  the  problem  of  flap-lag  stability  in 
forward  flight.  It  was  shown  that  periodic  coefficients  were  important  for 
flap-lag  stability  even  at  low  advance  ratios.  It  was  also  concluded  that  the 
variation  in  rotor  trim  and  in-flow  with  forward  speed  had  a  significant 
effect  on  blade  stability. 

In  the  service  life  of  a  helicopter,  numerous  encounters  with 
thunderstorm  and  clear  air  turbulence  can  be  expected.  Because  of  the  very 
nature  that  the  lift  is  generated  by  blade  rotation,  some  level  of  self-created 
turbulence  is  also  unavoidable.  Thus,  to  make  the  analysis  realistic, 
random  turbulence  in  the  atmosphere  is  included  in  many  recent 
investigations,  and  considerable  progress  has  been  made  towards 
understanding  mean  square  flap-lag  stability.  The  inquiry  into  the  effects 
of  atmospheric  turbulence  on  the  stability  of  a  rotor  blade  was  initiated  by 
Lin,  Fujimori,  and  Ariaratnam  [8].  Here,  the  stability  of  uncoupled  flap 
motion  in  the  presence  of  random  turbulence  was  investigated  for  the 
operating  condition  of  hover  and  forward  flight. 

The  coupled  flap-lag  equations  of  motion  for  forward  flight  previously 
derived  by  Peters  [6]  were  generalized  by  Prussing  and  Lin  [9]  to  include 
random  turbulence  in  the  airspeed  components  of  the  blade.  Here,  a  mean 
square  stability  analysis  was  performed  only  for  the  case  of  hover  with 
vertical  turbulence.  It  was  shown  that  vertical  turbulence  appears  to 
increase  the  stability  of  the  coupled  flap-lag  motion  for  realistic  rms 
turbulence  velocities. 

In  this  work,  almost  sure  stability  for  coupled  flap-lag  motion  is 
examined  using  the  well  developed  concept  of  Lyapunov  exponents.  To  this 
end,  the  equations  of  motion  derived  by  Prussing  and  Lin  [9]  are  used,  i.e,. 
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where 


cs  =  Cjj  +  5  C§  +n  cj  +  «C8 

Kj  =  iqi  +  §K5+nKj  +  t>K^ 

Fi  =  5F?  +  TiP1n  +  uFi\  i  =  1,2, j  =  1,2. 

The  variables  p  and  £  are  the  blade  flap  and  lead-lag  angles, 
respectively,  and  {;,  T|  and  x>  are  the  non-dimensional  lateral,  longitudinal, 
and  vertical  random  turbulence  components  of  the  air  flow  relative  to  the 
blade.  The  matrices  C  and  K  consist  of  deterministic  (Cy,  Ky)  and  stochastic 
components  and  the  inhomogeneous  forcing  terms  Fi  and  F2  are  also 
functions  of  the  random  turbulence  components.  This  analysis  is  an 
extension  of  the  recent  work  reported  by  Sri  Namachchivaya  and  Prussing 
[10]  on  the  almost  sure  stability  of  decoupled  flap  motion. 

As  the  first  step  in  the  analysis,  a  set  of  approximate  Ito  equations 
are  derived  using  both  deterministic  and  extended  stochastic  averaging 
methods.  The  effect  of  turbulence  is  greater  near  the  critical  points  of  the 
dynamical  system.  In  the  vicinity  of  such  critical  points,  extended 
stochastic  averaging  yields  a  set  of  2-D  reduced  Ito  equations  as 


da  =  ma  (a,<{>,p,X)  dt  +  Gy  (a,|i,X)  dwjy 
d<f>  =  m^  (a,<|>,p,X)  dt  +  Gy  (a,<|),|i,^)  dw^- 

where  ma,  m^  are  the  drift  coefficients  and  Oy  and  Gy  are  the  diffusion 
coefficients  of  the  associated  Fokker-Planck  equation,  and  Waj,  W<pj  are 
Wiener  processes  of  unit  intensity.  It  is  worth  pointing  out  that  a  and  <{> 
processes  are  coupled  and  the  rate  of  growth  of  the  amplitude  process 
depends  on  the  evaluation  of  the  <j>  process  on  a  unit  circle.  The  Lyapunov 

/s 

exponent  X  defined  as 


4 


/ 

X  =  lim  -  In  1  I  x  (t;  x0,  t0)  I  I  ,  x  -  (dp,  8p} 

t— >  oo  t 


is  calculated  for  the  above  Ito  equations  using  the  approach  given  by 
Khasminskii  [3],  and  Sri  Namachchivaya  and  Prussing  [10].  The 
dynamical  system  is  said  to  be  almost  sure  stable  (stable  with  probability  1) 
•  if  X  <  0. 

Consider  an  n-degree  of  freedom  system,  which  represents  many 
physical  problems  including  the  flap-lag  motion,  i.e. 

xj  +  2ajXj  +  flf  xj  =  ey2AijXj^(t)  +  ep  BjjXj  sim)t,  ij  =  1,2 . ,n 

where  Ajj  and  Bjj  are  elements  of  constant  matrices.  In  this  equation,  we 
shall  assume  that  the  1st  mode  (lead-lag)  has  small  damping  compared  to 
the  rest  of  the  modes.  The  maximal  Lyapunov  exponent  X  can  be  calculated 
as 


X  = 


/lidM  +  IwaShH 
l  hW  I  ’ 


kjj  kji 
4kj 


where  Ii.j^(p)  represents  the  Bessel  function  of  imaginary  arguments,  and 
imaginary  order, 

if  -  ^ . .  a  =  j ; *  - + : kft  v<i>  - 1  ^  [v®5>  -  ¥»»>])  /  d 

D  =1  +  [kfx  (2S(0)  +  S(l))] ,  kj  =  oycoo ,  v  =  co0  (1  -  eX) 

Ja 

Py  =  Ai/0)2 ,  ky  =  By  CO2 


S  (s)  and  \j /(*)  are  cosine  and  sine  spectrums  of  the  stochastic  excitation 
£(t),  and  X  is  the  detuning  parameter. 
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Introduction  and  Background 

Currently  at  NASA  Langley  Research  Center,  there  is  an  effort  to  integrate  various 
disciplines  in  the  rotor  blade  design  process*.  An  investigation  at  partially  integrating  some  of 


1 


these  disciplines  is  reported  in  Ref.  2,  where  the  integration  of  aerodynamic  loads  effects  and  the 

dynamic  aspects  of  blade  design  was  addressed  by  coupling  a  comprehensive  helicopter  analysis 

3  4 

code,  CAMRAD  ,  the  nonlinear  optimization  algorithm,  CONMIN  ,  and  an  approximate 

analysis  technique.  A  combination  of  the  blade  root  4/rev  vertical  shear  and  the  blade  weight 

was  minimized  with  constraints  on  coupled  lead-lag  and  flapping  frequencies,  blade 

autorotational  inertia  and  centrifugal  stress.  The  paper  demonstrated  a  significant  reduction  in 

the  4/rev  vertical  shear  and  blade  weight,  which  were  objective  functions,  along  with  overall 

reductions  in  the  amplitudes  of  the  oscillatory  vertical  airloads  both  azimuthal  and  radial.  As  a 

byproduct,  it  was  shown  that  optimization  also  reduced  the  total  power  required  by  the  rotor 

while  maintaining  the  same  C^Jo,  C^,  being  the  rotor  thrust  coefficient  and  a  the  thrust  weighted 

solidity  of  the  blade. 

It  was  of  interest  to  extend  the  work  of  Ref.  2  by  including  other  sources  of  blade  vibration 
that  are  transmitted  to  the  hub.  For  a  four-blade  rotor  the  other  major  oscillatory  loads  are  the 
3/rev  and  5/rev  inplane  and  radial  shears,  the  3/rev  and  5/rev  flapping  and  torsional  moments  and 
the  4/rev  lagging  moment.  These  are  included  in  the  optimization  formulation,  either  in  the  form 
of  objective  functions  or  as  constraints,  in  the  present  work.  Also,  the  total  thrust  of  the 
optimized  blade  in  Ref.  2  was  lower  than  that  of  the  reference  blade  although  the  C^/a  was  held 

constant  during  optimization.  To  avoid  this  loss  in  the  rotor  thrust,  a  constraint  is  imposed  on  the 
total  thrust. 

Problem  Statement  and  Formulation 

Optimization  techniques  are  applied  to  minimize  the  vibratory  blade  loads  of  a  four-blade 
helicopter  in  forward  flight.  This  is  done  by  considering  all  the  major  sources  of  vibration 
(forces  and  moments)  in  the  rotating  frame  that  are  transmitted  to  the  nonrotating  hub  as  4/rev 
harmonics  of  the  longitudinal,  lateral  and  vertical  forces  and  the  pitch,  roll  and  yaw  moments. 

The  objective  functions  are  the  4/rev  vertical  shear  and  the  3/rev  inplane  shear.  Constraints  are 
imposed  on  the  3/rev  radial  shear,  the  3/rev  flapping  and  control  moments,  the  4/rev  lagging 
moment.  Constraints  are  also  imposed  on  natural  frequencies,  blade  weight,  centrifugal  stress 
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and  blade  autorotational  inertia.  A  constraint  is  incorporated  on  the  rotor  thrust  to  ensure  that  the 
optimum  blade  performs  the  same  rotor  task  as  the  reference  blade.  Side  constraints  are  also 
imposed  on  the  design  variables  to  avoid  impractical  solutions.  The  bounds  are  selected  from  an 
existing  advanced  rotor  blade  model  which  will  be  called  the  ’reference’  blade  hereafter. 

Blade  Model  and  Design  Variables 

A  linear  taper  is  allowed  along  the  blade  planform  (Fig.  1)  and  the  blade  taper  ratio,  X,  is  X 

=  c  /  c  ,  where  c  is  the  root  chord  and  c  is  the  tip  chord.  It  is  assumed  that  the  blade  stiffnesses 
r  t  r  t 

arise  solely  from  the  blade  structural  component  (i.e.,  the  contributions  of  the  nonstructural 
masses,  the  skin  and  honeycomb,  etc.  are  neglected. 

The  design  variables  are  the  blade  stiffnesses  at  the  root,  the  taper  ratio,  the  root  chord,  the 
radius  of  gyration  at  blade  root  and  the  nonstructural  weights  located  spanwise.  The  blade 
preassigned  properties  and  the  helicopter  performance  parameters  are  used  from  the  existing 
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reference  blade  data.  A  method  called  the  ’Global  Criteria  Approach’  is  used  to  formulate  the 
multiple  objective  optimization  problem. 

Analysis 

The  program  CAMRAD  is  used  for  the  blade  aerodynamic  and  dynamic  analysis  and  the 
program  CONMIN  is  used  for  the  optimization.  In  addition,  approximate  analyses  of  the 
objective  function  and  constraints  involving  Taylor  series  expansions  are  used  to  save 
computational  effort.  Sensitivity  analyses  for  computing  derivatives  of  the  objective  function 
and  the  constraints  are  part  of  the  procedure. 

The  blade  is  trimmed  at  each  step  of  design  optimization.  A  wind  tunnel  trim  option, 
provided  in  CAMRAD,  which  consists  of  trimming  the  rotor  lift  and  drag,  each  normalized  with 
respect  to  a,  and  the  flapping  angle  with  collective  and  cyclic  pitch  and  shaft  angle  is  used. 
Results 

Optimum  designs  have  been  obtained  with  and  without  the  thrust  constraint  and  the  two 
results  are  compared.  For  ease  of  reference  the  following  notation  will  be  used 
Case  1:  optimum  design  with  thrust  constraint, 
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Case  2:  optimum  design  without  thrust  constraint. 

Case  2  signifies  a  design  where  the  hub  shears  can  be  minimized  at  the  expense  of  reducing  the 
rotor  thrust.  The  remaining  constraints  are  the  same  for  both  the  cases.  Note  that  this  would 
mean  keeping  the  blade  thrust  weighted  solidity  fixed  during  optimization. 

Some  typical  results  obtained  are  presented  in  Figs.  2-4.  Fig.  2  presents  a  comparison  of 
the  4/rev  vertical  and  the  3/rev  inplane  shears,  which  are  the  objective  functions.  The  reductions 
are  larger  in  Case  2  due  to  the  fact  that  the  rotor  has  much  lower  thrust,  however,  the  optimum 
blade  in  Case  1  also  has  lower  values  of  these  shears.  It  was  of  interest  to  investigate  whether 
optimization  causes  any  changes  in  the  blade  section  airloads.  Therefore,  the  normal  and  inplane 
complonents  of  the  total  aerodynamic  forces  on  the  section,  the  lift  L  and  the  drag  D,  resolved 
with  respect  to  the  hub  plane  and  denoted  by  F  and  F  respectively,  are  plotted.  Fig.  3  presents 

Z  X 

the  azimuthal  distribution  of  F^  at  a  radial  station  of  y/R  =  0.99.  The  figure  indicates  that 

optimization  does  not  change  the  magnitude  of  F  in  case  1,  but  reduces  its  amplitude 

z 

significantly  in  Case  2.  Fig.  4  presents  a  comparison  of  the  rotor  power  coefficient,  Cp  , 


normalized  with  respect  to  solidity  a,  for  the  reference  and  the  optimum  blades.  Even  without 
any  consideration  of  the  rotor  power,  i.e.,  the  power  not  included  in  the  optimization  formulation, 
the  optimum  blades  in  case  1  does  not  suffer  any  increase  in  the  total  power  and  in  case  2  the 
power  requirement  reduces  significantly. 

The  final  work  will  contain  more  details  of  the  optimized  blade.  For  example  a  complete 
history  of  the  constraints  (flapping,  lagging,  control  moments,  etc.)  will  be  studied.  More  results 
will  be  included  to  demonstrate  the  effect  of  reducing  the  4/rev  vertical  and  the  3/rev  inplane 
shears  on  the  blade  section  airloads.  Since  a  comparison  of  Cases  1  and  2  indicate  a  trade-off 
between  reductions  in  the  harmonics  of  the  loads  and  maintenance  of  a  certain  thrust  level, 
results  aic  being  obtained  with  several  other  thrust  conditions  to  establish  a  pattern  which  would 


be  of  interest  to  a  blade  designer.  Additional  results  ascertaining  the  effect  of  such  optimization 
on  blade  performance  will  also  be  obtained. 
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Fig.  2.  Comparison  of  objective  functions 


c 

'co 

UD 

c 


D 

o 

sz 


£ 


C\J 

0 

C/D 

0 

O 


0 

E 

i— 

o 


c 

0 

c 

o 

Q. 

c 

l— 

o 

o 


"D 

0 

O  li_n  JD 


0.00403  0.00403 


INTEGRATED  ROTORCRAFT  OPTIMIZATION 


Acoustics  Structures 


FOCUS  OF  CURRENT  WORK 


<Z) 

Eh 

Oh 

H 

O 

£ 

O 

O 

!* 

H 

« 


W) 


00 

•  r-H 

a 

•  pH 

u 

a 

bo 

00 

Jh* 

3 

"d 

•  H 

00 

3 

a 

a 

an 

o> 

bo 

a 

ns 

a 

•  pH 

00 

a 

OO 

"0 

£ 

p3 

© 

a 

> 

<D 

*H 

3 

n3 

0 

a 

b© 

JD 

o 

•  pH 

d 

!> 

CO 

•  pH 

GO 

pJ-H* 

44? 

a 

pa 

4J 

a 

a 

pd 

o 

a 

a 

— • 

«H 

o 

0 

00 

ns 

p3 

44? 

a 

a 

a 

a 

tH 

a 

a 

0 

£ 

y 

•  pH 


0 
•  pH 

4-? 

d 

„  >»  .a 

ft  ^  s 

g  T3  *43 
r9  a  ft 

U  c«  o 


O 

U  g 

©  5 

d  2 


0) 


a 

•ph 

d 

© 

d 

a 

W> 


+5 

a 


a 

o 

•  pH 

+3 

a 

S3 

•  pH 


oo 

jg 

3 

cC 


Q)  •»H  *pH 


&  a 

O  5> 


jg 

IS 

u 

c 

0 

•  pH 

+» 

a 

S3 

•pH 

a 

•  •H 

4-? 

a 

o 

pa 

o 

a 

o> 

+3 

a 

v 

ns 

pin 

3 

o> 

pd 

44? 

a 

•  pH 

• 


I 


TEST  PROBLEM  (” Reference  Blade”) 


OPTIMIZATION  PROBLEM 


0Q 

P 

cd 

03 

A 

OQ 

P 

0 

o 

p 

03 

fl 

cd 


•pH 

> 

03 

P 


CO 


P 

0 

0 

p 


0 
•  pH 

P 

cd 

Pi 


P 

0 

0 

P 

V 

u 

cd 


P 

cd 


Pi 
03 

’  a 

~  cd 
p 


p 

o 

o 

Pi 

V 

13 

cd 


«  ^ 
0  d 


■"o  13 

0 

-  ^ 
^  o 

Q  *rH 

*  rt 
N  P 

N 

PH 

W 


o 

•H 

p 

o 

fl 

a 


«  ii  J 

a  I  I 


•s  O 


cd 

o 

•pH 

p 

v 

> 


03 

> 

•pH 

p 
u 
© 

2s  Hi 


> 

03 

P 


cd 

•»H 

P 

> 


0 Q 

S3 

-  w  *pH 

S  13 

cd 

p 


a  a 


w  • 


W> 

•H 

03 

03 

H 


03 
03 

03  _ 

d  ? 

p 

o 

■p  2 

CO  tJ 


•  pH 


•  • 

OQ 

d 

•pH 

cd 

p 

-p 

OQ 

d 

o 


p 

cd 


to 

d 

•  pH 

to 

to 

cd 

pH 

> 

03 

p 

13 

d 

cd 

P*H 

cd 

d 

o 

•  pH 

03 

P 

0 

-p 

13 

d 


03  cd 

* 


to 

d 


•pH 


OQ 

3  ^ 

?3  a 

s  * 

p  cd 

>  > 

03  03 

P  P 

CO  CO 


to 

jd 

"a 

cd 

cd 

o 

•pH 

P 

03 

Jh 

p 

S3 

a 

p 

OQ 

P 

cd 


OQ 

03 

•pH 


OQ 

P 

d 

03 

a 

a  ^ 


o 

d 

03 

d 

0< 

03 


cd 

•H 

P 

0  OQ 

3  2 

•  PH  g 

l-H  4J 

Cd  OQ  'P 

d  ^ 

*H  H  yj 

P  gl  .pH 

cd  £P  03 

P  S3  !> 

O 

•  *pH 

O  P 

S3  03 


03 

13 

cd 


QQ 

03 

3 

cd 

•pH 

P 

s 

a 

•pH 

OQ 

03 

-8  * 

d 

p 

p 

p 
o 

p 
Q 


d 

o 

OQ 

13 

d 

d 

o 


«!  O  PQ  ti  0 


FORMULATION  OF  MULTIPLE  OBJECTIVE 


subject  to  <  0 


FLOW  CHART  OF  OPTIMIZATION 

PROCEDURE 


analysis  analysis 


Structural  analysis: 


GO 

GO 

© 

U 

-M 

QQ 


GO 

U 

•  1-4 

d 

0 

cd 

fl 

1—4 

1—4 

X> 

cd 

0 

*-4 

-P 

GO 

u 

cd 

•1-4 

s 

X> 

cd 

cd 

fl 

0) 

>> 

-p 

CO 

"d 

0 

0 

& 

z 

• 

• 

■ 

I 


ANALYSIS,  CONTINUED 


TRIM  CRITERION  USED 


RESULTS 


COMPARISON  OF  OBJECTIVE  FUNCTIONS 


0 

O 

C 

0 


0 

DC 


lO 

CO 


o 

CO 


in 

■ 

CM 


© 

c\i 


in 


T- 

CM 

0 

0 

0 

C/> 

CO 

CO 

O 

O 

SW«' 

E 

E 

3 

3 

E 

■  MB 

E 

■  Mi 

Q. 

Q. 

o 

o 

0  M  tmmm 


0  I 


O 

CM 


in 


m 

o 


4/rev  3/rev 

vertical  shear  (fz)  inplane  shear  (fx) 


NONSTRUCTRUAL  WEIGHT  DISTRIBUTIONS 


Nondimensional  radius,  y/R 


CENTRIFUGAL  STRESS  DISTRIBUTIONS 


p 

c\T 

o 

CD 

<0 

<0 

CO 

CO 

O 

o 

s— ^ 

£ 

£ 

3 

3 

E 

IM 

£ 

mmm 

O  CL  Q.  £ 

oc  o  o  <o 


lU 

i  <*> 


CO  CN 
<0  ^ 

<J>  (0 
<*-»  -Q 

co  J-: 


Nondimensional  radius,  y/R 


AZIMUTHAL  DISTRIBUTIONS 
F  NORMAL  COMPONENT  OF  AIRLOAD,  Fz 

y/R  =  0.75,  |A=  0.3 


O 

CO 

CO 


o 

Is-  o> 
04  S 
TJ 


0) 

U) 

o  c 

00  CO 


o  o  o  o  o  o 

lO  ^  CO  CM  r 


ir\ 


O 

O 


AZIMUTHAL  DISTRIBUTIONS 
OF  INPLANE  COMPONENT  OF  AIRLOAD,  F 

y/R  =  0.75,  (i=  0.3 


p 

cm" 

o 

(1) 

(0 

(/) 

CO 

CO 

O 

O 

s>— 

E 

E 

3 

3 

E 

■  MM 

E 

■  MM 

a>  a.  o. 
COO 


<D  IT)  CO  CM 


c 

o  o 

c  c  _ 

03  O  X.Q 
“  Q_LL  =1 

cE 

—  O 
o 


Azimuth  angie,  deg 


COMPARISON  OF  NORMALIZED 
TOTAL  POWER  COEFFICIENT 


e> 

**  ***»*. 

a 

O 


Reference  Optimum  Optimum 

(Case  1 )  (Case  2) 


COMPARISON  OF  NORMALIZED 
INDUCED  POWER  COEFFICIENT 


Reference  Optimum  Optimum 

(Case  1 )  (Case  2) 


COMPARISON  OF  NORMALIZED 
PROFILE  POWER  COEFFICIENT 


■ 


(Casel)  (Case  2) 


VARIATION  OF  LOCAL  ANGLE  OF  ATTACK 

WITH  MACH  NUMBER 
y/R  =  0.99 


0  G 
CO  (0 
03  CO 


o  o 


E  E 

D  =3 

E  E 


x 


CO 

E 


LO^rOOCNr-OT-CMCOrT 


Mach  number 


VARIATION  OF  LOCAL  ANGLE  OF  ATTACK 

WITH  MACH  NUMBER 
y/R  =  0.75 


Mach  number 


VARIATION  OF  LOCAL  POWER  COEFFICIENT 

WITH  MACH  NUMBER 
y/R  =  0.750 


r-  CM 

0  0 
if)  if) 
CO  CO 


o  o 


(0 

OCDO 

o  o  o> 
~  II 

0  o 

N  Ow 

o 


CO 

E 


0 

g 

o 


a 

o 

a 

o 


Mach  number 


CONCLUDING  REMARKS 


13  to 

>>  a 

"0  60 


•  pH 


CD 

13  a> 
cd  X3 


© 

*3  d 

d  *— « 

3  43 
<0  . 
d  *■* 

1  ° 
1  0 

© 
cd 


tS 


•o  ® 

5  | 

2  *2 

¥  8 
43  O 

J3  S. 


> 

© 

© 


CO 

13 


13  | 

.d  © 
t  -SI 

5  CD 
>  © 
^  d 

$  A 

a 

A 

N1  #rH 


CO 

d 

O 

«H 

43 

© 

d 

<2 

© 

> 

•H 

43 

© 

2* 

O 


60 

.S  13 

to  to  d 
d  to  § 

•s  d  *d 

A  ^  43 

Oi  ±  d  _ 

^  ©  o  pfS 

cd  2  2  3 


© 

d 


§  fc  * 


CQ 


d 

d 


g  d  .5  «-< 

&  5  «  d 

I  a  § 

«  13 

5  fl 


M  ■§ 

f!  CD 

©  ^ 
•rH 

© 

d 
© 

d 
a* 
© 

13  H  ,© 

a  o  ^ 

rt  g  +s 

»h  £  a 

>  -2  §. 

2  -2  | 

s  a  § 

n  <i  a 

•  • 

CQ 

d 

•N 

a 

43 

CD 

d 

0 

O 


»> 

CQ 

CD 

© 

© 

43 

CD 


§0 

<2 

•pH 

u 

'S 

© 

© 

€\ 

cd 

43 

u 

© 

d 

»fH 


43 

CQ 

d 

u 

X 

43 

U 

O 

43 

o 

© 


43 

X 

to 

IN 

I 


Global  criteria  approach  used  for  multiple 
objective  functions 


GO 

Tj 

>> 

u 

1C 


o 

rH 

3 

•  rH 

TJ 

0) 

•  pH 

3 

43 

X 

0 

GO 

a 

to 

•  rH 

GO 

OP 

H3 


•  rH 

43 

A 

O 


GO 

O 

00 

3 

u 

I 

43 

§ 

I 

43 

0) 

0) 

TJ 

3 


GO 

s 

o 

GO 

•m 

P 

3 

ft 


0 

O 


"S 

•  rH 

3 

p 

43 

00 

s 

o 

v 

43 

GO 

3 

p 

X 

43 


3  xi 

“  i  ^ 


”8 
•  rH 

3 

P 

43 

GO 

3 

0 

o 

43 

00 

3 

p 

X 

43 


43 

3 

O 

r£ 

43  +3 

•^  *rH 

£  £ 

•  •  •• 

H  M 

0)  0) 

GO  GO 

3  3 

n  n 


0) 

Pi 

d  « 

§  rt 

>  H 

fl)  GO 

W  0) 

CO 

3 

O 


•  rH  ^ 
43 


GO 

3 

O 

•  rH 
+3 

O 

3 

”0 

0) 

P 

3 

U 

cS 

•m 

3 

to 

•  rH 

CO 


0 

•O 

p 

cS 

GO 

Pi 

3 

a? 

,3 

GO 

a? 

3 

j3 

"ft 

3 

•  rH 


>» 

43 

Id 

3 

0) 

ft 

5s 

W) 

•  rH 

<a> 


3 

O 

X 

43 
•  rH 


43 

3 

GO 

Is 

o 


O 

s 

TJ 

3 

3 

GO 

0) 

a 

Pi 

cS 

p 

3 

Q) 

r3 

GO 

Id 

a 

•  rH 

43 

•  rH 

P 

a 

p 

a? 


43 

O 


blade  root  remained  within  bounds 


CONCLUDING  REMARKS 


Q 

H 

£ 

z 

M 

H 

Z 

0 

o 


0) 

a 

a 

r—H 

ft 

a 

#f«*! 

"0 

3 

3 

13 

a 

© 

o 

3 

s 

o 

•  pH 

43 

© 

© 

CO 

©H 

0 

CO 

© 

3 

43 

**-h 

fH 

a< 

a 

◄ 


M 

0) 

05 

a 

O 


-p 

a 

3 

© 

43 

CO 

3 

0 

© 

43 

CO 

3 

Vi 

p3 

+3 

O 

43 

© 

3 

"d 

rH 

© 

CO 

3 

o 

3 

•  pH 

3 

0 

•ph 

+3 

© 

3 

13 

© 

v 

/•* 

w 

jz: 


a 

a 

a 

"•I-* 

43 

a 

o 

3 

•  pH 

13 

© 

13 

J3 

"© 

3 

•ph 


43 

0 

3 


3 

© 

© 

43 

•  pH 
© 
© 

© 

© 

3 

3 


Vi 

<S 

V 

© 

Ph 


CO 

3 

^3 


iH 

© 

CO 

3 

o 


3 

•  pH 


© 

© 

3 

© 

Vi 

,© 

©H 

© 

V 


© 

p3 


43 

CO 

3 


© 

CO 

3 

O 

3 

•  pH 


i 

| 

1 

I 

I 


Application  of  DYNOPT  Optimization  Program  For  Tuning 
Frequencies  of  Helicopter  Airframe  Structures 


T.  Sreekanta  Murthy 

Lockheed  Engineering  and  Sciences  company 
Hampton  VA  23666 


ABSTRACT 

Among  the  various  approaches  which  are  employed  in  attempts  to 
reduce  vibrations  in  helicopters,  the  method  of  airframe  structural 
tuning  has  begun  to  receive  renewed  attention.  In  the  airframe  tuning 
approach,  the  natural  frequencies  of  the  airframe  are  separated  from 
the  frequencies  of  the  rotor-induced  loads  to  avoid  resonances  and  to 
reduce  dynamic  responses  via  appropriate  modification  of  the  airframe 
structure.  This  approach  to  the  design  of  airframes  has  the  potential 
for  producing  inherently  low  vibration  helicopters.  Research  work  on 
this  approach  .is  relatively  new  and  is  contained  primarily  in  refer¬ 
ences  1  through  5.  The  cited  papers  all  address  the  use  of  local  struc¬ 
tural  modifications  aimed  at  tuning  airframe  natural  frequencies.  Sev¬ 
eral  methods  have  been  used  in  those  papers  to  identify  airframe  mem¬ 
bers  for  modification.  Sciarra(l)  used  a  strai;*  energy  method;  Done 
and  his  co-workers  (2)  and  Sobey  (3)  used  the  Vincent  circle  trace  of 
responses  as  obtained  through  variation  of  stiffness  in  airframe  mem¬ 
bers;  Hanson  (4)  did  a  comparative  study  of  the  above  two  methods;  King 
(5)  used  a  modal  tuning  method.  Although  the  term  "structural  optimi¬ 
zation"  is  used,  for  the  most  part  the  papers  do  not  apply  structural 
optimization  in  the  usual  way.  Rather  the  term  is  used  to  indicate  that 
any  local  structural  modifications  which  have  been  made  are  the  best 
based  on  ad  hoc  considerations  such  as  reduction  of  responses  or  strain 
energy  in  a  member.  The  use  of  optimization  techniques  for  airframe 
structural  modifications  is  rather  recent.  Indeed,  references  6  through 
8  describe  what  are  apparently  the  only  published  reports  on  the  appli¬ 
cation  of  such  techniques  to  simple  models  of  airframes. 

As  part  of  the  on-going  j  search  on  helicopter  optimization  for 
vibration  reduction  at  LaRC,  a  study  was  initiated  to  develop  computa¬ 
tional  procedures  for  optimization  of  practical  airframe  structures 
under  dynamic  constraints.  One  of  the  objectives  of  the  study  was  to 
develop  sensitivity  analysis  procedures  for  constraints  on  the  steady- 
state  dynamic  response  displacement  of  airframes  under  rotor-induced 
loads.  Research  work  in  this  regard  involved  development  of  a  solution 
sequence  based  on  direct  matrix  abstraction  program  (DMAP)  of 
MSC/NASTRAN  to  compute  the  sensitivity  coefficients  for  the  dynamic 
response  constraints.  The  sensitivity  results  from  the  application  of 
the  solution  sequence  to  an  elastic  line  model  of  a  helicopter  airframe 
structure  are  discussed  in  reference  (9).  In  the  research  study,  a  com¬ 
putational  procedure  based  on  the  nonlinear  programming  approach  of 
optimization  was  developed  which  incorporates  the  dynamic  response  sen¬ 
sitivity  solution  sequence.  The  procedure  has  the  capability  to  solve 
optimization  problems  with  frequency  and  static  constraints  in  addition 
to  forced  response  constraints.  Implementation  c  the  procedure 
resulted  in  a  computer  code,  designated  DYNOPT,  for  optimization  of 
airframes  under  dynamic  constraints.  This  paper  describes  the  features 
of  the  DYNOPT  code  and  demonstrates  its  application  to  optimization  of 
the  finite  e1*.  ent  model  of  a  Bell  AH-1G  helicopter  airframe.  The 
results  of  tl ‘  .  study  are  described  and  discussed  in  this  paper. 
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ABSTRACT 


Very  little  is  known  about  isolated  rotor  stability  under  dynamically  stalled 
conditions.  This  is  particularly  true  of  lead-lag  modes  that  are  at  best  weakly  damped 
and  for  wmch  current  methods  of  predicting  damping  levels  merit  considerable 
refinement.  Given  the  complexity  of  the  flow  field  of  a  rotor  blade  in  dynamic  stall,  it 
is  virtually  imperative  that  analytical  results  are  compared  with  a  comprehensive 
data  base  of  test  results.  Accordingly,  the  present  study  investigates,  both 
experimentally  and  analytically,  the  effects  of  dynamic  stall  lift  and  dynamic  stall 
drag  on  the  flap-lag  stability  of  an  isolated  hingeless  rotor.  The  emphasis  is  on  the 
correlation  with  tests  performed  on  a  three-bladed  isolated  model  rotor  in  the  U.S. 
Army  Aeroflightdynamics  Directorate’s  7  by  10  foot  (2.1  by  3.0  m)  wind  tunnel.  The 
data  were  collected  in  two  tunnel  entries.  Model  configuration  and  test  procedures 
were  essentially  the  same  for  both  entries.  As  shown  in  Figure  1,  the  correlation 
covers  a  wide  range  of  test  conditions  for  various  values  of  rotor  speed  (fi),  collective 
pitch  ( 0o ),  shaft  angle  (as),  flap-lag  structural  coupling  ratio  (R)  and  advance  ratio 
(/i).  It  includes  cases  that  vary  from  near  zero  thrust  conditions  in  hover  and  forward 
flight  to  highly  stalled  forward  flight  conditions  with  advance  ratios  as  high  as  0.55 
and  shaft-tilt  angles  as  high  as  20°. 

The  experimental  model  closely  represents  a  simple  model  of  a  hingeless 
rotor  with  rigid  blades  having  spring  restrained  fla  g  hinges.  This  intentional 
simplicity  facilitates  isolation  of  dynamic  stall  aspect  the  s'  'ty  problem.  The 
structural  details  of  the  analytical  model  were  chosen  to  correspond  as  closely  as 
possible  to  the  experimental  model.  The  aerodynamic  representation  is  based  on  the 
ONERA  unified  dynamic  stall  lift  and  drag  models.  The  nonlinear  equations  of  blade 
motions  and  stall  dynamics  are  perturbed  about  a  periodic  orbit  or  forced  response, 
and  the  damping  values  (<r,  1/sec)  are  evaluated  from  the  Floquet  eigenanalysis. 


-2- 


The  predictions  are  based  on  four  theories  :  1)  linear  theory,  a  constant 
airfoil  lift-curve  slope  and  a  constant  profile  drag  coefficient,  2)  quasi-steady  stall 
theory,  which  is  the  linear  theory  refined  to  include  quasi-steady  nonlinear  lift  and 
nonlinear  drag  characteristics  of  the  airfoil  section,  3)  dynamic  stall  lift  theory,  with 
quasi-steady  drag  and  4)  complete  dynamic  stall  theory,  which  includes  dynamic  stall 
lift  and  dynamic  stall  drag.  Figure  2  shows  the  correlation  in  hover  for  a  collective 
pitch  angle  0o=  8°(Figure  2a),  and  in  forward  flight  for  0o=  6°  and  os=  20°(Figure 
2b).  Here,  the  dynamic  stall  effects  are  negligible,  as  determined  from  the  stall  angles 
based  on  trim  values  (|  angle  of  attack  aj  <  12°).  Nevertheless,  the  correlation 
facilitates  an  improved  appreciation  of  the  dynamic  stall  effects  on  lag  damping.  The 
predictions  cluster  into  two  categories  of  linear  theory  vis-a-vis  the  other  three 
theories.  In  hover,  as  seen  from  Figure  2a,  the  linear  theory  is  fairly  adequate  in  the 
prediction  of  increasing  and  decreasing  trends  of  damping  level,  including  the 
reduction  in  damping  level  around  400  rpm,  where  the  flap  and  lag  frequencies 
coalesce.  However,  the  discrepancies  between  the  data  and  the  linear  theory 
consistently  increase  for  fi  >  400  rpm.  By  comparison,  the  remaining  three  theories 
provide  excellent  correlation  throughout.  To  isolate  the  effects  of  nonlinear  quasi¬ 
steady  section  drag  characteristics  in  substall,  additional  correlations  (not  shown) 
were  conducted  between  the  data  and  the  linear  theory  refined  to  include  nonlinear 
quasi-steady  section  drag  in  substall.  It  was  seen  that  these  substall  nonlinear  drag 
characteristics  significantly  improve  the  correlation,  as  is  typical  of  low-Reynolds- 
number  effects  of  test  conditions  (in  hover,  Reynolds  number  was  close  to  240,000  at 
the  blade  tip).  In  forward  flight,  as  seen  from  Figure  2b,  the  linear  theory  provides 
poor  correlation  by  predicting  rapidly  increasing  damping  level  with  increasing 
advance  ratio  y,  whereas  the  data  show  nearly  constant  or  slowly  decreasing  damping 
level  with  increasing  y.  By  comparison,  the  other  theories  provide  excellent 
correlation.  As  in  hover,  the  improvement  is  primarily  due  to  nonlinear  profile  drag 


effects. 


The  stalled  high-thrust  conditions  are  presented  in  Figure  3  for  the  collective 
pitch  angles  of  0°  and  3°,  both  for  a  shaft-tilt  angle  of  16°.  The  predictions  of  linear 
and  quasi-steady  stall  theories  contrast,  each  providing  good  correlation  for  one  case 
and  qualitatively  inaccurate  prediction  for  the  other.  For  example,  the  linear  theory 
gives  reasonably  good  correlation  only  for  the  case  of  zero  collective  pitch  angle 
(Figure  3a).  But  for  the  second  case  of  3°  collective  pitch  angle  (Figure  3b),  it 
predicts  increasing  damping  with  increasing  advance  ratio  —  a  qualitatively 
inaccurate  prediction  in  comparison  with  the  test  data.  In  contrast,  quasi-steady  stall 
theory  gives  good  correlation  for  collective  pitch  angle  0o=3°(Figure  3b),  and 
qualitatively  inaccurate  predictions  for  collective  pitch  angle  0o=O°(Figure  3a).  On 
the  other  hand,  the  dynamic  stall  lift  theory  with  quasi-steady  drag  characteristics 
gives  fairly  accurate  prediction  in  both  cases.  In  other  words,  quasi-steady  stall 
theory  refined  to  include  dynamic  stall  lift  consistently  improves  the  correlation. 
Furthermore,  inclusion  of  dynamic  stall  drag  generally  provides  further  quantitative 
improvement,  except  for  minor  degradation  in  the  3°  collective  pitch  angle  case,  as 
seen  from  Figure  3b. 

The  preceding  correlation  covers  a  very  small  section  of  the  correlation 
covering  the  entire  data  base  of  nearly  1200  test  points  shown  in  Figure  1.  The  full 
presentation  with  about  20  figures  is  planned  to  cover  the  following  topics  :  dynamic 
stall  lift  and  dynamic  stall  drag  in  Floquet  stability  analysis,  trim  analysis  and 
perturbation  about  a  periodic  forced  response,  computational  reliability,  test 
procedures,  generation  of  data  base  and  correlation  of  theory  with  data  to  isolate 
dynamic  stall  effects.  Future  research  areas  uncovered  during  the  correlation,  such  as 
possible  refinements  in  the  ONERA  dynamic  stall  lift  and  drag  model  parameters,  are 
identified  as  well. 


4- 


The  detailed  experimental  and  analytical  investigations  lead  to  the  following 

major  findings : 

1)  Conventional  Floquet  analysis  with  the  ONERA  u...  „d  dynamic  stall  lift 
and  drag  models  provides  a  viable  approach  for  predicting  lead-lag  mode 
damping  in  stalled  forward  flight  conditions. 

2)  For  low-advance-ratio  cases  (n  <  0.20),  the  nonlinear  substall  drag 
characteristics  of  the  airfoil  section  significantly  improve  the  correlation  due 
to  low-Reynolds-number  effects  of  the  model  test  data,  as  investigated 
earlier  on  the  basis  of  analytical  and  experimental  results  restricted  to 
hovering  conditions. 

3)  In  comparison  with  the  linear  and  quasi-steady  stall  theories,  the  theory 
with  dynamic  stall  lift  and  quasi-steady  drag  consistently  improves  the 
prediction  for  the  entire  database. 

4)  Dynamic  stall  drag  generally  provides  further  quantitative  improvement. 
Minor  degradation  in  correlation  for  some  isolated  cases  as  seen  in  Figure  3b 
merits  further  investigation. 
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Fig.  2a  :  Lag  Mode  Damping  Correlation  in  Hover  (  6 
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MEASUREMENTS  AND  PREDICTIONS  CONCERNING  ROTOR  STABILITY 


I.  Cafarelli  and  J.J.  Costes 

ONERA  -  BP  72 
92322  CHATILLON  Cedex 
France 


The  necessity  to  better  understand  the  principal 
phenomena  underlying  stability  problems  in  rotors  has 
recently  led  the  ONERA  {Office  National  d'Etudes  et  de 
Recherches  A^rospatiales)  to  build  a  test  rotor  specifically 
designed  for  this  purpose. 

This  rotor  is  simple  enough  to  permit  a  physical 
insight  into  stability  behaviour  as  a  function  of  various 
fundamental  parameters  in  both  hover  and  forward  flight.  It 
is  also  expected  to  yield  a  data  base  for  validating 
existing  and  future  predictive  codes. 

The  first  series  of  tests  was  carried  out  on  an 
isolated  rotor  configuration.  The  blades  were  hinged  in 
flap  and  cantilevered  in  plane.  The  natural  frequency  in 
torsion  of  each  blade  was  varied  so  as  to  create  different 
coupling  conditions  between  •  the  flap  and  torsion  modes. 
Moreover;  the  relative  importance  of  other  parameters  was 
measured  such  as  rotor  speed,  collective  pitch,  rotor  tilt, 
forward  speed  and  63  coupling. 

Three  different  predictive  codes  are  beign  studied 
and  validated  with  respect  to  these  results.  Two  of  these 
are  restricted  to  hover.  They  differ  in  the  modeling  of  the 
rotor  dynamics.  Forward  flight  conditions  are  predicted 
using  a  Floquet  mode  resolution  code. 

First  results  show  the  sensitivity  of  the  modal 
coupling  to  the  accurate  and  detailed  modeling  of  the  rotor 
dynamics . 


Experimental  Investigations  of  Helicopter  Structural  Dynamics 
and  the  Interaction  with  Vibration  Reduction  Systems 


Manfred  Degener 


DLR  -  German  Aerospace  Research  Establishment 
Institute  for  Acroclasticity 

Bunsenstrassc  10,  D-3400  Goettingen,  Fed.  Rep.  Germany 


Abstract: 

Vibrations  are  one  of  the  major  problems  concerning  the  performance  and  ride  com¬ 
fort  of  helicopters.  The  vibration  behaviour  depends  both  on  the  rotor  induced  oscil¬ 
lating  loads  and  on  the  structural  dynamic  characteristics  of  the  helicopter  itself. 
Therefore,  vibration  treatment  must  take  into  account  both  aspects. 

The  design  of  modern  helicopters  is  a  combination  of  a  light-weight  airframe  struc¬ 
ture  and  large  concentrated  masses  such  as  gearbox  and  engines.  This  results  in 
complex  three-dimensional  vibration  behaviour  with  many  normal  modes  lying  in'  the 
same  frequency  range  as  the  oscillating  rotor  loads,  which  requires  a  thorough  ana¬ 
lytical  and  experimental  investigation.  In  Germany,  MBB  and  DLR  have  a  long 
cooperation  in  the  field  of  helicopter  structural  dynamics.  Ground  vibration  tests 
have  been  performed  by  DLR  on  the  complete  MBB  helicopter  family. 

This  paper  gives  a  survey  of  experimental  investigations  of  helicopter  structural  dy¬ 
namics.  It  describes  the  features  of  the  DLR  ground  vibration  lest  facility  with  special 
emphasis  on  applications,  experiences,  and  improvements  in  testing  helicopter  struc¬ 
tures.  This  relates  to  realistic  boundary  conditions  during  testing,  sensor  instrumen¬ 
tation,  dynamic  excitation,  and  data  reduction.  An  important  topic  is  the  correlation 
of  experimental  results  with  the  mathematical  model.  Concepts  for  correlation,  such 
as  orthogonality  check  and  modal  assurance  criterion,  arc  discussed. 

More  stringent  vibration  requirements  have  caused  research  organisations  and  the 
helicopter  industry  to  develop  special  vibration  reduction  systems.  The  effectiveness 
of  these  systems  can  be  strongly  influenced  by  interaction  with  the  structural  dyna¬ 
mics  of  the  helicopter.  This  requires  structural  tailoring  of  the  airframe  to  minimize 
vibration.  Experimental  results  of  a  recent  ground  vibration  test  on  a  helicopter  with 
and  without  a  vibration  reduction  system  arc  reviewed  to  illustrate  this  effect. 
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Frequency  Characteristics  of  Rotor  Excitation 


Aspects  of  Helicopter  Structural  Dynamics 

Design  Features: 

-  leight-weight  airframe 

-  large  concentrated  masses  (gearbox,  engines) 

-  combination  metallic/composite  materials 

-  vibration  reduction  systems 

-  different  payload  configurations 

Vibration  Characteristics: 

-  complex  three-dimensional  vibration  modes 

-  high  modal  density 

-  many  local  modes 

-  nonlinearities 

Requirements  for  Vibration  Analyses: 

-  large  finite  elemente  models 

-  sophisticated  test  methods 


Phase  Resonance  Method 


Equations  of  Motion  (physical  coordinates): 

m  u(P.t)  +  d  u(P,t)  +  k  u(P,t)  =  f(t) 

Modal  Transformation: 

n 

u(P,t)  =  £<D,(P)  qj(t)  =  <I>(P)  q(t) 

i=0 

Equations  of  Motion  (generalized  coordinates): 

M  q(t)  +  D  q(t)  +  K  q(t)  =  d>T  f(t) 

with 

M  =  <J>r  m  :  generalized  mass  matrix  (diagonal) 

K  =  <J>Tk<I>  :  generalized  stiffness  matrix  (diagonal) 
D  =  <I>Td<I>  :  generalized  damping  matrix 


Helicopter  Ground  Vibration  Test  Set-up 


Sensors  (Accelerometers) 


DLR  Ground  Vibration  Test  Facility 
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BO  108  -  Measurement  Point  Plan 


DLR 


Measurement  of  Normal  Modes 

Identification  of  Normal  Modes 

-  Measurement  of  MIF's  for  various  exciter  configurations 

Isolation  of  Normal  Modes 

-  optimizing  MIF-value  (  >  0.9  ) 

-  frequency  adjustment 

-  variation  of  exciter  forces  and  locations 

-  on-line  graphic  display 

Measurement  of  Modal  Parameters 

-  normal  frequency 

-  normal  mode  shape 

-  generalized  mass 

-  damping 

-  linearity  check 


Identification  of  Normal  Modes  by  MIF's 


Orthogonality  Check 


Generalized  Mass  Matrix: 

M  =  <I>Tm  d) 


Normalized  to  Mm  =  1: 


with 


M  =  <I)Tm  cl) 
0  <  Mik  <  1 


0<  MAC(i,k)  <  1 


Orthogonality  Check 


Modal  Assurance  Criter 


BK  117:  FEM  mode  shapes 


DLR  Ground  Vibration  Tests  on  MBB  Helicopters 


C\J 


DLR 


Antiresonance  Vibration  Isolation  System  (ARIS)  [3] 


Hydraulic  System  Mechanical  System 


Antiresonance  Vibration  Isolation  System  (ARIS) 
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Hydraulic  Force  Isolator  (MBB  Design)  [3] 


Frequency  /  Hz 


Vertical  Vibrations  at  Pilot's  Seat  (Fz  =  400  N) 

without  ARIS,  with  Rotor  Excitation  Device  (ca.  200  I 
with  ARIS,  with  Rotor  Excitation  Device 
with  ARIS,  without  Rotor  Excitation  Device 


Mode  Indicator  Function 


BO  108  -  Shake  Test  Results 


25  28  27  28  29  30 

Frequency  /  Hz 


Mode  Indicator  Function  (Vertical  Rotor  Excitation) 

without  ARIS,  with  Rotor  Excitation  Device  (ca.  200  kg) 
with  ARIS,  with  Rotor  Excitation  Device 
with  ARIS,  without  Rotor  Excitation  Device 


BO  108  -  Shake  Test  Results 
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Frequency  ✓  Hz 


Nonlinearity  investigations 
-  (Fz  =  560  N) 


—  ■  (Fz=  1200  N) 

—  (Fz  =  2000  N) 


Conclusions 


Helicopter  vibrations  are  influenced  by  aerodynamic  rotor  ex¬ 
citation  as  well  as  structural  dynamic  properties 


Advanced  analytical  and  experimental  methods  are  required  to 
evaluate  complex  helicopter  dynamics 


Interaction  of  structural  dynamics  with  vibration  reduction  sy¬ 
stems  requires  airframe  structural  tailoring  to  minimize  vibra¬ 
tion  levels  ’ 


Vibration  reduction  systems  can  reduce  vibration  levels  signifi¬ 
cantly  as  confirmed  by  ground  and  flight  tests. 
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EXPERIMENTAL  IDENTIFICATION  OF  HELICOPTER  ROTOR  DYNAMICS 

USING  KINEMATIC  OBSERVERS: 

An  Abstract 

Robert  McKillip,  Jr. 

Mechanical  and  Aerospace  Engineering  Dept. 

Princeton  University,  Princeton,  NJ  08544 


Recent  investigations  using  blade-mounted  instrumentation  has  revealed  the 
potential  for  identification  of  rotor  dynamics  as  a  two  step  process  [1].  Estimation  of  rotor 
modal  state  variables  is  performed  using  knowledge  of  the  kinematics  of  the  sensor  signal 
content,  and  then  these  estimated  state  time  histories  are  used  to  identify  the  dynamic 
equation  coefficients  for  a  particular  linearized  mathematical  model.  Such  an  approach  has 
the  advantage  that  it  decouples  the  assumed  rotor  dynamic  model  structure  from  the  state 
estimation  phase,  allowing  a  wide  variety  of  assumed  model  forms  to  be  identified  from  the 
same  set  of  state  estimates.  Additional  uses  for  the  separated  state  and  parameter 
estimation  scheme  include  adaptive  rotor  control  appli cations  as  well  as  on-line  stability 
analysis.  Fundamental  to  the  technique  is  the  use  of  blade-mounted  accelerometers  to 
provide  modal  acceleration  information  to  the  state  observer,  taking  the  place  of  an 
assumed  dynamic  model  as  is  normally  used  for  more  conventional  instrumentation  sets. 

A  series  of  hover  and  forward-flight  model  rotor  experiments  will  be  described  that 
used  blade  accelerometer  data  to  reconstruct  both  rotor  modal  state  variables  and  linear 
equation  parameters.  The  model  tests  were  conducted  at  Princeton’s  Rotorcraft  Dynamics 
Laboratory  using  a  four-foot  diameter  rotor  having  the  capability  for  both  Individual 
Blade  Control  and  Higher  Harmonic  Control  [2j.  Tne  data  processing  technique  will  be 
detailed,  and  results  of  the  parameter  identification  analysis  will  be  compared  to  linearized 
theory  (see  Fig.  1).  Examination  will  be  made  of  the  identified  model’s  ability  to  predict 
closed-loop  blade  response  experimental  data  for  both  hover  and  forward-flight  regimes. 
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Parameter  Estimation  Using 
Low-Pass  Filtered  Accelerometers 


Fig.  la:  Signal  Processing  Flow  Diagram 


Fig.  lb:  Flap  Mode  State  Identification  Results 


Fig.  lc:  Flap  Equation  Parameter  Identification  Convergence  History 


Hover  Test  Results 


ABSTRACT 


VALIDATION  OF  ROTOR  VIBRATORY  AIRLOADS  AND 
APPLICATIONS  TO  HELICOPTER  RESPONSE 

Jing  G.  Yen 

Bell  Helicopter  Textron,  Inc. 


The  prediction  of  rotor  vibratory  airloads  and  structural  loads  has 
been  a  challenge  to  rotorcraft  engineers.  The  development  of  advanced 
theories  and  computer  codes  requires  extensive  correlations  with 
experimental  data.  The  validation  is  to  demonstrate  that  the  theories  do 
in  fact  correctly  represent  the  physics  of  the  problem  and  that  the 
computer  code  is  adequate  to  be  used  as  a  design  tool. 

Bell  Helicopter  Textron,  Inc.,  has  recently  completed,  as  part  of  a 
program  to  improve  rotor  wake  representations,  a  major  cooperative  program 
with  the  Army/NASA  researchers  at  Langley  Research  Center.  The  program 
consisted  of  the  testing  of  a  scaled  model  with  pressure  instrumented 
blades.  The  model  was  a  0.2  Mach  scale  of  Bell  Advanced  Light  Rotor  (ALR). 
The  ALR  was  a  four-bladed  bearingless  rotor  with  advanced  airfoils. 
Pressure  transducers  were  installed  at  five  blade  spanwise  locations.  The 
testing  was  conducted  at  the  Langley  VSTOL  wind  tunnel  in  June-July,  1989. 

The  experimental  vibratory  airloads  (the  normal  force)  are  obtained  by 
integrating  the  individual  pressures.  The  airloads  data  are  presented  in 
time  histories  as  well  as  harmonic  amplitude  and  phase  angles.  The  data 
are  depicted  as  functions  of  the  advance  ratio  and  the  rotor  thrust 
coefficient.  Comparisons  of  airloads  are  made  between  the  experimental 
data  and  theoretical  values.  The  theoretical  values  are  calculated  using 
three  rotor  wake  models:  Drees  wake,  CAMRAD/Scul ly  free  wake,  and 
CAMRAD/ Johnson  free  wake.  Computed  blade  and  hub  loads  for  the  ALR  using 
the  free  wake  models  are  correlated  with  those  measured  in  flight. 

Effects  of  the  advance  ratio  and  the  rotor  thrust  coefficient  on  the 
vibratory  airloads  will  be  discussed.  Similarities  of  the  ALR  with  those 
reported  in  the  literature  will  also  be  addressed.  Validation  of  various 
rotor  wake  models  in  calculating  the  vibratory  airloads  and  the  rotor 
structural  response  will  be  assessed  based  on  the  ALR  experimental  data. 
Finally,  conclusions  from  this  study  and  recommendations  for  future  work 
will  be  made. 
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